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1.0 INTRODUCTION 


1 . 1 Background 

Interference from the presence of wind tunnel walls in 
model test data has been studied for many years, and correc- 
tions have been devised for the influence of solid and straight 
walls in works by Glauert (Reference 1) and Theodorsen 
(Reference 2). More recently the idea of adaptable walls that 
flex to conform to the wall-free stream lines about the model 
has been pursued actively both in the USA and in Europe 
(References 3-7). See also the bibliography prepared by 
Tuttle and Mineck (Reference 8). The general flexing of the 
walls of a wind tunnel introduces a complex mechanical problem 
and efforts have been made to find simple but still effective 
ways to reduce residual wall interferences to negligible 
values. An innovative scheme using solid side walls and 
flexing rods for upper and lower walls was described by Harney 
(Reference 9) at the Air Force Wright Aeronautical Laboratory. 
This report discusses in some detail the problem of proper 
choice of contouring for hybrid adaptation. At the NASA Ames 
Research Center, Shairer and Mendosa (Reference 10) describe 
research in which controlled air flow through porous walls is 
used in lieu of flexible walls to create the proper outer 
streamline shape. In England, at the University of 
Southhampton, a technique was described (Reference 11) in which 
only the upper and Lower walls of the tunnel were made 
flexible, resulting in a "hybrid adaptation." In this work the 
influence coefficients of each of the positioning jacks was to 
be determined experimentally and these data input to a computer 
for estimation of best wall positioning. At the Arnold 
Engineering Development Center of the U.S. Air Force, a 
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segmented variable porosity scheme has been employed using 
sixty-four individually controlled segments (Reference 12). 

Smith (Reference 13) has investigated the case of wind 
tunnels with solid side walls and 2-D flexing upper and lower 
walls. He presents a method for shaping walls that sharply 
reduces the model centerline upwash interference and also 
reduces the axial gradients. 

The present report also deals with solid wall tunnels 
having only upper and lower walls flexing. An algorithm for 
selecting the wall contours for both two and three dimensional 
wall flexure is presented and numerical experiments are used to 
validate its applicability to the general test case of three 
dimensional lifting aircraft models in rectangular cross- 
section wind tunnels. The method requires an initial approxi- 
mate representation of the model flow field at a given lift 
with walls absent. This representation should be at best a 
solution of the non-linear transonic flow equations to allow 
use of the method up to Mach numbers where wall speeds approach 
unity. The numerical methods utilized are derived by use of 
so-called Green's source solutions obtained using the method of 
images. First order linearized flow theory is employed with 
Prandtl-Glauert compressibility transformations. In addition 
to two dimensional flexing of the upper and lower walls, three 
dimensional flexing is also considered for cases in which three 
rows of jacks could be used or in which only a single center- 
line row of jacks would be used with fixed unclamped side 
edges. Equations are derived for the flexed shape of a simple 
constant thickness plate wall under the influence of a finite 
number of jacks in an axial row along the plate centerline. As 
a final task, the Green's source methods are developed to 
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provide estimations of residual flow distortion (interferences) 
with measured wall pressures and wall flow inclinations as 
inputs . 

1 . 2 Scope of this Report 

This report presents the following Fortran codes developed 
using a VAX /VMS operating system and essentially incorporates 
the Fortran 77 notations and conveniences: 

(1) PHIXZM - Green’s source representation of flow 
distortion at the model location in the wind tunnel produced by 
a model flow field in combination with wall boundary shapes. 
Flow field input can be from numerical computations or from 
computations based on measured tunnel wall values of pressure 
and slope. 

(2) AFMODL - An approximate representation of an AEDC 
model, see Reference 13, page 109, using pointed horseshoe 
vortices for lift of sweptback wings, doublets on the fuselage 
axis to represent fuselage lift, swept source and sink lines on 
the wings to represent thickness and sources on the axis for 
fuselage and balance sting thickness. It also applies the 
Tracor blockage algorithm and provides an input file for PHIXZM 
to estimate residual flow distortion at the model due to the 
hybrid (incomplete) adaptation. 

(3) NONLVN - Applies Tracor blockage algorithm using 
input data from a non-linear code and provides an output for 
introduction in PHIXZM. 
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(4) VEEXPHINO - Computes the wall-free normal velocity 
field at the walls from measured wall slope and pressure. 

(Wall Slope may result from flexure or boundary layer growth). 
Provides input to PHIXZM for estimating residual flow 
distortion at the model region. 

(5) JACK_DISPL - Computes the jack displacements and 
resulting residual normal velocities at the wall for 3-D 
flexible plate upper and lower walls controlled by a single 
control row of jacks. The Tracor Algorithm is applied for 
input into PHIXZM. 

A modified version of AFMODL called AFMODLJ is appended to 
provide input of wall-free normal velocities at the panels and 
at jack stations for the 0.3 Meter Tunnel. 
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2.0 ESTIMATION OF TUNNEL FLOW DISTORTION AT MODEL 
2 . 1 Tracor Algorithm for Reduction of Blockage 

The presence of fixed straight walls in a rectangular 
tunnel can be represented in inviscid flow by a doubly infinite 
set of images of the model as described clearly in Reference 1. 
This is general and true even for transonic and supersonic 
flow. However, when the model image fields produce sizable 
pressure gradients and flow angularity at the model position it 
will clearly introduce errors in drag, moments, and lift forces 
that require correction. The ideal case represented by 
flexible walls contoured to the free air streamlines would 
permit testing with no corrections for walls and having only 
Reynolds number mismatches to be corrected for. The boundary 
layer growth on the walls must of course be included in 
correctly positioning the walls. Harney (Reference 9) and 
Wolf, et. al (Reference 11) have shown the possibility of 
minimizing tunnel flow field distortion by flexing only the 
upper and lower walls leaving the side walls straight and 
thereby reducing mechanical complexity and cost of the wind 
tunnel installation. In both references cited, however, no 
clear scheme for positioning the flexed 2-D walls was 
demonstrated for the general case of an aircraft model having 
both thickness and lift. The Tracor algorithm to be described 
does in fact show remarkable ability to negate the axial 
pressure gradient and the upwash on the centerline of the 
tunnel and to reduce markedly the spanwise upwash variation 
(washin) normally associated with flat side walls. The 
algorithm concept is simple. Consider Figure 1 showing an 
axial view of a lifting model in the tunnel with the streamline 
traces projected in a plane normal to the axis of the tunnel. 
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The sketch on the left shows the wall free conditions and that 
on the right the fixed wall condition. In the lower half of 
the tunnel there is for the case shown a lack of freedom for 
the flow to expand beyond the walls hence a blockage exists 
which requires the flow to speed up abnormally in that section 
to satisfy continuity of mass flow in the tunnel. At the same 
time there is a crowding of upward flowing streamline traces 
near the wing tip leading to an induced upwash over the wing. 

In Figure 2 where the walls are permitted to effectively move 
outward as one moves axially the streamline traces can resume 
their more nearly free air patterns near the model as 
illustrated in the sketch. A logical manner to determine the 
wall slope thus appears to be as follows; allow the wall in 
each half of the tunnel (upper and lower) to permit a net out- 
flow through the boundary equal to the integral of the free-air 
outflow in the half tunnel. This is essentially the concept of 
the Tracor Algorithm treated in this report. It requires, of 
course, a free air estimation of the normal velocity components 
at the wall control surfaces. It also appears logical that a 
non-linear transonicf ree-air solution at a given lift 
coefficient can be effectively used so long as the disturbance 
velocities at the wall are truly small. It is also clear that 
the lift coefficient of the computed data and the lift 
coefficient measured or set experimentally should be matched as 
the lift is the primary driver of the outer flow field 
independent to first order of its distribution. As will be 
shown, application of the Tracor algorithm for setting the flex 
walls reduces the axial pressure gradients to negligible 
levels. It also sharply reduces the upwash at the model 
centerline and reduces the lateral upwash gradients to 
generally acceptable levels. 
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2 • 2 Derivation of PHIXZM - A Green's Source Code for 
Estimation of Residual Tunnel Flow Distortion 


Once the residual normal velocities at the tunnel walls 
(control surface) are given from a free-air calculation and 
application of the Tracor algorithm, input into PHIXZM yields 
the tunnel x- and z-wise disturbance velocities near the 
model. PHIXZM utilizes a Green's source concept; that is, a 
normal wall velocity can be cancelled over a small panel by a 
solution of the governing fluid differential equations that 
produces a uniform normal velocity over the panel and zero 
normal velocity at all other points on the tunnel walls. The 
derivation of this Green's source solution is done using the 
method of images for a tunnel of rectangular cross section and 
infinitely long in the axial direction. The image system for a 
Green's source on a wall panel is illustrated in Figure 3. A 
similar figure of paired sources can be drawn for sources on 
the floor. Note that symmetry about the centerline is assumed 
throughout this report. The paneling system used is as 
follows: an even number of panels is taken vertically on the 

wall, the wall panels are square but the floor panels can be 
slightly rectangular depending on the tunnel width to height 
ratio. The program is written to make them as square as 
possible. Distances to panel centers in the X direction (down- 
stream) are governed by an integer index, IX, upstream values 

of IX are negative. Panel coordinates on the walls at a given 

IX location are made dependent on an integer index, IZ, 
starting at 1 just to the side of the centerline on the tunnel 
floor. Thus, panel center coordinates can be specified by two 

indices IX and IZ. The coordinate system and numbering system 

are shown in Figure 4. Positions of image sources are defined 
by indices IZ, IX, and I and L, I being the number of image 


- 8 - 


pairs away from the tunnel vertically and L the corresponding 
index laterally. At a value of I and L equal to M (arbitrary 
integer) the calculations are hastened by replacing the outer 
discrete sources by smearing them uniformly in strength to 
infinity on the plane being used. This permits the effect of 
all the additional sources to be integrated in closed form. If 
the value of M selected is sufficiently large the variation of 
disturbance velocities across the tunnel by the outer sources 
is small enough so that both axial and upwash velocities need 
only be computed at the tunnel center, y & z equal to zero. 


The equations for the axial and vertical velocities at the 
centerline for any sources are as follows: 


<J> (upwash velocity )= 1 5 - 7 *- } 

Z 4n8 [ (xj-x) 2 /S 2 +(y i-y ) 2 +( z j-z) 2 ] 3//2 ' 

( 1 ) 


Z 1 - z 


-Or X 1 “ X 

<j> ( axial velocity)* — { zr-p? 

4rr 6 3 [ ( x j -x ) 2 /B 2 + ( y j-y ) 2 +( z !~z ) 2 ] ' 


( 2 ) 


where 0 is the source strength, 6 is \1 - (Mach) 2 , Xj, yj z lf 
and x, y, z are the coordinates of the source and the field 
point, respectively. These disturbance velocities are derived 
from the potential function, $, the unit source solution to the 
linear first order compressible differential equations of 
motion. The source strength Q, is related to the normal 
velocity, VN, at its panel by the relation 


0 


2» VN • Ax* Az or 2 • VN • AX • Ay 


(3) 
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for wall or floor panel locations. AXAZ is the area of panel. 
Making use of the above considerations the program PHIXZM was 
written to provide $ x and <fr z at points at and near the test 
model with input normal velocity distributions obtained either 
from estimates (linear or nonlinear) for the model under test 
at a given lift coefficient, or from the estimated wall-free 
values calculated by the program VEEXPHINO using measured test 
values of axial and normal velocities at the wall control 
surfaces. Fortran listing of the program PHIXZM is given in 
Appendix A. 

2 . 3 A Simplified ABDC Model Flow-Field Code 

In order to exercise the PHIXZM code and assess the 
beneficial effect of flexwalls set according to the Tracor 
Algorithm, a wall-free flow field computation was needed. The 
following describes a simplified modeling approach applied to 
the AEDC Wind Tunnel Model described in Reference 12. The 
model itself is simple consisting only of a circular body, with 
swept back wings and horizontal tail. The body thickness was 
represented by a single source located behind the nose and a 
sink located at the discontinuous base sting intersection. The 
body lift was represented by two semi-infinite doublet lines 
one originating at the body source and the second (negative) at 
the body sink. Wing and tail thickness were represented by 
swept source and sink lines lying respectively at the leading 
and trailing edges of the airfoils. Wing and tail lift were 
represented by a swept line vortex, in effect a pointed horse- 
shoe vortex. For a given total lift, the angle of attack was 
estimated from simple swept wing theory and the main wing was 
estimated to contribute 80% of the total lift, the tail 20%. 

The model was assumed to be mounted so that it moves upward 
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with angle of attack on an arc about the center of rotation 
located well behind the model. The equations for the flow are 
given in a Fortran listing of the code named AFMODL presented 
in Appendix B. This program also applies the Tracor Algorithm 
and produces a file of the displacements of upper and lower 
surfaces for the Langley 0.3-raeter TCT with flexwalls for the 
AEDC model at arbitrary lift coefficient, CL, and Mach number. 
This file is called ZDISPL.DAT. The program also produces a 
file called PHINWALL.DAT containing the estimated wall-free 
normal velocities at the center of each wall panel for input 
into PHIXZM . 

2.4 Extension to Transonic Nonlinear Wall-Free Codes 


The use of small disturbance-linearized theory limits the 
applicability of the method to moderate combinations of Mach 
number, model/tunnel size ratio, and lift coefficient. Some 
extension can be obtained by using a nonlinear transonic code 
for computing the wall-free flow field up to the condition for 
which wall speeds approach the speed of sound or deviate 
substantially from the main flow Mach number. The code used 
should provide a file called VNZERONL.DAT that can be 
introduced into a modified version of AFMODL to replace the 
simplified model representation calculation of the wall-free 
normal velocities. This has been done and is presented in 
Appendix c as NONLVN . This program generates a file 
PHINWALL.DAT for input to PHIXM and a file called ZDISPL.DAT 
giving the coordinates of the upper and lower walls as dictated 
by the Tracor blockage algorithm. 
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2 . 5 Comparison of Green's Source Method with Closed Form 

Solutions to Assess Numerical Accuracy 

The PHIXZM code for estimating wall induced flow 
distortion at the model location was checked for accuracy by 
comparison with Glauert's (Reference 1) nearly closed form 
equation for the upwash produced by a horseshoe vortex in a 
square tunnel, namely: 

Upwash at the tunnel center = 0.137 SCl/D 2 (4) 

Here S is the wing area, C L the lift coefficient and D 2 the 
tunnel cross sectional area. For the comparison, the program 
AFMODL was modified to produce a horseshoe vortex by setting 
sweep, fuselage diameter, wing and tail thickness and tail 
lift, all to zero. The wall-free normal velocities obtained 
were then input to PHIXZM to obtain the tunnel center upwash. 
For a ratio of S/D 2 of 0.021615 and Cl = 0.5, the Glauert 
value was 0.0014806 and the Green's Source Method 0.0014799. A 
total of 1640 double panels was used for this computation. To 
check the accuracy of the axial velocity computations, a direct 
calculation was made of the axial velocity distribution in the 
tunnel center produced by a unit strength source at the 
origin. This numerical solution was calculated using the method 
of images in a square tunnel. For the comparison, the AFMODL 
program was modified to produce a unit source by setting lift 
coefficient, wing and tail thicknesses and rear end sink 
strength all to zero. Once again the wall-free normal 
velocities were input to PHIXZM to obtain the axial disturbance 
velocity ratios due to the walls. The following table shows 
the result: 
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Table 1 

Axial Disturbance Velocity Ratios 


X/D 

PHIXZM 

DIRECT SOURCE 

0.1 

0.07137 

0.07129 

0.2 

0.13928 

0.13912 

0.3 

0.20095 

0.20072 

0.4 

0.25470 

0.25441 

0.5 

0.29999 

0.29964 


Again 1640 double panels were used. Clearly the numerical 
accuracy is more than adequate considering the general 
approximations in the basic theory. 
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3.0 APPLICATION OF THE CODES TO THE AEDC MODEL IN THE LANGLEY 

0.3-METER TCT FLEXWALL TUNNEL 

3 . 1 Basic Model Results 

The AFMODL program was used to compute the wall-free 
normal velocities for the AEDC model at a lift coefficient of 
0.55 and Mach number of 0.77. The input of these data into 
PHIXZM after modification by the Tracor blockage algorithm 
provides an estimate of the residual flow field distortion 
given in the following table computed using 3240 panel pairs: 

Table 2 


NX= 40 NZ= 20 MACH=0 .770 A=0.330 D*0.330 M= 20 



PHIXM 

PHIZM0 

PHIZMl 

PHIZM2 

-10 

-0 . 64892E-04 

0 . 22807E-04 

0 . 13517E-04 

-0 . 28321E-04 

-9 

-0 . 55712E-04 

0 . 12977E-04 

0 . 15917E-05 

-0.49112E-04 

-8 

-0 .40475E-04 

-0.77779E-06 

-0 . 14702E-04 

-0 . 74822E-04 

-7 

-0 . 19946E-04 

-0 . 19454E-04 

-0 . 36323E-04 

-0 . 10582E-03 

-6 

0 . 42184E-05 

-0 .44735E-04 

-0 . 64688E-04 

-0 . 14249E-03 

-5 

0 . 29777E-04 

-0.79225E-04 

-0.10174E-03 

-0 . 18482E-03 

-4 

0 . 54287E-04 

-0.12646E-03 

-0 . 14982E-03 

-0 .23150E-03 

-3 

0 . 75375E-04 

-0 . 19052E-03 

-0.21123E-03 

-0 . 27874E-03 

-2 

0 . 90839E-04 

-0 . 27524E-03 

-0 . 28765E-03 

-0 . 31974E-03 

-1 

0 .98676E-04 

-0 . 38297E-03 

-0 . 37951E-03 

-0 . 34630E-03 

0 

0 . 97263E-04 

-0 . 51315E-03 

-0.48547E-03 

-0 . 35270E-03 

1 

0.85857E-04 

-0.66133E-03 

-0.60220E-03 

-0 . 33981E-03 

2 

0.65291E-04 

-0 . 81894E-03 

-0 . 72458E-03 

-0 . 31709E-03 

3 

0 .38460E-04 

-0 . 97462E-03 

-0 .84623E-03 

-0 . 301 58E-03 

4 

0. 99236 E- 05 

-0.11170E-02 

-0.96049E-03 

-0 . 30923E-03 

5 

-0.15650E-04 

-0 .12378E-02 

-0 . 10617E-02 

-0 . 34 333E-03 

6 

-0 . 35440E-04 

-0 . 13340E-02 

-0 . 11464E-02 

-0 . 39444E-03 

7 

-0 . 49124E-04 

-0 . 14070E-02 

-0 . 12140E-02 

-0 .44952E-03 

8 

-0 . 57658E-04 

-0 . 14609E-02 

-0 . 12660E-02 

-0 . 49944E-03 

9 

-0 .61693E-04 

-0 . 15005E-02 

-0 . 1 3055E-02 

-0 . 54 104E-03 

10 

-0 . 61263E-04 

-0 . 15304E-02 

-0 . 13358E-02 

-0 . 57531E-03 
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Here IX values of ±10 represent distances of one half the 
tunnel height. They lie just ahead of the fuselage nose and 
just behind the fuselage base. It can be seen that the axial 
values and gradients of the disturbance pressure are negligible 
and the upwash angles are small, the largest being near the 
base of approximately 0.1 degrees. 


The values that would be incurred without flexing the 
walls are presented in the following table: 


Table 3 



NX= 40 NZ= 

20 MACH*0 , 770 

A=0 .330 D=0 . 330 

M= 20 


PHIXM 

PHIZM0 

PHIZMl 

PHIZM2 

-10 

0.20055E-02 

-0 . 32051E-03 

-0 . 32980E-03 

-0 . 37164E-03 

-9 

0 . 33860E-02 

-0 . 16585E-03 

-0 . 17724E-03 

-0 . 22795E-03 

-8 

0 . 489 50E-02 

0 . 68470E-04 

0 . 54544E-04 

-0 . 55792E-05 

-7 

0 . 65023E-02 

0 . 39917E-03 

0 . 38230E-03 

0 . 31280E-03 

-6 

0 . 81652E-02 

0 . 840 16E-03 

0.82021E-03 

0.74240E-03 

-5 

0.98299E-02 

0 . 14003E-02 

0 . 13778E-02 

0.12947E-02 

-4 

0 . 11435E-01 

0 . 20817E-02 

0 . 20583E-02 

0 . 19766E-02 

-3 

0 . 12915E-01 

0 . 28776E-02 

0 . 28569E-02 

0 . 27894E-02 

-2 

0 . 14204E-01 

0 . 37719E-02 

0.37595E-02 

0 . 37274E-02 

-1 

0 . 1 5244E-01 

0.47391E-02 

0 . 47426E-02 

0.47758E-02 

0 

0 . 15989E-01 

0 . 57455E-02 

0 . 57731E-02 

0 . 59059E-02 

1 

0 . 164 13E-01 

0 . 67526E-02 

0 . 68 1 17E-02 

0 . 70741E-02 

2 

0 . 16517E-01 

0 . 77220E-02 

0.78164E-02 

0.82239E-02 

3 

0 . 16328E-01 

0 . 86205E-02 

0.87489E-02 

0.92935E-02 

4 

0 . 15891E-01 

0.94239E-02 

0 . 95804E-02 

0 . 10232E-01 

5 

0 . 15263E-01 

0.10119E-01 

0 . 10295E-01 

0.11013E-01 

6 

0 . 14500E-01 

0 . 10703E-01 

0 . 10890E-01 

0 . 11642E-01 

7 

0 . 13656E-01 

0 . 11178E-01 

0 . 11372E-01 

0 . 12136E-01 

8 

0 . 12775E-01 

0 . 11554E-01 

0 . 11749E-01 

0 . 12516E-01 

9 

0 . 11897E-01 

0.11840E-01 

0 . 12035E-01 

0 . 12799E-01 

10 

0 . 11054E-01 

0 . 12044E-01 

0 . 12239E-01 

0 . 12999E-01 


These disturbances would clearly alter the flow over the model 
to an extent that would make test data only marginally 
correctable . 
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The effect of using a finer panel grid can be seen by 
comparison of the first table of results with the following 
table computed with a total of 7260 panel pairs: 


Table 4 



NX= 60 NZ= 

30 MACH=0 .770 

A=0 .330 D=0 .330 

M= 20 


PHIXM 

PHIZM0 

PHIZM1 

PHIZM2 

-15 

-0 . 41970E-04 

-0.98716E-04 

-0 . 10551E-03 

-0 . 13814E-03 

-14 

-0 . 37653E-04 

-0 . 11470E-03 

-0 . 12213E-03 

-0 . 1 5837E-03 

-13 

-0 . 32167E-04 

-0 . 13305E-03 

-0.14097E-03 

-0 . 18022E-03 

-12 

-0 . 25633E-04 

-0 . 15408E-03 

-0 . 16223E-03 

-0 . 20350E-03 

-11 

-0 . 18248E-04 

-0 . 17815E-03 

-0 . 18615E-03 

-0 . 22797E-03 

-10 

-0 . 10275E-04 

-0 . 20568E-03 

-0.21296E-03 

-0 . 25325E-03 

-9 

-0 . 20334E-05 

-0 . 23715E-03 

-0 . 24292E-03 

-0 . 27886E-03 

-8 

0.61165E-05 

-0 . 27306E-03 

-0 . 27630E-03 

-0 . 30413E-03 

-7 

0 . 13792E-04 

-0 . 31395E-03 

-0 . 31330E-03 

-0 . 32817E-03 

-6 

0 . 20614E-04 

-0 . 36030E-03 

-0 . 354 10E-03 

-0 . 34982E-03 

-5 

0 . 26217E-04 

-0 . 41253E-03 

-0 . 39877E-03 

-0 . 36771E-03 

-4 

0 .30280E-04 

-0 . 47088E-03 

-0 . 447 27E-03 

-0 . 38027E-03 

-3 

0 . 32538E-04 

-0 . 53537E-03 

-0 . 49940E-03 

-0 . 38598E-03 

-2 

0 . 32802E-04 

-0 . 60568E-03 

-0 . 55480E-03 

-0 . 38359E-03 

-1 

0 . 30978E-04 

-0 . 68114E-03 

-0 . 61293E-03 

-0 . 37249E-03 

0 

0 . 27074E-04 

-0.76066E-03 

-0.67309E-03 

-0 . 35316E-03 

1 

0 . 21216E-04 

-0 . 8427 5E-03 

-0 . 73445E-03 

-0 . 32740E-03 

2 

0 . 1 3655E-04 

-0.92560E-03 

-0.79608E-03 

-0 . 29841E-03 

3 

0.47493E-05 

-0 . 10072E-02 

-0 . 85697E-03 

-0 . 27047E-03 

4 

-0 .50464E-05 

-0 . 10854E-02 

-0.91611E-03 

-0. 24813 E- 03 

5 

-0 . 15231E-04 

-0 . 11584E-02 

-0 . 97246E-03 

-0 . 23523E-03 

6 

-0 . 253 14E-04 

-0 . 12244E-02 

-0.10251E-02 

-0 . 23399E-03 

7 

-0 . 34866E-04 

-0 . 12822E-02 

-0 . 10731E-02 

-0 . 24455E-03 

8 

-0.43557E-04 

-0 . 13313E-02 

-0 . 11158E-02 

-0 . 26522E-03 

9 

-0 . 51165E-04 

-0 . 13713E-02 

-0 . 11528E-02 

-0 . 29317E-03 

10 

-0 . 57555E-04 

-0 . 14028E-02 

-0 . 1 1840E-02 

-0 . 32524E-03 

11 

-0 . 62632E-04 

-0 . 14264E-02 

-0 . 12096E-02 

-0 . 35869E-03 

12 

-0 . 66298E-04 

-0.14433E-02 

-0 . 12299E-02 

-0 . 39 1 53E-03 

13 

-0 .68437E-04 

-0 . 14548E-02 

-0.12458E-02 

-0 . 42255E-03 

14 

-0 .68900E-04 

-0.14620E-02 

-0 . 1 2580E-02 

-0 . 45121E-03 

15 

-0 . 67540E-04 

-0 . 14665E-02 

-0 . 12675E-02 

-0 . 47738E-03 


Only modest changes in the already small distortion values 
are evident. 
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3 . 2 Sensitivity to Span Load Shift 

Naturally, the flow over the test model will not be 
exactly the same as that computed and used in setting the flex 
walls. The effect of any differences can be estimated by 
considering possible alterations in the local distributions for 
a given fixed lift coefficient. In the basic calculations the 
span of the wing trailing vortex pair was set to correspond to 
an elliptic span load distribution. To determine sensitivity 
to span load alteration the vortex pair spacing, set by the 
variable, SV, in AFMODL was reduced from it/4 times the wing 
span to 2/3 the wing span; however, the wall contours were held 
at the values set for the basic calculation. The resulting 
flow distortion for comparison with that of the basic case is 
presented in the following table: 

Table 5 



NX= 40 NZ= 

20 MACH=0 .770 

A=0 .330 D=0 . 330 

M= 20 


PHIXM 

PHIZM0 

PHIZM1 

PHIZM2 

-10 

-0 . 19797E-03 

0 . 19535E-03 

0 . 19034E-03 

0 . 16217E-03 

-9 

-0 . 17327E-03 

0 . 22410E-03 

0 . 21795E-03 

0 . 18447E-03 

-8 

-0 . 14063E-03 

0 . 25440E-03 

0 . 24677E-03 

0.20811E-03 

-7 

-0 . 10121E-03 

0 . 28463E-03 

0.27513E-03 

0.23196E-03 

-6 

-0 . 57348E-04 

0 . 31170E-03 

0 . 30007E-03 

0 . 25387E-03 

-5 

-0 . 12184E-04 

0 . 33057E-03 

0 . 31702E-03 

0 . 27072E-03 

-4 

0 . 30909E-04 

0 . 33428E-03 

0 . 31992E-03 

0 . 27920E-03 

-3 

0 .68808E-04 

0 . 31449E-03 

0 . 30186E-03 

0.27712E-03 

-2 

0.98958E-04 

0 . 26286E-03 

0.25630E-03 

0 . 26470E-03 

-1 

0 . 11951E-03 

0 . 17313E-03 

0.17863E-03 

0 . 24434E-03 

0 

0 . 12944E-03 

0 . 43387E-04 

0 . 67660E-04 

0 . 21832E-03 

1 

0 . 12876E-03 

-0 .12227E-03 

-0.73454E-04 

0 . 18599E-03 

2 

0 . 11903E-03 

-0 . 31349E-03 

-0 . 23731E-03 

0 . 14193E-03 

3 

0 .10359E-03 

-0 . 51545E-03 

-0 . 4 1 344E-03 

0 . 75798E-04 

4 

0 . 87 1 18E-04 

-0.71241E-03 

-0 .59019E-03 

-0 . 20730E-04 

5 

0 .74140E-04 

-0 .89178E-03 

-0.75705E-03 

-0 . 14418E-03 

6 

0 . 67138E-04 

-0 . 10466E-02 

-0.90657E-03 

-0 . 28007E-03 

7 

0 .66036E-04 

-0 . 1 17 56E-02 

-0 . 10352E-02 

-0.41225E-03 

8 

0 . 69496E-04 

-0 . 12810E-02 

-0 . 11427E-02 

-0.53011E-03 

9 

0 .76572E-04 

-0 . 1 3665E-02 

-0.12312E-02 

-0 . 63020E-03 

10 

0 . 87040E-04 

-0 . 14364E-02 

-0 . 13038E-02 

-0 . 71393E-03 
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Comparison of the two relevant Tables 2 and 5, shows 
differences of low order and within acceptable limits. 


3 . 3 Sensitivity to Fore and Aft Load Shift 

To test sensitivity to fore and aft shifting of the lift, 
the sweptline vortex was moved from the wing quarter chord to 
the wing half chord. The wall flexure remained at the setting 
for the basic flow calculation. The flow distortion under 
these conditions are presented in Table 6 for comparison with 
the basic case in Table 2. 

Table 6 


NX* 40 NZ* 20 MACH*0 .770 A*0.330 D=0.330 M* 20 



PHIXM 

PHIZM0 

PHIZMl 

PHIZM2 

-10 

-0 . 34504E-03 

-0.89885E-04 

-0 . 97977E-04 

-0 . 13622E-03 

-9 

-0 .38631E-03 

-0 . 13871E-03 

-0 . 14840E-03 

-0 . 19393E-03 

-8 

-0 . 42445E-03 

-0 . 20122E-03 

-0 . 21285E-03 

-0 . 26589E-03 

-7 

-0 . 45633E-03 

-0 . 279 38E-03 

-0 . 29338E-03 

-0.35379E-03 

-6 

-0 . 47747E-03 

-0 . 37502E-03 

-0 . 39170E-03 

-0 . 45894E-03 

-5 

-0 .48176E-03 

'-0 . 48939E-03 

-0 . 50875E-03 

-0 . 58144E-03 

-4 

-0.46157E-03 

-0 . 62255E-03 

-0 . 64386E-03 

-0 . 71881E-03 

-3 

-0 . 409 10E-03 

-0.77266E-03 

-0 . 79392E-03 

-0.86419E-03 

-2 

-0 . 31875E-03 

-0 . 93543E-03 

-0.95285E-03 

-0 . 10051E-02 

-1 

-0 . 19024E-03 

-0 . 11040E-02 

-0.11116E-02 

-0.11242E-02 

0 

-0 .31089E-04 

-0 . 12692E-02 

-0 . 12594E-02 

-0 . 12025E-02 

1 

0 . 14291E-03 

-0 . 14210E-02 

-0 . 13856E-02 

-0.12251E-02 

2 

0 . 31042E-03 

-0 . 15496E-02 

-0 . 14823E-02 

-0 . 11892E-02 

3 

0 . 45004E-03 

-0 . 16477E-02 

-0 . 1 5459E-02 

-0 . 1 1097E-02 

4 

0 . 54637E-03 

-0 . 171 18E-02 

-0 . 15783E-02 

-0 . 10141E-02 

5 

0 . 59356E-03 

-0 . 17435E-02 

-0 . 15851E-02 

-0.92637E-03 

6 

0 . 59499E-03 

-0 . 17493E-02 

-0 . 15740E-02 

-0.85758E-03 

7 

0 .56030E-03 

-0 . 17380E-02 

-0 . 15530E-02 

-0 . 80778E-03 

8 

0 . 50211E-03 

-0 . 17182E-02 

-0 . 15283E-02 

-0 . 77220E-03 

9 

0 . 43324E-03 

-0 . 16964E-02 

-0 . 15047E-02 

-0.74657E-03 

10 

0 . 36436E-03 

-0 . 16769E-02 

-0 . 14846E-02 

-0 . 72886E-03 


Once again the variations are seen to be minimal 
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4.0 CONSIDERATION OF THREE DIMENSIONAL FLEXING OF THE UPPER 

AND LOWER WALLS 

The two dimensional flexing previously considered in this 
report appears to lack the ability to prevent some residual 
aerodynamic wash-in (induced twist of the wings). It seemed 
therefore that three dimensional flexing might improve the 
situation. One practical scheme considered envisioned a single 
row of jacks along the flex wall centerline with the lateral 
edge of the wall plate held fixed but not clamped. The jacks 
would conform the plate at each axial jack location to satisfy 
the Tracor blockage algorithm. To provide input to the flow 
distortion program PHIXZM the plate slopes in the axial 
direction were required at every panel center. For this a 
program called JACK_DISPL was developed using classical thin 
plate theory. A description of the method used and a Fortran 
listing are presented in Appendix D. 

Application of the analysis to the Langley 0.3-Meter TCT 
flexwall tunnel with the AEDC model at Cl * 0.5 indicated as 
expected the usual reduction in blockage and induced upwash; 
however, the spanwise variation of upwash was not improved at 
all even though it was basically small, e.g., only about 0.03 
degrees in a tunnel quarter width away from the center. The 
output data are presented below in Tables 7 and 8. 
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Table 7 - Rigid 

Walls 




MACH=0 CL-0 

.5 


IX 

PHIXM 

PHIZMO 

PHIZM1 

PHIZM2 

-10 

0 . 193E-02 

0 . 105E-02 

0 . 107E-02 

0 . 114E-02 

-8 

0 . 268E-02 

0 . 147E-02 

0 . 149E-02 

0 .158E-02 

-6 

0 . 344E-02 

0 . 195E-02 

0 . 198E-02 

0 . 211E-02 

-4 

0 . 414E-02 

0 . 248E-02 

0 . 253E-02 

0 . 270E-02 

-2 

0 . 473E-02 

0 . 305E-02 

0.311E-02 

0 . 335E-02 

0 

0 . 513E-02 

0 . 362E-02 

0.370E-02 

0 . 403E-02 

2 

0 . 531E-02 

0 . 417E-02 

0 . 427E-02 

0 . 469E-02 

4 

0 . 527E-02 

0 . 467E-02 

0 . 479E-02 

0.530E-02 

6 

0 . 504E-02 

0.511E-02 

0 . 524E-02 

0 . 582E-02 

8 

0 . 467E-02 

0 . 546E-02 

0 . 562E-02 

0 . 625E-02 

10 

0 . 424E-02 

0 . 574E-02 

0 . 591E-02 

0 . 659E-02 




Table 8 - Flexed 

1 Walls 




MACH=0 CL=0 . 

5 


IX 

PHIXM 

PHIZXO 

PHIZM1 

PHIZM2 

-10 

0.532E-04 

-0 . 738E-04 

-0.946E-04 

-0 . 172E-03 

-8 

-0 .137E-03 

-0.715E-04 

-0 . 104E-03 

-0.220E-03 

-6 

-0 . 274E-03 

0 . 287E-04 

-0 . 789E-04 

-0 . 251E-03 

-4 

-0 . 329E-03 

0.649E-04 

-0.988E-05 

-0 . 264E-03 

-2 

-0 . 285E-03 

0 . 204E-03 

0 . 978E-04 

-0 .264E-03 

0 

-0.133E-03 

0 .369E-03 

0 .225E-03 

-0 .261E-03 

2 

0 . 119E-03 

0 . 536E-03 

0 . 357E-03 

-0 . 255E-03 

4 

0 . 436E-03 

0 . 707E-03 

0 .495E-03 

-0 . 234E-03 

6 

0 . 752E-03 

0 . 892E-03 

0 . 649E-03 

-0 . 183E-03 

8 

0 . 991E-03 

0 . 108E-02 

0.815E-03 

-0 . 102E-03 

10 

0 . 110E-02 

0 . 126E-02 

0.971E-03 

-0 . 161E-04 

It 

is probable ! 

that the lack of 

displacement 

of the plate 


edges reduces the effectiveness of the 3-D flexure. 
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5.0 ESTIMATION OF FLOW DISTORTION FROM TEST DATA 


As mentioned in the introduction the concepts for using 
measured data on a suitable control surface during a test have 
been developed for some time. For application to hybrid 
conditions of partial adaptation the following analysis is 
presented to identify some functions used in the estimation of 
distortion that have the property of being calculable without 
iterative or inverting procedures that often introduce errors 
when applied to large matrices such as those needed for panel 
methods using large numbers of panels. 


Consider a solid wall tunnel or a ventilated tunnel where 
wall pressures and corresponding normal velocity can be measured 
or determined during a test. It is assumed that nothing is 
known of the model flow field. Disturbance velocities, axial 
and normal to a control surface coincident with original walls 
are denoted $ x and $ n respectively and the following cases 
are defined: 


o , o 
<p x and $ n 

$ x and 

o i 

6 $ x ci n c t 5 x 


Hn 


n 


are values at the wall for the wall-free 
condition (fully adapted), 
are the values measured during the test at 
the control surface (partially adapted), 
are the outer and inner change in $ x produced 

by the residual wall presence in its partially 
or non adapted state. 

is change in $ n from partial to fully adapted 
state . 

is defined positive inward. 


- 21 - 


Thus we may write 



(5) 


and 


**n m *n~*Z (6) 

Our goal is the computation of $ n or 6$ n for use in 
estimation of the residual wall interferences. The flow 
produced by the model in the region outside the control surface 
can also be considered to be produced by a distribution of 
normal velocity, $ n over the control surface and we may write 
the matrix formula: 

♦SCil s *°<j> G°(i,j) (7) 

where G° is a Green's source function for the outer flow and i 
and j are field point and source location indices. 

The change in the internal flow, 6$*, can be related to the 
normal velocity change, 6$ n as 

6**(i) = 6 * n (;j) Ql(i ' ^ • (8) 

where G* is a Green's source function. 

Introducing Equations (7) and (8) into (5) yields 
6 * n (j) G i (i,j) = -♦£ (j) G°(i,j) + ♦J(i) 


(9) 



- 22 - 


and using (6) to eliminate 6$ n we obtain 

♦°(j) - ~ ♦n (j) tG°(i,j) - G i (i,j)]- 1 (10) 

This wall-free value is not exactly correct because the 
measurements are made with the model in a partially adapted 
state. It represents the wall-free flow about a model tested in 
a slightly non-uniform stream. However, the residual errors for 
practical cases can be very small. 

The Green's function, G 1 , can be generated directly by 
means of the method of images; however, the function, G°, 
cannot. Fortunately as shown below the inverse function 
[G° - G*-] -1 can be generated directly. 

The wall presence in its partially adapted state can also 
be represented by a distribution of doublets or vortex elements 
over the tunnel wall control surface. This representation 
produces a flow change that is discontinuous in tx and 
continuous in $ n at the control surface and we may write: 

- r( j ) (11) 

where r is the local vortex intensity on the control surface. 
Also 

6 *n “ 5 6 ♦n (12 > 

Outside the control surface Equation (7) applies also to the 
increment, 6$ x and 

«*°(j) * «* n (i) G°(i,j) (13) 
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Substitution of Equations (8) and (13) in (11) leads to 

r(j) * 6$ n (i) tG i (i,j) - G°( i , j ) ] (14) 

and by matrix inversion 

6 * n (i) * r( j) [G i ( i , j ) - G°(i,j)r 1 f (15) 

but a relationship between r and 6$ n can be written as 

6* n (i) » - r ( j ) A (i,j) (16) 

where A is a function uniquely determined by the control surface 
geometry. It is indeed the normal velocity produced at a point 
on the wall, i, by a unit vortex element at j . It can be 
computed directly without images. 

Comparison of Equations (15) and (16) yields 

[G 1 - G°] “ 1 = -A (17) 

Finally from (10) and (17) 

♦„ - <♦" " ♦“ G 1 ) tA] (18) 

Equation (18) shows that the free air normal velocity 
distribution can be determined from the measurements and two 
directly calculable arrays, G 1 and A. No iteration nor matrix 
inversion is required, thus a high degree of accuracy can be 
assured when using large numbers of elements. It is fast and 
straightforward to use elements or panels with sizes as small as 
one or two inches in an eight foot square wind tunnel. 
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To obtain the flow distortion at the position of the model, 
i.e., pressure and flow angle and their gradients, we may use 
the following expressions: 

♦ x * U n " ♦n 1 ^ (19) 

♦ n = ( *n " ♦S ) l * m] (20) 

where G m and H m are Green's functions for wall sources. 

Both functions are developed by direct calculation using the 
method of images. 

The complete code for computation of $ n has been 
developed; it is called VEEXPHINO and a Fortran listing is 
presented in Appendix E. It requires the input functions 
measured wall normal velocity, MEASVN.DAT, and measured wall 
axial disturbance velocity, MEASVX.DAT. These functions must be 
developed to provide values at the centers of the panels chosen 
to correspond to the panel system chosen for VEEXPHINO. The 
output of VEEXPHINO called PHINWALL.DAT serves as input to 
PHIXZM for calculation of the residual flow distortion at the 
model location. If the system discussed is utilized in a 
practical case the pressures and normal velocities will probably 
be measured at a relatively small number of points. It is 
probably best to extrapolate the measured values to the more 
numerous chosen panel centers and then proceed rather than to 
limit the number of panels. The arrays G*- and A become very 
large if many panels are chosen, therefore, VEEXPHINO calculates 
smaller arrays from which G* and A are generated by using the 
geometrical similarities of the panel arrangements in a tunnel 
with a constant cross-sectional shape. 
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The first part of VEEXPHINO computes the first term in 
parenthesis in Equation (18). It uses the same technique of 
images used for PHIXZM to derive the constituent array GO used 
to generate the function GMi,j). The second part of the code 
develops the matrix array A(i,j) by use of the equations for 
horseshoe vortices given in Reference 1. Each panel is assumed 
to contain a centrally located bound vortex and two trailing 
vortices that extend downstream to infinity. The equations are 
shown and described in some detail in the listing. 

As mentioned previously, when the walls are not fully 
adapted the procedure described is only approximate in that the 
model produces a flow field influenced by the residual wall 
presence. However, the numerical experiments performed with 
theoretical models shows that partial adaptation in practical 
cases reduces the flow distortions to negligible values. 
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6.0 CONCLUSIONS 

The results of the computations indicate the power of 
flexing upper and lower walls of a rectangular wind tunnel in 
reducing the flow distortions at a model test location. Axial 
pressure gradients can be reduced to negligible values and 
upwash and upwash gradients can be sharply reduced. Once the 
walls have been set for a given lift coefficient and^ using a 
good calculated approximation to the model flow field , the 
present results show that the residual flow distortions are 
insensitive to variations in span loading and fore and aft 
loading. Thus once set the data obtained can be quite accurate 
even though the flow about the model is somewhat different from 
that computed. The calculated residual flow distortion at the 
model offers a good measure of the quality of the data and when 
not too large can be used as a basis for corrections, i.e., 
small angle of attack corrections for the fuselage, an 
indication of the induced aerodynamic twist of the wings and 
correction for the induced tail angle of attack. 
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FIGURE 3 - SOURCE IMAGE SYSTEM FOR M = 3 
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FIGURE 4 - COORDINATE AND PANEL NUMBERING SYSTEM 



FICURE 5 - AEDC MODEL AND DIMENSIONS 
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APPENDIX A - PHIXZM FORTRAN LISTING 


Primary Symbols 

A 

D 

GO 


HO, HI, H2 


MACH 

PHIXM 


PHIZMO 

PHIZMl 

PHIZM2 

VN 


Tunnel height 
Tunnel breadth 

Axial disturbance velocities at 

y + z = 0 produced by a Green’s source 

in the ring at IX = -NX 

Upwash velocities at z * 0, y * 0, 

y = A/8 and y = A/4 produced by the 

Green's source 

Mach number 

Tunnel centerline axial disturbance 
velocities produced by a distribution of 
normal velocities, VN, over all panels 
Centerline upwash due to the VN 
distribution 

Upwash at z * 0 and y • A/8 
Upwash at z * 0 and y * A/4 
Difference of calculated wall free 
velocity ratios and wall streamline 
slopes 



C-S o 


-35- 


C PROGRAM PHIX2M 

C THIS PROGRAM COMPUTES THE AXIAL VELOCITY INCREMENT 

C AND THE UPWASH IN THE MODEL VICINITY CAUSED BY HALLS 

PHIXM IS (+) FOR FLOWS DOWNSTREAM, AND PHIZM(UPWASH) 
IS (4) FOR INDUCED UPWARD FLOW AT THE MODEL, 

C A RECTANGULAR TUNNEL IS ASSUMED OF HEIGHT=A AND 

C WIDTHS 

C SELECT AN EVEN' ' NUMBER OF PANELS ON THE VERTICAL 

C WALL(NZ> .PANEL HEIGHT IS THUS A/NZ. PANEL LENGTH IS 

C SET AT A/NZ ALSO. TUNNEL LENGTH IS SET BY NX, THE 

C NUMBER OF PANELS FOREWARD AND AFT OF THE ORIGIN. 


J PANEL SIZE ON THE WALL IS A/NZ BY A/NZ, SQUARE! 

C ON THE FLOOR IT IS A/NZ BY D/NY 

C PRANDTL-GLAUERT COMPRESS I BILTY CORRECTIONS ARE USED 

C PROGRAM REQUIRES AN INPUT FUNCTION VN REPRESENTING 

C THE NORMAL VELOCITIES THAT MUST BE CANCELED AT THE 

HALLS FOR ALL THE PANEL CENTER POSITIONS. THE VN ARRAY 
C CAN BE CALCULATED FROM THE PROGRAM PHINO OR FROM DATA 

C COMPUTED FROM AN AIRPLANE CODE LIKE THE BOPPE CODE OR 

C OTHER SIMPLER REPRESENTATIONS. 


DIMENSION VN(80»-80580)»X(-80580)» Y(80)»Z(80) ,60(40,-80520), 
1 H0(40, -80520), Hl(40, -80520), H2(40, -80520), PT(-40540, 40), 

1 PJ(-40540, 40), Q<0!40, 40), 01(0540,40), 02(0540, 40), 

I R(0:40,40),R1(0540,40),R2(0:40,40),S1(0540,40), 

1 S2(0540,40),T1(0:40,40),T2(0540,40), PHIXM(-20520), 

1 PHIZM0(-20!20),PHIZM1 (-20520), PHIZM2(-20520) 

C NOTE! VALUES IN THE DIM STATEMENT ABOVE CORRESPOND TO 

C A SQUARE TUNNEL WITH NZ=20, SOME '20' VALUES ARE THE LATERAL 

C EXTENT OF IMAGES (M) TAKEN BEFORE SMEARING 

C SOME MUST CHANGE WITH THE VALUE SELECTED 

C FOR NZ. 


PARAMETER (PI=3. 14159) 

REAL MACH 
TYPE 10 

]:) FORMATUOX, 'ENTER NX NZ MACH A D AND M') 

READ(5,t) NX, NZ, MACH, A, D,M 
BETA=SQRT ( 1-MACHtMACH ) 


NZ2=NZ/2.0 

NY2-NINT(D<NZ2/A40,1 ) 

NY=NY2*2.00 

DYDZ S D/A$NZ/NY 


C’XViOli Ihw «1‘ .{ 7.' JlLi Ll.-' 

OF, POOR QUALITY 


C NOW SET UP THE COORDINATES OF THE PANELS. 


DO IX=-NX»NX 
X(IX)-IXtA/NZ/BETA 
END DO 
DO IZ=l»NY2 
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Y(IZ)=(IZ-0.5)tD/NY 

hid -M2 

END DO 

DO IZ=(NY2+1 ) » (NY2+HZ) 

Y(IZ)=D/2 

Z(IZ)=<0,5t<NZ+NY+l)-IZ)*A/NZ 
END DO 

DO I 1- ( NZ+NY2+1 ) * ( NZ+NY ) 
Y<IZ)=((NZtNY+0.5)-IZ)ID/NY 
Z< IZ<--A/2 
END DO 

C SET UP SOHE CONHON REPEATING TERNS TO SAVE TINE 


DO JZ=1»NZ2+NY2 
YD=Y < JZ) 

YDl=YD-A/8 

YB2=YP-A/4 

YPl=YD+A/8 

YP2-YDIA/4 

DO I=-N»H 

PT(I>JZ)"(2tAII+Z( JZ))t(2FA>I+Z(JZ)) 
PJU»JZ)=<2*A*I+A-Z(JZ))M2*AH+A-Z<JZ)) 

END 30 
DO L=0»N 
OL-DIL 
DF'DML+l) 

';(L..tZ)=<DL+YD)*(DL+YD) 

Q1 <L» JZ )=( DL+YD1 )*< DL+YD1 ) 

02 ( L > JZ ) s ( DL+ YD2) * ( DL+ YD2) 

R(L> JZ)= < DP- YD ) * ( DP-YD ) 
R1<L»JZ)=(DP-YP1)*(DP-YP1) 

R2 ( L r JZ ) = ( DP-YP2 ) < ( DP-YP2 ) 

SI (L» JZ)=< 0L+YP1 >> ( DL+YP1 ) 

S2 ( L > JZ ) = ( DL+YP2 ) t < DL+ YP2 ) 

TI(L) JZ) r (DP-YDl)*(DP-YDl) 
T2(L.JZ)=(DP-YD2)»<DP-YD2) 

END DO 
END DO 

C THE PRIMARY COMPUTATION BEGINS HERE 

DO KX=-NX»NZ2 

C THESE CALCULATIONS BELOU GIVE VALUES FRON THE FIRST 

C RING OF SOURCE PANELS AT IX=-NX 

XX=X(-NX)-X(KX) 

xs=xx*xx 

C FIRST UE CALCULATE THE SHEARED SOURCES 

EE=A/PI/D/NZ/NZ/BETA 

DH=(N+1)*D 

HN=<H+0.75)*2*A 

HN2=(N+0»25)*2tA 

GS=EE/BETA* < AT AN ( - XXtSQRT < XS+DNtDH4HNtHN ) / ( HMtBH ) )+ 
1 AT AN ( -XX4SQRT ( XS+ DM*DH+HH2*HM2 ) / ( HM2*DH ) ) ) 
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GS IS THE PHIX VALUES FROM THE SHEARED SOURCE 
IMAGES* 

G=<XS+HM2I*2> 

F =BH+SGRT ( DMM2+G ) 

G2 i (XSTHHtt2) 

F2=BMLSQRT<DH«2+G2> 

HS=-£E* ALOG ( FtSORT ( G2 ) / ( F2ISQRT ( S ) ) ) 

HS IS THE PHIZ VALUES FROM THE SHEAR. 


00 JZ= 1 >NY 2 
00 I=-HfH 
P 1 =PT<I,JZ)+XS 
P 2 =PJ(I»JZ)+XS 
HO L=0«H 
TT 1 =P 1 +Q(L» JZ) 

TT 1 =TT 1 *SQRTCTT 1 ) 

TT 2 -PltQl<L»JZ) 

TT 2 ^TT 2 *S 0 RT<TT 2 ) 

TTJif H 02 <L«JZ) 

TT 3 =TT 3 >SQRT(TT 3 ) 

TT4=P2+R(L»JZ) 

TI 4 =TT 4 *SQRT<TT 4 ) 

f! 5 -P 2 +Rl(L>JZ) 

TT 5 =TT 5 IS 0 RT(TT 5 ) 

TTc=P2+R2<L.JZ) 

TT 6 =rT 6 tSQRT<TT 6 ) 

IT 3 =PHR 1 (L»JZ) 

TT 6 =TI 8 ISQRT(TT 8 ) 

IT?=PHS 2 (L.JZ) 
rrR-rr?<soRT(TT?) 
i Ili=P 2 +Tl(LiJZ) 

TT 1 1 =TT 1 1 *SQRT ( TT 1 1 ) 
fT 12 =P 2 +T 2 (L» JZ) 

TT 12 s TT 12 tSQRT(TT 12 ) 

QQ=A»A/( 2 *PI*B£TA*NZINZ) 

TTT= < 2 /TT 1 + 2 /TT 4 ) #QQ 
RR=<l/TT 2 +l/TT 5 +i/TT 8 +l/TTll)IQ 0 
SS=( 1 /TT 3 + 1 /TT 6 + 1 /TT 9 + 1 /TT 12 )* 0 Q 

G 0 ( JZ * KX ) *-XX$TTT/BETA+GO< JZ * KX ) 
H 0 (JZrKX)=At( 2 II+O» 5 )tTTTtH 0 (JZ»KX) 
Hl(JZrKX)=A*( 2 II+ 0 » 5 )IRR+Hl(JZ»KX) 
H 2 (JZ,KX)=A*< 2 *I+ 0 , 5 )tSS+H 2 <JZ.KX) 

£ND DO 
END DO 

THE TERM DYDZ APPEARS BELOW TO ACCOUHT FOR THE 
DIFFERENCE IN PANEL WIDTH IF ANY OF ROOR AND WALL PANELS 


G0<JZ*KX)=<G0<JZ»KX)+6S)*DYDZ 
H0( JZ»KX)=(HO(JZ>KX)+HS)IDYDZ 
H1(JZ.KX)=(H1(JZ»KX)+HS)>DYDZ 
H2(JZ»KX)=<H2(JZ>KX)FHS)IDYDZ 
END DO 
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DO JZ=NY2+1»NY2+NZ2 
DO l=-M,« 

P1=PT(I»JZ)+XS 
P2=PJ(I>JZ)tXS 
DO L=OfH 
rTl-Pl+(ML»JZ) 

TTl=TTl*SGRT(TTl) 

TT2=P1+CI1<L>JZ) 

TT2=TT2*SQRT(TT2) 
f T3=F*l+Q2<L* JZ) 

T T3-TT3*SORT(TT3) 

IM=P2tR<l,JZ) 

IT 4=TTA*SQRT ( TT4) 

!T5=P2+R1(L>JZ) 

TT5=TT5*S0RT(TT5) 

TT6=P2+R2(L»JZ) 

TT6=TT6*SQRT(TT6) 

TTB=P1+S1(L» JZ) 

TT8=TT8*SQRT(TT8) 

TT9=P1+S2(L»JZ> 

TT9=TT9*SQRT<TT9) 

Tfll-P2m(LtJZ) 

TT11=TT1USQRT(TT11) 

ITl2=P2+T2(L» JZ) 

TT12=TT12*SQRT(TT12) 

T[T=(2/TTlf2/TT4)*00 
RRl=< l/TT2+l/TT8)tQQ 
PR2=( 1 / TT5+1/TT11 ) >QQ 
SS1=(1/TT3+1/TT9)*0Q 
SS2= < 1 . T T6+1 /TT12) *BQ 

GO ( JZ f KX ) =-XXtTTT/BETA+GO < JZ r KX ) 
H0(jZ.KX)=QQ*((2*A*I+Z(JZ))»2/TTl+(2tA»(I+.5)-Z(JZ))»2/TT4)+ 
1 HO( JZjKX) 

Hl(JZrKX)=<2*A«+Z(JZ))IRRH(2*A*(I+,5)-Z(JZ))*RR2+Hl(JZfKX) 

H2(JZ»KX)=(2*A*IH(JZ))*8Sl+(2*AI(I4*5)-Z(JZ))*SS2+H2(JZfKX) 

END DO 
END DO 

GO(JZ»KX)=GO(JZiKX)+GS 

HO<JZ»KX)=HO(JZ»KXHHS 

H1(JZ»KX)=H1<JZ*KX)+HS 

H2(JZ»KX)=H2(JZ>KX)+HS 

END DO 
END DO 

C THIS COMPLETES THE CALCULATION FOR SOURCES LOCATED 

C AT JX=-NX IN THE LONER HALF OF THE TUNNEL 

0PEN(UNIT=2*NANE =, PHINHALL.DAT'»STATUS='0LD') 

READ(2»*) ((UN(IfJ)» I=1»NY+NZ) * J=-NX»NX) 

CLOSE (2) 


jIwal pag? 

° P POOR Qu £ 


«r JS 

T 


DO IX=-NZ2»NZ2 
DO JX=-NX»IX 
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KX=-NX-< JX-IX) 

DO JZ=1*(NY2+NZ2) 

PHIXM(IX)=PHIXM(IX)-VN( JZ* JX)tGO(JZ*KX) 
F'HIZMO(IX)=PHIZMO(IX)-MN(JZf JX)tHO(JZfKX) 
PHIZM1(IX>-PHIZM1(IX>-VN(JZ*JX)IH1(JZ*KX) 
PHIZM2(IX)=PHIZM2UX)-VNUZ,JX)»H2(JZiKX> 
END DO 

DO JZ=(NY2+NZ2+1 ) * (NZ+NY) 

K2=MY+N2+1-J2 

HIXMm)=PHIXM(IX)-VN<JZ*JX>*GO<KZ*KX) 
r"H!?MO(IO=PHIZMOaXHW(JZfJX)*HO(KZ»KX) 
FHIZhiaX)=PHIZMl(IX)4VN(JZ»JX)*Hl(KZ*KX) 
F'HIZM2( IX >=PHIZH2< IX ) +VN( JZ» JX)*H2<KZ»KX) 
END DO 
END DO 

DO JX=(IX+1)»NX 
LX--HX+JX-IX 
DO JZ=1 » < HY2+NZ2) 

PHIXM' IX)=PHIXN< IXWNI JZ* JX)*GO< JZ»LX) 
PHI7M0<IX) S PHIZM0(IX)-VN( JZ*JX)>HO(JZ*LX) 
pHizpiax)=PHizmax)-VN(jZfjx)*Hi(jz»LX) 
PHIZ«2<:X)=PH:Z«2(IX)-W(JZ»JX)*H2(JZ»LX) 
END DO 

DO JZ=IN¥+NZ42)/2» (NZ+NY) 

KZ s NY+NZ+l-JZ 

PHIXN(IX)=PHlXn-;iX)tVH«JZ»JX}*GO(XZ»LX) 

PH I ZrtO < I X > =PM I r^OI X >♦ J2 . jx ) *HO( KZ »LX ) 

PH I ZN1 ( I X .» =PH ! Zh t f I X H VN ( jZ * JX ) *H1 < KZ • LX ) 
PHIZM2 X ) =PH! ZK < * y ; U’Ni D7. > JX ) *H2< KZ * LX ) 
END DO 
END DO 
END DO 


OPEN<UNIT=l <Z*.D hT -T;Tlj?='NEU') 
WRITE* 1 *47 J <CR» NX* « 2* a , D* 


47 

1 

FORMAT <5X* 'CR=' >F5*3* NX=' * E.i« «,>' -13* MAC*)*’ • 
'A*' *F5.3*'D=' *F5.3* 'M=' *1 
WRITE(1*50) 

•F5.3* 

50 


FORMAT ( lOX'PHIXM MOX * • PHIZMO •• OX* PH EZ;*1 • 1 C-< 

- FHIZM5M 



WRITE* 1*60 aX*PHIXMtIX**PHIZnOXlX)*P4K^Ox 

•*?m:M 2 «ix). 


1 

(X=-NZ/2*NZ/2> 


60 


FORMAT ( 7Xr I5» 4E15-5) 



CLOSE(UNIT=l) 


END 

POOR QUAXJCga 
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APPENDIX B - AFMODL FORTRAN LISTING 


Primary Symbols 
A 

AR & ART 
C & CT 

CIRC & CIRCT 

CL 

DX 

DZXL & DZXU 
MACH 
S, ST 
SV, SVT 
VN 

VNA 

VNS 

VNN 


WT , WTT 
XJ 


Tunnel height and breadth (assumed equal 
here) 

Wing and tail aspect ratios 
Wing and tail chords 

Circulation of wing and tail vortices 
Lift coefficient 

Extension of panels downstream of NX 
X-wise slope of flexwalls 
Mach number 
wing and tail spans 

Wing and tail vortex spans (spacing) 
Computed wall-free normal velocities at 
panel centers due to test model 
Anti symmetrical component of VN 
Symmetrical component of VN 
Difference between wall-free normal 
velocities and wall-slopes due to 
flexing or boundary layer growth 
Wing and tail air foil thicknesses 
Location of Langley 0.3-meter TCT 
flexwall jacks 
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PROGRAM AFMOBL *«■ 

COMPUTES THE MALL- FREE VELOCITIES(VN) NORMAL TO THE 
HALLS OF A SQUARE HIND TUNNEL OF AN AEDC MODEL 
TO BE USED IN A LANGLEY EXPERIMENTAL TEST. 

IT ALSO COMPUTES THE RESIDUAL NORMAL VELOCITIES AFTER 
APPLICATION OF THE TRACOR ALGORITHMFOR BLOCKAGE (VNN) 

IT PRINTS OUT A FILE CALLED VNZERO.BAT* THE HALL-FREE 
NORMAL VELOCITIES. FOR COMPARISON HITH THE VALUES FROM 
A NONLINEAR CODE SUCH AS 'TUNCOR'. 

THIS PROGRAM USES 2KNX+1 PLUS DX PANELS OVER THE LENGTH OF 
THE HIND TUNNEL AND NV+NZ PANELS AROUND THE HALF 
PERIMETER. TO INCREASE THE NUMBER OF PANELS 
MEH VALUES MUST BE PUT INTO THE DIMENSION STATEMENT 
AND THE LIMITS NX. DX AND NZ , 

THIS PROGRAM REPRESENTS THE MODEL BY THO SOURCES* BODY) 

VEE UNE-VORTICES<UINGiTAIL LIFT) AND +1- LINE SOURCES 
(HINGSTAIL THICKNESS). BODY LIFT BY LINE DOUBLETS 
IT ASSUMES THE CP IS LOCATED AT THk ZERO POSITION OF THE 
X COORDINATES 

THE HING IS ASSUMED TO CARRY 80Z OF THE TOTAL LIFT 

DIMENSION X(-80J120).XJ!-20IO).Y(80).Z<80).VNS<80.-80U20)» 
VNA<80>-80il20)»VN<80.-80!120)»DZXL(-B0*120)»DZXU(“80{120). 
Zl(-80tl2Q).ZU<-80:i20).VNN<80.-8Qll20) 

XJ(IX) ARE THE 1/3 METER TUNNEL JACK LOCATIONS 
X(IX).Y(IZ) AND Z(IZ> ARE THE CENTERS OF THE PANELS 

THE TRACOR ALGORITHM IS APPLIED IH THIS PROGRAM 
THE NUMBERS BELOH ARE FIXED BY THE DIMENSIONS OF THE 
AEDC MODEL 

PARAMETER <PI=3.H159.S=0.1O668* SV=. 0837863. AR*3. 5. 

D-0 .03556 . C= . 060198 » CT=0 . 03556 » W T s . 0072238 » WTT= . 004267 » 

ST* . 05588 > SVT* . 0X3888 > A= . 33 * ART=3.H29) 

REAL MACH 
TYPE 10 

FORMAT! 1 OX. 'ENTER CL MACH NX NZ') 

READ(5.t) CL. MACH. NX. NZ 
8ETA=SQRT ( 1-MACHIMACH) 

THE FOLLOWING RESULT FROM THE MODEL DIMENSIONS 

XRVC*-. 03759 ‘LOCATION OF VORTEX tC/L INTERSECTION(HfcTERS) 
XRV=XRVC/BETA 

XRVTC =0.0889 .'SAME FOR THE TAIL 
XRVI*XRVTC/BETA 

<SL=-. 10795 ! LOCATION OF FOREWARD SOURCE 

XST=0.1524 ! METERS. SANE FOR THE REAR SOURCE 

TG=. 57735/BETA ! TANGENT OF THE STRETCHED SWEEPANGLE 
REAL SWEEP IS 30 DEGREES! 


TS*1+TGITG 
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DX=0 

N72=NZ/2 

NY=NZ 

NY2=NY/2 


C NOU SET UP THE COORDINATES OF THE POINTS AT WHICH 

C VELOCITIES ARE TO BE CALCULATED. 


r 


DO IZ=1*NZ2 
Z(IZ)=A/2 

Y<IZ'-A/NZXIZ-0.5) 

END DO 

HO IZ=N22+1 » NY2+NZ 

mZ)=A/2 

Z < ! Z > = A/NZ4 < NZ+0 . 5- 1 Z ) 

END DO 

DO IZ=NZ+NY2+l,NZfNY 
Z t IZ >=-A/2 

,'UZ>=<NZ+NYF0.5-IZ)4A/NZ 
END DO 

00 IK--NX.NX+DX 

*UX)=A/NZ4IX ! FOR PANEL CENTERS 
END DO 

X J ! -20 ) =-27 . 75$ . 0254 

XJ<-19)=-224.0254 

<J(-18)=-174.0254 

XJ<-17)=-134.0254 

xj<-16)=-104,0254 

XJ(-15>=-84.0254 

XJ<-14)=-6. 54.0254 

XJ’! -13)=— 54.0254 

<J<-12)=-3.54t.0254 

XJ(-ll>=-24.0254 

<J{-10)=-. 54.0254 

XJ<-9)=1, 04.0254 

XJ(-8)=34.0254 

XJ(-7)=54.0254 

xj (-6 >=74,0254 

XJ(-5i=104.0254 

XJ:-4>=14*.0254 

XJ'-3>=194.0254 

XJ< -2) =244.0254 

XJ(-1)=294,0254 

xj(0)=34. 54.0254 ! ALL FOR JACK LOCATIONS 


PAG-v 


1 1 ? 


OF POOR QUALITY 


C!RC=0.84CL4S4S/(SV4AR> 

AL PHA=0 . 84CL/ ( 3 . 21 1 . 7554MACH446 ) !IN RADIANS 

THIS ASSUMES AN APPROXIMATE MACH DEPENCENCE FOR Cl/ALPHA 


Q=WT 

QT=WTT 

C I RCT= . 24CL4ST4ST/ ( SVT4ART ) ! TAIL CIRCULATION 

BBD=l+.54D4D4BETA4BETA/( (XST-XSL)4(XST-XSL) ) ! APPROX STREN6TH CORRECTOR 
DO IX=-NX*NX+DX 

STRETCHED DISTANCES FROM FIELD POINT TO INTERSECTION OF C/L AND? 

XV= < XRVC-X (IX)) /BETA ! WING VORTEX 

«L=<XRVC-C/4-X(IX))/BETA ! WING LEADING EDGE 

XT=( XRVCE34C/4-X ( IX ) ) /BETA ! WING TRAILING EDGE 
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XLS=IXSL-X(IX>)/BETA 
XTS=< XST-X ( IX ) ) /BETA 
XVT=(XRVTC-X(IX))/BETA 
XLT=XVT-CT/4/BETA 
X TT^XVT+ J*CT/4/BETA 


> BODY SOURCE 
1 BODY SINK (AFT) 

' TAIL VORTEX 
' TAIL LEADING EDGE 
• TAIL TRAILING EDGE 


ORIGINAL PAGE IS 
OF POOR QUALITY 


C RH IS THE ARH LENGTH FROM THE POINT OF ROTATION OF 

C THE BALANCE/STING SYSTEH TO THE CP ORIGIN OF COORDS. 

C THE MODEL IS DISPLACED UPWARDS BY A DISTANCE ALPHAtRM 

RN=O.5012o9 ! METERS 


DO IZ=1>NY2 
F=(XV+Y<IZ>*TG> 

F1=(XV-Y( 1Z)ITG) 

FT=< Xv‘T+y( IZiITG) 

FTT-«xVT-rdZ)*TG) 
n='M (Z>M2 

Z Z I Z ) F ALPHAtRM ) *12 

-4»'FtF+TStZZ) 
lK=4t(FltFlFTSIZZ) 
M=4«.‘,<L+YUZ)tTG)tt2+TS*ZZ) 

04*4* U XL-Y ( IZ ) tTG) W2FTSIZZ ) 

( XT +Y ( IZ)tTG)«2+TSIZZ) 
Q6MI((XT-Y(IZ)tTG)tt2FTStZZ) 
Q1T=4*<FT*FT+TS*ZZ) 

02T=4*< F 1 T*F 1 T+TS*ZZ ) 

03T=4*( < XLT+Y ( IZ )*T6)**2+TS*ZZ) 

Q4T a 4t< (XLT-Y ( IZ)*TG) W2+TStZZ) 
iT5TMt(<XTTFY(IZ)tTG)t*2+TS*ZZ) 

Q6T«4* ( ( XTT-Y< IZ)tTG) **2+TS*ZZ) 
Rl=S0RT(XVtt2FZZW> 

RESORT ( ( X VFTGtSV > H2FZZF < Y < IZ ) -SV ) **2 ) 

R21 -SORT ( ( XV+TGtSV) tt2+ZZ+( Y( IZ)+SV)tt2) 
R3»SQRT ( ( XLFTGtD/2 )**2+ZZ+< Y( IZ)-D/2)**2) 

R3 1 =SORT ( ( XL+ TGtD/2 ) t*2f ZZ+ ( Y ( IZ ) +D/2 ) 112 ) 
R4=SQRT ( ( XLFTGtS ) t*2FZZ+( Y ( IZ ) -S ) W2) 

R4 1 =SGRT ( < XLFTGtS) t*2FZZ+ < Y ( I Z ) f S ) tt2 ) 
P5=SQRT < < XT+TG*D/2 ) **2+ZZ+ < Y ( I Z ) -D/2 ) **2 ) 
R51*SQRT( ( XTFTGtD/2 ) W2FZZF ( Y< IZ ) FD/2 ) *12) 
R6 s SQRT((XTFTG*S)*t2FZZF<Y(IZ)-S)tt2) 

P6 1 =SQRT ( < XTF TG*S ) tt2FZZF ( Y ( IZ ) FS ) **2 ) 

R7»ZZFYYFXLStXLS 

R8=ZZFYYFXTStXTS 

R1T=SQRT ( XVTM2FZZFYY ) 

R2 T =SGRT ( < XVTF TGtSVT ) tt2FZZF ( Y ( IZ ) -SV I ) <12 ) 
R2 1 T *SQRT ( ( XVT FTGtSVT ) tt2FZZF ( Y ( IZ ) FSVT ) **2 ) 
WT»SOR T ( < XL TFTGtD/2 ) U2FZZF ( Y ( IZ )-D/2) *«2) 
R3 1 T =SQRT ( ( XLWGtB/2 ) t*2FZZF ( Y < IZ ) FD/2) «2 ) 
M T»Stt T < < XLT+TG*ST ) **2FZZF ( Y( IZ )-ST) **2) 

R41 T =SQRT ( ( XLTFTGtST ) t*2FZZF ( Y( IZ) FST ) 1*2 ) 
R5T=SQRT ( < XTT+TG*D/2 ) **2+ZZ+ ( Y ( I Z) -D/2 ) **2) 
R5 1 T=SQRT( ( XTTFTGtD/2 ) **2+ZZ+( Y ( IZ) FD/2)»*2) 
R6T=SQRT ( ( XTTFTGIST ) **2+ZZ+( Y ( IZ)-ST)**2) 
R61T=3URT( (XTTFTGtST)t*2FZZF<Y(IZ)FST)t*2) 


VNS(IZmX)=-D*D*DDD*(Z(IZ)+ALPHA*RH)/16/BETA*(1/(R7*S0RT(R7)) 

-0,61/<R8tSORT(R8>>> 
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THESE ABOVE ARE THE BODY SOURCE TERMS 

W 9 1 1 ; . IX >=-G*< 2(IZ) +ALPHA*RM) /PI/BETA* ( ( TS*S+XL*TG-Y< IZ ) ) / ( Q3*R4)- 
<fS*B/2fXL*TG-Y(IZ))/(Q3*R3>+<TS*S+XL*TG+Y<IZ))/(Q4*R41)- 
(TS*D/2+XL*TG+Y<IZ))/<G4*R31MTS*S+XT*TG-Y(IZ))/<Q5*R«>+ 
.TSID/2+XTITG-Y(IZ)>/'Q5tR:.)-aS*Sm*TGUaZ))/(Q6*R61)+ 

< rs*D/2f X T*rG+ Y < IZ > > / <U6*R5i » ) +VNS < IZ » I X > 

THESE WERE THE WING LINE SOURCE TERMS 

VNS ( IZ » I X ) =-OT* < Z < IZ > +ALPHA*Rrt > /PI /BETA* ( ( TS*ST+XLT*TG-Y ( IZ) ) / 
(Q3T4R4T i-( TS*D/2+XLT*TG-Y(IZ> )/(Q3T*R3T)+(TS*ST+XLT*TG+Y(IZ))/ 
(G4T*R41T)-(TS*D/2+XLT*TG+Y(IZ))/(Q4T*R31T)-(TS*ST*XTTirG- 
t'.I7)>/«35T$R6T»+<TS*0/2+XTT*TG-Y(IZ))/(Q5T*R5T)- 

< TS*3HXTT*TG+Y( IZ ) ) / ( Q4T *R61T ) + ( ISID/2+XTT *TG+Y ( IZ ) ) / 

< 06T*R51T > HVNSt IZfIX) 

THESE WERE THE TAIL LINE SOURCE TERMS 

VNA < I MX ) =D*D* ALPHA/8* < ( YY-ZZ )/ < YY+ZZ ) *121 ( 1 -XLS/SQRT (R7 ) ) 
tZZ*XL3/<YYRZZ)/(R7*S0RT(R7>)-0,61*UYY-ZZ)/(YY+ZZ)«2* 

< 1 -XTS,'SQRT(R8) ) f ZZ*XTb7( YYRZZ )/(R8*S0RT(R8) ) ) ) 

T HESE ARE THE DOUBLET ( BODY LIFT) TERMS 


T T =4 < Q 1 * '. ( TSISU+XU* TG- Y < I Z ) )/R2- ( X V*TG- Y 4 IZ) >/Rl) 
rO=4/H2*(<TS*SURXV*TG+Y(IZ) )/R21-(XV*TG+Y(IZ))/Rl) 

VNA( IMX)=VNA( IZ»IX)+CIRC/4/PI*<FtTT+Fl*TO+(Y(IZ)-SV)/ 
(ZZF(r<IZ)-SV)**2)*(t-(XV+SM*TG)/R2MY<IZ>*SV)/(ZZ+ 
<r<!Z>+SV)**2)*<l-(XVRSU*TG)/R21>> 

THE ABOVE ARE THE WING LINE VORTEX TERMS 

TTT-4/QiT*UTS*SVT+XVT*TG-Y(IZ))/R2T-(XVTtTG-Y(IZ))/RlT) 

T0T=4/Q2T*UTS*SMT+XVT*TG+Y(IZ))/R21T-<XVT*TG+Y(IZ))/R1T> 

VNA( IZ* IXi=VNA< IZ? IX)+CIRCT/4/PI*(FT*TTT+FlT*TOT4<Y(IZ)-SVT)/ 

1 ZZ+tY' IZ)-SVT)**2>*(1-(XVT+SVT*TG)/R2T)-(Y(IZ)+SVT)/ 
(ZZt«Y(IZ)+SVnt*2)*(l-(XVT+SVT*TG)/R21T)) 

THE ABOVE ARE THE TAIL VORTEX TERMS 

VNX.IZ? U)=VNA( IZ»IX)+VNS(IZ>IX) 

END DO 

DO IZ=NY2+NZ+1>NY4NZ 'UPPER WALL 
F=(XV+Y(IZ)*TG) 

F1=(XV-Y(IZ)*TG) 

FT=(XVT+Y< IZ)*TG) 

FTl=<xVT-Y(IZ)*TG) ORIGINAL PAGE IS 

nsT<IZ) ** 2 OF POOR QUAJUEQfl 

Z Z* i Z ( I Z ) FALPHAIRM ) *12 

Gl=4*<F*F+TS*ZZ) 

N2=4*’F1*F1+TS*ZZ) 

03=4«((XL+Y<IZ)*TG)**2+TS*ZZ) 

!Ct=4*U/.L-Y(IZ)*TG)**2+TS*ZZ) 

)j=4*( (XTFY(IZ)*TG)**2+TS*ZZ) 

< 6=4*’(XT-Y<IZ)*TG)*»2+TS*ZZ) 

01T=4*(FT*FT+TS*ZZ) 

IJ2T=4*(F1T*F1T+TS*ZZ) 

03T = 4* < ( XLT4Y ( IZ ) *TG ) **2+TS*ZZ ) 
i/4T-4*(ULT-Y(IZ)*TG)*<2+TS»ZZ) 

GZ ■ =4* ( ( XTHY ( I Z ) *TG ) **2+TS*ZZ ) 
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io T =4*iO(Tr-y< 1Z)ITG>**2+TS*ZZ) ORIGINAL’ PAGE IS 

R1-SGRT' <'J**2+ZZ+YYi OE POOR QUALITY 

'< 2 - S5R T < ( XV+TGISV ) * *2+ ZZ+ ( Y ( IZ ) -SV ) **2 ) 
R21=S0RT((XV+TG*SV)**2+ZZHY(IZ)+SV)M2) 

R3=$GR f « ( XL+TG*D/2 ) M2+ZZ+ ( Y < IZ ) -D/2 ) *12 ) 

R31-SQRT ( ( XUTGID/2 ) W2+ZZ+ ( Y ( IZ HD/2 ) **2 ) 

R4- SORT < ( XL+ TG*S ) **2tZZ+ ( Y( IZ ) -S ) **2) 

R4 1 -SORT < UL+TG*S > M2+ZZ+ < Y < IZ > +S > **2 ) 

R3*SIKT< <XWG*B/2)M2tZZ+< Y( IZ)-D/2)**2) 

R5 1 =S0RT< ( XTITGID/2 ) 1*2 IZZ+ < Y < IZ) ID/2 ) t<2 ) 

R p=SuR T < < XTTTGlS ) **2+ZZ+ ( Y < I Z ) -S ) 1*2 ) 

R6l=SGRT(<xT+TG*S)*t2+ZZHY(IZHS)**2) 

R/=::+YUXLS*XL5 

R8=ZZ+YY+XTS*XTS 

R1 I =S0RT ( X VTII2+ZZ+Y Y ) 

R2T=SQRT ( ( XVT+TGISVT ) W2+ZZ+ < Y (IZ ) -SVT ) **2 ) 

R21 r=SQR T ( ( X VT+TG*SVT > **2+ZZ+ < Y< IZ) +SVT ) **2 ) 

R>T--GRT< <XlT+TG*D/2)**2tZZ+(Y(IZ)-D/2)**2) 
r3lT--i)RT( (XLT+TG*D/2)**2+ZZ+(Y(IZ)+D/2)l*2) 

R4T =S0RT < < XL T FTG*ST ) II2+ZZT ( Y ( I Z ) -ST ) *12 ) 

R4 1 T -5QRT ( > XLT+TG*ST)**2+ZZ+(Y( IZ)+ST) **2 ) 

R3!=SGRT< ( XTT+TGtD/2 )**2+ZZ+ < Y( IZ)-D/2) **2) 

S5 1 T -SORT ( ( XTTITGID/2 ) * I2+ZZ+ < Y ( IZ ) +D/2 ) 1*2 ) 

R6T=SQRT < ( XTT+TG*ST ) II2IZZH Y f IZ ) -ST ) **2) 
R6lT=SGRT<<XTW6*ST)«2+ZZ+(Y<IZHST)**2) 


VNSUZ.iX)=+B*BtDDD*<Z(IZ)+ALPHA*RH)/16/BETA*U/<R7tSQRT(R7)) 

-0.6i/(R8*S0RT(R8))) 

THESE ABOVE ARE THE BODY SOURCE TERNS 

VNS( IZ*lX)=G*(Z(IZ)+ALPHAtRN)/PI/RETAt((TStSFXLtTG-Y<IZ))/(03IR4)“ 
i TS*D/2+XL*T6-Y < IZ ) ) / < Q3*R3 )+< TS*S+XL*TG+Y ( IZ ) ) / ( Q4IR41 >- 
( T SID/ 2TXLITGTY X IZ) )/< 04IR31 ) - ( TSIS+XTITG-Y ( IZ ) ) /(Q5IR6) + 
i TS*t</ 2+XTITG-Y < IZ ) ) / < Q5IR5 ) -( TS*S+XT*TG+Y< I Z ) ) /(Q6*R61)+ 

( TS*D/2+XT*TGtY( IZ ) ) / ( Q6IR51 ) ) +VNS< I Z » IX ) 

TfCSE HERE THE WING LINE SOURCE TERNS 

VNS i IZ > I X ) =0T* ( Z < IZ ) +ALPHA*RN ) /PI/BETA* ( ( TS*ST+XLT*TG-Y< IZ ) ) / 

< Q3TIR4T ) - ( TSID/2+XLTITG-Y ( IZ ) ) / ( Q3TIR3T ) + ( TS*STtXLT*TG+Y< IZ ) ) / 
ul4T*R41T)-<TStD/2+XLTITG+Y(IZ)>/<04T*R31T>-<TS*ST*XTT*T6- 
Y ( I Z ) ) / ( G5T IR6T ) +< TSID/2+XTTITG- Y ( I Z ) ) / ( Q5T IR5T ) - 
( TSIST WTITG+Y < IZ ) ) / ( 0&TIR61T ) + ( TSID/2 +XTT*TG+Y< IZ ) )/ 

( Q6T IR51T ) ) +VNS ( I Z r IX ) 

THESE WERE THE TAIL LINE SOURCE TERNS 

VNA< IZiIX)=-D*D*ALPHA/8*( (YY-ZZ>/<YY+ZZ)**2*U-XLS/SfiRT(R7) > 
+ZZ*XLS/(YY+ZZ)/(R7*SURT(R7))-0,61*((YY-ZZ)/(YY+ZZ)**2* 

( 1-XTS/SQRT (R8) ) +ZZ*XTS/( YY+ZZ)/(R8*SQRT<R8) ) ) ) 

THESE ARE THE DOUBLET 1 ! BODY LIFT) TERNS 


TT=4< 01 *( ( TS*SV+XV*TG~Y< I Z ) ) /R2- ( XV*TG-Y< IZ))/R1» 
T0=4/02*( < TS*SV+XV*TG+Y< IZ) > /R21- l XVITGFY < IZ ) ) /R1 ) 

VNA< IZ? IX)-VNA(IZ«IX)-CIRC/4/PI*(F*TTtFl*T0+(Y(IZ)-SM)/ 
-ZZHY(IZ)-SV)**2)*(1-(XVFSM*TG)/R2)-(Y(IZ)*SV)/<ZZ+ 
U(IZHSV)**2)*(1-<XV+SV*TG)/R21)) 

THE ABOVE ARE THE WING LINE VORTEX TERNS 
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TTT=4/fUT*<(T3*SVT+XVT*TG-Y<IZ))/R2T-<XVTtTG-YUZ))/RlT) 
T0T-4.02T*'. iTS*SVTfxVT*TG+Y<IZ))/R21T-<XVTtTG+Y(IZ))/RlT) 

VNA < 1 Z • I X * = VNA < I Z » I X ) -C I RC T/4/P I * ( FT*TTT+F 1T*T0T+<T(IZ) -SVT ) / 
■ZZ+tY(IZ)-$VT)**2>*U-<XVT+SVT*lG)/R2T)-<Y<IZHSVT)/ 

< ZZ * <■ Y ( I Z ) +SVT ) **2 > * ( 1 ■ - < XVT+SVWG ) /R2 1 T ) ) 

THE MOVE ARE THE TAIL VORTEX TERMS 


VN( IZ< IX»=TVNA< IZ'IXHVNSI 12 • IX > 


r.N[l DO 


DO :Z=Nr2+l.MY2+MZ 'COMPLETE SIDEWALL 

fY=Y t !)«2 

ZZ=<Z<IZ)+ALFHA*Rh*«2 

F-<XV+Y(IZUTG) 

F!=<XV-Y(IZ)*TG> 

FT=(XVT+YUZ»*TGi 

i-Tl*»XVT-Y«.lZ>tTG) 


t)l=4tlF*F+n*ZZ< 

V2=4*tFl*Ft+TS*ZZ) 

53=4* i < ti + r ( 1 z ) *TG) **2+TS*ZZ) 

P4= 4 * • ■ xl-y< IZUTG) **2+TS*ZZ ) 

55=4* < < XT+Y < IZ ) *TG ) **2ETS*ZZ ) 

D6=4*< (xt-y;iz)*tg)«2+ts*zz) 

rtiT=4*(FT*FT+TS*ZZ) 

iiZT-4*(FlT*FlT+TS*ZZ) 

G3T=4*'(XLT+Y(IZ)*TG)**2+TS*ZZ) 

<M T=4* ( ( XLT-r < IZ ) *TG )**2+TS*ZZ) 

Q5T=4*(( XTTIY < IZ) *TG) W2+TSIZZ ) 
t)6T-4*((XTT-Y(IZ)*TG)**2+TS*ZZ) 

R 1 =SOR T < X VW2+ZZ + Y < I Z ) «2 ) 

R2=S0RT(IXV+TStSV)«2+ZZ+(Y<IZ)-S»)«2) 

R2l=SQRTt(XV+TG*SV)**2+ZZW(IZ)+SV)**2) r/ , 

R3=SQRr((XL+rG*Ii/2.'**2+ZZE(Y(IZ)-D/2)W2) OF ^OCm '' ’ ' " T? 

R31=SGPT<(XL+TG*D/2)**2+ZZ+<Y(IZHD/2)«2) '* iUOi ^ 

R4=SGRT((XL+TG*S)**2+ZZHY(IZ)-S)**2) 

■?41=5QRT((XLTTG*S)**2tZZE(Y(IZ)+S)*l2) 

F 5=SORT < ( XT HGID/2 ) **2+ZZ+ < Y ( I Z ) • -D/2 ) **2 ) 

R51 =S0RT ( ! XT +TGID/2 ) **2+ZZ+ ( Y ( IZ ) +D/2) **2 ) 

Ro=SGF'T ( ( XT+TG*S ) W2+ZZ+ < Y( IZ)-S)**2) 

R61=SQRT< < XT+TG*S)**2+ZZ+< Y ( IZ )+S)**2) 

R7=ZZ+YY+XLS*XLS 

R8=ZZ+YY+XTS*XTS 

R1T=$QRTUVT**2+ZZ+YY> 

R2T =SGRT i. ( XVT+TGtSVT > W2+ZZ+ < Y ( I Z ) -SVT ) 1*2 ) 

R21 T-^prk ;xvr+TGtSVT)*l2+ZZHY(IZHSVT)f*2) 

R3T=S0RT ( < XLT+TGID/2 ) «2+ZZ+ ( Y( IZ ) -D/2) 1*2) 

R3;T=S')RThXLT+TG*D/ 2 ‘**2+ZZ+(Y(IZ)+D/2)t*2) 
R4T=SQRT((XLnTG*ST)W2EZZE<Y(IZ)-ST)«2) 

R4lT=S0Rri (XLT+TG*ST)**2+ZZ+(Y(IZ)+ST)*I2) 

R5T =SORT ( < XTT +TGID/2 ) **2+ZZt< Y ( IZ) -D/2 )*t2 ) 

F:5 1 T=SQR T < < XT T+ T6*D/2 ) **2+ZZ+ < Y ( IZ ) +D/2 )**2) 

R6 T=SQRT ( ( < TT+TG*ST ) H2EZZH Y < IZ ) -ST ) «2 ) 

RolT=SGRT( ( <TTTTG*ST)**2EZZ+(Y(IZ)+ST)lt2) 

VNS(IMX)-D*D*DDD/16/DETA*Y(IZ)*(1/(R7*SQRT(R7))- 
O.ol< TP8*S0RT(R8) > > 



THESE ARE THE BODY SOURCE TERMS OH THE WALL. 


0 £ 


'•NS< IZ« IX*=VN$( IZ»IXi-Q'PI/BETA<(((XL+Y(IZ)tTG)t(XL+SITG)+ZZ)/ 
Q3/R4- < < XL+Y ( IZ ) *TG ) t < XL+D/2*TG) +ZZ ) /Q3/R3- ( ( XL-Y ( IZ ) *TG) * 

< XL+S*TG>+ZZ >/Q4'R41+(<XL-Y< IZ) *TG) *( XL+D/2*TG)+ZZ)/ 

Q4/R31-( (XT+Y(IZ)*TG>l(XT+StTG)EZZ)/Q5/R6+( (XT+Y<IZ)tTG)l 
i < T+ 0/ 2*TG ) +ZZ ) / G5/R5+ ( (XT-Y(IZX16)I(XIESITG)+ZZ)/Q6/R61 
-<'XT-Y(IZ>*TG)*<XT+B/2ITG)+ZZ)/Q6/R5t) 

THE ABOVE ARE THE DING THICKNESS (LINE SOURCE) TERMS 

VNSi IZ. IX >=VNS( IZ. IX# -QT/PI.'BETAN ( <XLT+Y( IZ ) *TG ) * < XLT+ST*TG ) 
*U ) 03T /R4 T-: ( XLTW I Z ) *TG ) * < XLT+D/2*TG ) +ZZ ) /Q3T/R3T- 
<<XLT-Y(IZ)*TG»*(XLT+ST*TG)+ZZ)/Q4T/R41Tt((XLT-Y(IZ)*TG)* 
'XLT+D/2*TG>+ZZ)/Q4T/R31M(XTW(IZ)*TG)l<XTT+STtTG)+ZZ)/ 

Q3T /R6T *( ( XTT4-Y C I Z > TTG ) * < XTT+D/2*TG ) +ZZ > /Q5T/R5T+ ( ( XTT- 
i(IZ!*TG)»(XTT+ST*TG)+ZZ)/Q6T/R61T-((XTT-Y<IZ)ITG)*(Xm 
9/2 ft TQ » +ZZ » /Q6T/R51 T > 

THESE ARE THE TAIL THICKNESS TERMS 

W A ( 1 1 . 1 X i *- Dm ALPHA/8XY < IZ ) * < Z ( I Z ) +ALPHAIRM ) / < YY+ZZ ) *< < 1-XLS/ 
SrtRT ( R 7 > > *2/ < Y Y+ZZ ) - XLS/R7/SQRT <R7)-0.61*((I -XTS/SGRT (R8) > *2/ 
i YY+ZZ ' -XTS/R8/SQRT ( R8) ) ) 

THE ABOVE ARE THE DOUBLET TERMS (BODY LIFT) 

TTM./QIK <TS*SVEXV*TG-Y( IZ) )/R2-(XV*TG-Y< IZ) )/Rl) 
IOM/Q2*((TS*SV+XV*TGn<IZ))/R21-(XVtTG+Y<IZ)>/Rl) 

TTT=4/G1TK(TS*SVT+XVT*TG-Y(IZ))/R2T-(XVT*TG-Y(IZ))/R1T) 
fOT s 4/R2TI( ( TS*SVT+XVT*TG+Y(IZ))/R21T-<XVTITG+Y(IZ))/R1T) 

VNA(IZfIX)=VNA(lZfIX)+CIRC/4/PI*(Z(IZ)+ALPHA*RM)*(TG*(TT-TO) 

-l/(ZZ+(Y(IZ)-SV)tt2)*(l-(XV+SV*TG)/R2)+l/(ZZHY(IZ)+SM)ll2) 

t(l-(XV+SV*TG>/R21)> 

THESE ARE THE TERMS FROM THE SWEPT LINE VORTEX 

VNA ( I Z 1 1 X ) =+C I RCT/4/P 1 1 ( Z ( IZ ) EALPHAIRM ) ♦ ( TGI ( TTT-TOT ) - 1 / 

(ZZE< Y(IZ)-SVT)tt2)#(l-(XVT+SVTtTG)/R2T)+l/(ZZ+(Y(IZ)+SVT) 
»«)*(1-(XVT+SVT*TG)/R21T))+VNA(IZ»IX) 

THESE ARE FOR THE TAIL VORTEX SYSTEM 

VN(lZfIX)=VNA(IZ»IX)+VNS(IZ»IX) 

END DO 'COMPLETE TUNNEL IS NOW DONE 

END DO 

WE NOW APPLY THE TRACOR ALGORITHM AND COMPUTE THE 
SHAPE OF THE WALL TO ILLIMINATE BLOCKAGE. DZXL AND DZXU 
ARE THE SLOPES OF THE LOWER AND UPPER WALLS. 

DO IX=-NX.NX 
DO IZ s lfNY2+NZ2 

07XL(IX)=DZXL(IX)-VN(IZjIX)/NY2 
END DO 

THIS INTEGRATES THE NORNAL VELOCITIES OVER THE LOWER HALF 
OF THE TUNNEL AND DIVIDES BY THE LOWER WALL WIDTH 

DO IZ=NY2+NZ2+1»NY+NZ 
DZXUI IX)=DZXU( IX)+VN< IZ> IX)/NY2 


DO 12=1, NY2 

VNN<IZ.IX>=+VN(IZ,IX)+BZXL(IX) 
END DO 


DO !Z=NY2H.NY2+NZ 
VNN<!Z,IX)=FVN<IZ,IX; 

END DO 

DO t:=NY2TNZM.NY+NZ 
vMK< IZ » IX)-+VM< IZ » IX)-DZXU< IX) 

HD DO 
E-T DO 

this calculates the residual normal velocities on the tunnel 

4 ALL THAT MUST BE NEGATED BY THE GREEN'S SOURCES 


UE NON USE SIMPSONS RULE TO GET THE DISPLACEMENTS OF THE HALLS 


ZL<-NX)=0 

IL < -NX+ 1 >=< DZXL ( -NX ) FDZXL < -NX+1 ) ) IA/NZ/2 
ZU(-NX>=0 

ZU!-NX+1 ) =( DZXU< -NX ) FBZXU! -NX+1 ) 1IA/NZ/2 
DO IX=-NX+2»NX 

ZL<IX>=ZL<IX-2)+A/3/NZt(DZXL!IX-2H4tDZXL<IX-m 
DZXL(IX) ) 

ZU<IX>=ZU(IX-2)+A/3/NZ*<DZXU<IX-2H4*DZXU(IX-m 

DZXU(IX)) 

END DO 

OPEN<UNIT=?»NAME= ZDISPL .DAT' ,STATUS='NEH' ) 

WRITE!?, 200) 

FORMAT! 10X, '.33 METER TUNNEL AND AIRFORCE MODEL') 

WRITE<7,210) CL. MACH 

FORMAT <20X»'CL=' .F4.2.5X, 'MACH='F4»2) 

WftITE<7» 220) 

FORMAT! 10X, 'I '»8X»'ZL(I)'»8X,'ZU<I)') 

HRI*E<7,230) (I,ZL(I),ZU<I)»I=-NX,NX) 

F0RHAT!7X,I5»2E15,5) 

CLOSE (7) 

OPEN(UNIT=l»NANE='VNZERO.DAT'»STATUS s 'NEH' ) 

WRITE! 1,240) 

FORMAT! 10X, 'HALL FREE NORMAL VELOCITIES. LOWER HALF') 

URITE! 1.245) CL. MACH. NX. NZ 

FORMAT ( 30X» 'CL=',F4.2,5X»'HACH S ',F4.2»5X, 'NX='I4.5X» 'NZ='I4) 


DO IX=-NX.NX 

WRITE! 1.250) IX. ! VN< IZ. IX) . IZ=1 .NY2+NZ2) 

FORMAT < 2Xf I3» 2X r9E13» 4./* < 7X.9E13, 4) ) 

END DO 

WRITE! 1-255) 

FORMAT! lOX.'HALLFREE NORMAL VELOCITIES, UPPER HALF') 


DO IX=-NX,NX 

WRITE! 1.250) IX , ! VN( IZ, IX) , IZMNY2INZ2+1 ) , (NY4NZ) ) 



END DO 
CLOSE! D 


' f’EN ( UN I T r 2 > NAME = ' F'H I NWALL . DA T ' » S T A TUS- ' NEU ' ) 
URITE(2il) ( < VNN( IZ, IXb IZ=1 , (NY+NZ) ) »IX=-NX»NX> 
CLOSE! 2) 
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APPENDIX C - NONLVN FORTRAN LISTING 


Primary Symbols 


A 

AR & ART 
C & CT 

CIRC & CIRCT 

CL 

DX 

DZXL & DZXU 
MACH 
S, ST 
SV, SVT 
VN 

VNA 

VNS 

VNN 


WT , WTT 
XJ 

VNZERO 


Tunnel height and breadth (assumed equal 
here) 

Wing and tail aspect ratios 

Wing and tail chords 

Circulation o£ wing and tail vortices 
Lift coefficient 

Extension of panels downstream of NX 
X-wise slope of flex walls 
Mach number 
Wing and tail spans 

Wing and tail vortex spans (spacing) 

Computed wall-free normal velocities at 
panel centers due to test model 
Antisymmetrical component of VN 
Symmetrical component of VN 
Difference between wall-free normal 
velocities and wall-slopes due to 
flexing or boundary layer growth 
Wing and tail airfoil thicknesses 
Location of Langley 0.3-meter TCT 
flexwall jacks 

Wall-free normal velocities at panel 
centers from nonlinear flow computation 



PROGRAM NONLVN 


ORIGINAL PAGE 13 
QE POOR QUALITY 


THIS PROGRAM COMPUTES THE WALL SHAPES FOR THE 0.3 METER 
TUNNEL WITH THE USE OF A NONLINEAR CODE SUCH AS TUNCOR 
THAT PROVIDIES THE WALL -FREE NORMAL VELOCITIES IN A FILE 
CALLED VNZERONL.DAT. THE VALUES ARE COMPUTED FOR 
A MODEL SHAPE AND LIFT COEFFICIENT. COMPARISONS 
OF RESULTS WITH TEST DATA SHOULD BE MADE AT THE 
SAME LIFT COEFFICIENT! 

THIS PROGRAM USES 2INXF1 PLUS DX PANELS OVER THE LENGTH OF 
THE WIND TUNNEL AND NY+NZ PANELS AROUND THE HALF 
PERIMETER. TO INCREASE THE NUMBER OF PANELS 
NEW VALUES MAY BE NEEDED IN THE DIMENSION STATEMENT 
AND THE LIMITS NX, DX AND NZ . 

DIMENSION X(-80:i20),Y(80).Z(80),VNZER0(80, -80J80), 

DZXL ( -80 ! 1 20 ) , DZXU ( -60 : 1 20 ) , 
ZL(-80:i20),ZU(-8C:i20),VNN(SO,-80!120) 

PARAMETER' PI=3. 14159) 

REAL MACH 
TYPE 10 

FORMAT (1 OX, 'ENTER CL MACH NX NZ') 

READ<5,») CL, MACH, NX, NZ 
BETA=SGRT ( 1 -MACHIMACH ) 


DX=0 

NZ2=NZ/2 

NY=NZ 

NY2-NY/2 

NOW SE’ UP THE COORDINATES OF THE POINTS AT WHICH 
VELOCITIES ARE TO BE CALCULATED, 

DO IZ-1,NZ2 
Z(IZ)=A/2 

Y(IZ)=A/NZ*(IZ-0.5) 

END DO 

DO IZ=NZ2+1,NY21NZ 
Y ( IZ , =A/2 

Z(IZ>=A/NZ*(NZ+0.5-IZ) 

END DO 

DO IZ=NZ+NY2+1'NZ+NY 
Z(IZ)=-A/2 

Y( IZ > =( NZ+N Y+0 » 5- 1 Z ) *A/NZ 
END DO 

ItO IX=-NX,NX+DX 

X(IX)=A/NZ»IX ' ALL DIMENSIONS FOR PANEL CENTERS 
END DO 


0PEN(UNIT=1,NAME=' VNZERONL.DAT' > STATUS- OLD’ ' 
READ(l-l) (<VNZER0(IZ,IX),IZ=1.NY+NZ).IX=-NX,HX) 
CLOSE ( 1 ) 


WE NOW APPLY THE TRACOR ALGORITHM AND COMPUTE THE 



SHAPE OF THE HALL TO ILLIHINATE BLOCKAGE. DZXL AND DZXU 
ARE THE SLOPES OF THE LOWER AND UPPER WALLS. 


BO IX=-NX>NX 
BO 12=1 »NY2+NZ2 

DZXL (IX) =BZXL < IX )- VNZERO( IZ » I X ) /NY2 
ENB DO 

THIS INTEGRATES THE NORMAL VELOCITIES OVER THE LOWER HALF 
OF THE TUNNEL AND DIVIDES BY THE LOWER WALL WIDTH 

DO IZ=NY2+NZ2+1>NY+NZ ' UPPER HALF 

DZXU ( IX ) =BZXU ( I X ) T VNZERO (IZ<IXUNY2 
END DO 

DO IZ=1»NY2 

VNN(IZ>IX)*+VNZERO(IZ>IX)+DZXL<IX) 

END DO 

BO IZ=NY2F1jNY2+NZ 
VNN ( IZ > I X ) =+VNZERO< IZ» IX ) 

END BO 

BO IZ=NY2+NZ+1 »NY+NZ 
VNN(IZ»IX)=rVNZERO(IZ»IX)-DZXU(lX) 

END DO 
END DO 

THIS CALCULATES THE RESIDUAL NORMAL VELOCITIES ON THE TUNNEL 
WALL THAT MUST BE NEGATED BY THE GREEN'S SOURCES 


WE NOW USE SIMPSONS RULE TO GET THE DISPLACEMENTS OF THE WALLS 
ZL(-NX)=0 

ZL ( -NX+ 1 ) = ( DZXL ( -NX ) +DZXL ( -NX+ 1 ) )*A/NZ2 
ZU(-NX)*0 

ZU < -NX+ 1 ) = < DZXU ( -NX ) +BZXU( -NX+ 1 ) ) * A/ NZ2 
DO IX=-NX+2»NX 

ZL ( I X ) =ZL ( IX-2 ) +A/3/NZ»< DZXL ( IX-2 ) +HIDZXL (I X- 1 >+ 

DZXL (IX)) 

ZU ( IX ) »ZU ( IX-2 ) +A/3/NZt ( DZXU (IX-2 )+4>DZXU (IX-l)t 
DZXU(IX)) 

END DO 

OPEN (UNIT*?* NAME* ' ZDI SPL . DAT ' > ST ATUS s ' NEW ' > 

WRITE(7»200) 

F0RMAT(10X» '.33 METER TUNNEL AND AIRFORCE MODEL') 

WRITE(7i210) CL» HACH 

FORMAT < 20X» 'CL=' »F4.2r5X» 'MACH='F4.2) 

WRITE(7»220) 

F0RMAT(10X»'I'»8Xf'ZL(I)'f8X»'ZU(I)') 

WRITE(7.230) ( I » ZL< I ) » ZU( I ) » I=-MX»NX ) 

FORMAT ( 7X» I5» 2E15*S) 

CLOSE (7) 


OPEN( UNI T=2> NAME= ' PH I NWALL . DAT STATUS* 'NEW) 
WRITE(2,») ( (VNN( IZ*IX)»IZ=1» • NYTNZ) ) , IX=-NX,NX) 
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CLOSE* 2) 
END 
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APPENDIX D - JACKDISPL ANALYSIS AND FORTRAN LISTING 
Determination of Wind Tunnel Wall Displacements 

In this Appendix, we outline the analysis upon which the 
wall displacement program JACK_DISPL is based. The two primary 
outputs of this program are the displacements of both the floor 
and ceiling jacks needed to relieve the blockage and an 
estimated residual interference velocity normal to the floor and 
ceiling. This is then used as an input to the routine PHIXZM to 
estimate residual flow distortion at the model. 

The starting point is the basic relation adopted for 
relating the streamwise slope of the wind tunnel flexible wall 
(here floor and ceiling) to the wall-free normal velocities 
induced by the model. Thus for unit free stream velocity 

S/2 V n dA * J-b/2 tx dy (D,1) 

where V n = Inflow velocity due to model normalized on free 
stream velocity. 

w = Wall displacement from flat initial position (taken 
to be positive outward); w is a function of x and y. 
d£ = Differential length along the perimeter S of the 

tunnel cross section at each streamwise station x. 
x = Streamwise position along wall; x * 0 is taken here 
as being located at the beginning of the flexible 
wall. 

y « Spanwise coordinate; y * 0 is located at center of 
flexible wall of total span b. 

The integration on the left hand side of (D.l) is performed over 
the half perimeter symbolically denoted S/2. When considering 
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deflection of the floor, the S/2 is taken as the lower half of 
the tunnel; S/2 is taken over the upper half of the tunnel to 
determine flexure of the ceiling. This is the application of 
the Tracor Blockage Algorithm to the problem. 

We define a spanwise integrated displacement as 

w*(x) = dy * 2 f^ 2 w(x,y) dy (D.2) 

and the integrated normal velocity as 

V*(x) = / s/2 v n d£ * (D.3) 


combining 


(D.l 


- D.3) and rearranging, 

V* = - - fk/2 w dv = - 

v 3x ' -b/2 w ay 

w*(x) * -/* V*(x' ) dx' 

J O 


3w* 

3x 


( D.4 ) 
(D.5) 


Thus the spanwise integrated wall displacements have been 
expressed in terms of the input wall-free velocities normal to 
the tunnel walls. 

We now relate these values of w* to adjustments that are 
obtainable with a given single streamwise series of jacks that 
are located on the floor and ceiling at y » 0. To do this, we 
will model each flexible wall as a simply supported rectangular 
plate subjected to a concentrated ("point") load at each jack 
location. Inspection of (D.5) shows that it is not the load, 
but the displacements that are important. Thus we seek the 
displacement at each jack location, w(y = 0), that enables the 
resulting plate shape to satisfy (D.5). Loads will be employed 
only as intermediate variables used to obtain displacements. 

The actual loads needed for a given displacement depend strongly 
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on the details of the plate construction - stiffness, thickness, 
ribbing, etc. No attempt is made to calculate actual loads as 
they are not needed. 

The solution for the displacement of a simply supported 
rectangular plate of length a and width b subjected to a point 
load P at (x,y) * (c, 0 ) is given in Reference 14 for y > 0 as 


w(x,y;c) * 


2 w D m*l 


l { [ ( l+o_ tanh a m ) sinh ( ^ (b-2y)) 


- rp ( b-2y ) cosh ( ^ (b-2y)) ] 


sin ( ) sin ( ) 

fl q 

m 3 cosh a_ 


(D.6) 


Here, D is the plate flexural rigidity. We can define a 
dimensionless displacement due to a load of unit force as 


w (x,y;c) 


w(x,y;c) » — - l {} ( D. 7 ) 

2 ir 3 m«l 


where the summation expression is the same as in (D.6). 


Since the differential equation whose solution is given by 
(D.6) is linear, we can superpose solutions. Thus the displace- 
ment due to a series of loads * P(x * 5^) can be 
expressed 

_ "* A 

w(x,y) * J P i • w ( x ,y ; 5 ^ ) (D.8) 


where 
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P i 


a 


2 


D 


(D.9 ) 


a normalized load that has units of length. 

If w is evaluated at a finite number of discrete points (xj, 
y^) , relation (D.8) represents a set of linear algebraic 
equations , 

A 

w(Xj,y k ) - l w r (D.10) 


or in matrix notation, 



P 


(D.ll) 


w is a column vector whose, elements are arranged in the order 


(*l»yi»)» (xj ,y 2 ) , . . . ( x ! ,y N ),(x 2 ,yi) 

k 


( V y k’ 




Its length is Nj • where Nj and N^ are the total 
number of discrete x and y locations, respectively, of interest. 
P is a column vector of length Ni, the total number of 
concentrated loads, w is a rectangular matrix of length Nj # 

N)< and width Ni. It should be noted that in general, the 
locations {xj} and { c i } need not be the same. 

If the integration of (D.2) is applied to (D.6), another 
matrix equation is similarly found 


w* * w* • P 


(D.12) 
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where P is as previously defined, w* is of length Nj , and w* 
is of size Nj by N^. The denotes spanwise integrated 
quantities analogous to those of (D.ll). 

The slope of the flexible walls, 3w/3x (x,y), is found by 
differentiating (D.10) leading to an analog of (D.ll): 

w x = w x P ( D. 13 ) 


where the elements of w x are dw/dx. The expressions for the 
elements of the matrices w* and w x have been obtained by 
analytically integrating or differentiating the elements of w 
given by (D.6,D.7). 

The algorithm implemented in JACK_DISPL and its subroutines 
is now outlined. Calculations are repeated for the floor and 
ceiling . 


1. Using as input the normal velocities in free air at the 
wall locations, epxression (D.3) is calculated at each jack 
location by a Simpson's integration. 

2. w* is found at each jack location using a trapezoidal 
integration of (D.5). 

3. The elements of the various matrices in (D.ll - D.13) 
are calculated using expressions (D.6,D.7) and the appropriate 
integrated and differentiated forms of (D.6,D.7) defined above. 
The infinite series is truncated at, typically, 150 terms. 
Numerical experimentation has shown this to be well converged. 
Simplified expressions for the hyperbolic functions are used 
where appropriate for large arguments. 



4. P is obtained from (D.12) using w* from (D.5) and IMSL 
Subroutine "LEQT2F". 

5. The needed displacements at the jack locations are 
determined from (D.ll) with (x^y^) * ( x jack'°)* 

Displacements at other locations can be similarly evaluated. 


6. The "residual" or adjusted normal velocities are 
calculated using (D.13), and the expression (based on the 
original normal velocities and the effect of the sloping walls) 


n , residual 


(x,y) = V 


n ,or iginal 


(x,y) + 


3 w 
3 x 


(x,y) . 


( D. 14 ) 


In general, the signs of the two right hand terms will be 
opposite such that the new "residual" V n will be much less 
than the original. 


Fortran Listings 

Listings for the programs JACK_DISPL and AFMODLJ follow. 
AFMODLJ is a modified version of AFMODL that provides the input 
wall-free normal velocities at the panel and jack locations of 
JACK DISPL. 



-60- 


PROGRAM JACOISPL 

C THIS PROGRAM READS THE DATA FILE CVNJ.BAT] FOR THE LOCATION OF THE JACKS 

C IN THE HIND TUNNEL AND THE MAGNITUDES OF THE NORMAL VELOCITIES TO BE USED 

C FOR CALCULATING THE JACK DISPLACEMENT NEEDED TO RELIEVE THE EXCESS OUTHARD 
C FLOH IN THE TUNNEL. 

C 

C THE DATA FILE [UNO. DAT] IS THEN READ FOR THE DIMENSIONS OF THE UlNO 
C TUNNEL » THE NUMBER OF PANELS IN THE X. Y. l Z DIRECTIONS AND THE 
C LOCATIONS OF THE CENTER OF THE PANELS IN THE X J ' DIRECTION* THE WIDTHS 

C OF THE PANELS IN THE Y W DIRECTIONS. AND THE MAGNITUDE OF THE NORMAL 

C VELOCITIES AT THE PANEL CENTERS ON THE FLOOR AND CEILING. 

C 

C THE PROGRAM THEN CALLS THE SUBROUTINE 'SIMPSON* TO INTEGRATE THE NORMAL 
C VELOCITIES FOR THE LOWER HALF AND UPPER HALF OF THE TUNNEL AT EACH JACK 
C LOCATION USING SIMPSON'S INTEGRATION METHOD. THE SUBROUTINE 'MATRIX' IS 

C THEN CALLED TO CALCULATE THE MATRICES FOR THE JNIT POINT LOAD DISPLACEMENT * 

C THE UNIT POINT LOAD DISPLACEMENT DIFFERENTIATED WITH RESPECT TO X. AND THE 
C UNIT POINT LOAD DISPLACEMENT DIFFERENTIATED WITH RESPECT TO X AND INTEGRATED 
C WITH RESPECT TO Y USING A SUM OVER AN INFINITE SFRlES 
C TRUNCATING THE CALCULATION AFTER THE FUNCTIONS ARE DETERMINED 10 HAVE 

C CONVERGED. THE RESULTS OF THE SIMPSON INTEGRATION AND THE INFINITE SERIES 

C INTEGRATION ARE USED TOGETHER TO CALCULATE THE NORMALIZED LOAD AT EACH JACK 
C LOCATION USING THE IMSL SUBROUTINE 'LEQT2F*. USING THE NORMALIZED LOADS. 

C THE ACTUAL DISPLACEMENT OF THE JACKS IN THE PL DOR AND CEILING OF THE iijMit 

C TUNNEL ARE CALCULATED USING MATRIX MULTIPLICATION. THE PROGRAM THEN 
C CALCULATES THE RESIDUAL FLOW BASED ON THE DIFFERENCE BETWEEN ACTUAL 
C MEASUREMENTS OR THEORETICAL CALCULATIONS AND THE CALCULATED VALUES FROM THIS 
C PROGRAM. 

C 

C 

C THE ASSUMPTIONS OF THIS PR06RAN ARE AS FOLLOWS: 

C 1) THE TUNNEL IS SQUARE AND THE NUMBER OF PANELS IN I HE 

C Y-DIRECTON IS THE SAME AS IN THE Z-BIRECTION. AND THAT 

C DELTA. Y = DELTA.Z 

C 

C 2) THE NORMAL VELOCITIES ARE POSITIVE WHEN THEY FLOW INTO THE 

C TUNNEL. AND THEY ARE LOCATED AT THE CENTER OF THE PANEL. 

C 

C 3) THE DEFORMATIONS OF THE FLOOR AND CEILING OF THE TUNNEL ARE 

C SMALL AND CAN BE TREATED AS LINEAR. 

C 

C 4) THE JACKS LIE ALONG THE CENTERLINE OF THE TUNNEL «Y = 0) 

C AND THE X POSITIONS OF THE CEILING JACKS ARE THE SAME AS THE 

C FLOOR JACKS. 

C 

C 

C THE FOLLOWING VARIABLES ARE USED IN THE PROGRAM: 

C 

C DELTA. Y = REAL. WIDTH OF PANEL IN Y DIRECTION 

C DELTA.Z = REAL. WIDTH OF PANEL IN Z DIRECTION 

C BISPL-LO = REAL. ACTUAL DISPLACEMENT OF THE LOWER JACKS AS CALCULATED 
C BY MATRIX MULTIPLICATION 

C BISPLJJO = REAL. ACTUAL DISPLACEMENT OF THE UPPER JACKS AS CALCULATED 
C BY MATRIX MULTIPLICATION 

C DUMMY = CHARACTER. USED TO READ COMMENT LINES IN THE INPUT C ILE 
C IA = INTEGER. INITIAL DIMENSION SIZE OF THE INTEGRATION MATRIX U.INTEGO 

C ID = INTEGER. SPECIFIES THE ACCURACY OF THE ELEMENTS IN THE MATRICES 

C SENT TO THE IMSL SUBROUTINE FOR AN ACCURACY CHECK 

C IB = 0 INDICATES THAT AN ACCURACY CHECK IS NOT WANTED 
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C 

C 

c 

c 

c 

c 

c 

c 

r 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c. 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

n 

L 

c 

c 

c 

c 

c 

c 

c 

r 

c 


IER = INTEGER. OUTPUT FROM IMSL SUBROUTINE WHICH HAT INDICATE AN ERROR 
IER = 129 INDICATES AN ALGOR I THU I CALL f SINGULAR MATRIX 
IER = 34 INDICATES THE RESULTS FAILED AN ACCURACY TEST 


IZ = INTEGER* COMBINED NUMBER OF PANELS IN THE LOWER OR UPPER HALF OF 

THE TUNNEL IN THE COMBINED T i / DIRECTIONS. USED AS AN 

INDEX NUMBER 

N = INTEGER. NUMBER OF COLUMNS IN THE MATRIX ESTABLISHED IN THE SIMPSON 
SUBROUTINE 

NX = INTEGER. NUMBER OF PANELS IN THE X-DIRECTION 
NXSI = INTEGER. NUMBER OF JACKS IN THE TUNNEL 
NT = INTEGER. NUMBER OF PANELS IN THE Y-DIRECTION 

NZ = INTEGER. NUMBER OF PANELS IN THE Z-OIRECTION 

RESIB.Li.) * REAL .2-0. RESIDUAL FLOW AFTER SHAPING FLOOR 
RESID.UI.) * REAL.2-D. RESIDUAL FLOW AFTER SHAPING CEILING 
RLENGTH = REAL. LENGTH OF WIND TUNNEL 

SIMS.LO = REAL. INITIALLY THE INTEGRATED DISPLACEMENT (1ST INTEGRATED 
ONLY SPANUISE. THEN ALSO STTREAMWISE) AT EACH JACK FROM 
THE SIMPSON SUBROUTINE FOR THE FLOOR. AFTER IT RETURNS 
FROM THE IMSL SUBROUTINE IT IS THE NORMALIZED LOAD AT 
EACH JACK 

SIMS.UO = REAL. INITIALLY THE INTEGRATED DISPLACEMENT • 1ST INTEGRATED 
ONLY SPANWISE. THEN ALSO STTREAMWISE > AT EACH JACK FROM 
THE SIMPSON SUBROUTINE FOR THE CEILING. AFTER RETURNING 
FROM THE IMSL SUBROUTINE IT IS THE NORMALIZED LOAD AT 
EACH JACK 

U_INF=REAL» THE FREESTREAM VELOCITY IN THE WIND TUNNEL 
VNJ( . ) = REAL. 2-D. THE NORMAL VELOCITY AT THE PANEL POSITIONS A) THE 
JACKS 

VNO(.) = REAL. 2-D. THE NORMAL VELOCITY AT THE PANEL POSITIONS ON FLOOR 
AND CEILING 

W. DIFF ( » » ) = REAL. 3-D. POINT LOAD DISPLACEMENT DIFFERENTIATED WITH 

RESPECT TO X, CALCULATED BY 'INFINITE SERIES' 
W„BISP-L<») = REAL. 2-D. POINT LOAD DISPLACEMENT CALCULATED BY 'INFINITE 

SERIES' 

W_DISF'.U(>) = REAL. 2-D. POINT LOAD DISPLACEMENT CALCULATED BY 'INFINITE 

SERIES' 

WSTARLM = REAL. 2-D. SPAN INTEGRATED MATRIX. LOWER 

WSTARUX . ) = REAL. 2-D. SPAN INTEGRATED MATRIX. UPPER 

WORKSPACE! ) = REAL. DIMENSIONED WORK SPACE FOR THE IMSL SUBROUTINE 

WX_L(f) = REAL. 2-D. SLOPE OF FLOOR OF WIND TUNNEL 

WXJJ(.) = REAL. 2-D. SLOPE OF CEILING OF WIND TUNNEL 

X. LOCO = REAL. LOCATION OF CENTER OF PANEL FROM THE 1/2 INCH POINT 
XS-LOCO = REAL. X-LOCATION OF NORMAL VELOCITIES USED TO CALCULATE USTAR 
XSI.LOCO = REAL. LOCATION OF JACK IN X-OIPECTIGN RELATIVE TO THE 

1/2 INCH POINT 

Y. LOCO = REAL. LOCATION OF CENTER OF PANEL IN Y-DIRECTION RELATIVE TO 

CENTER OF TUNNEL 


IMPLICIT REAL (A-H.O-Z) 

IMPLICIT INTEGER (I-N) 

CHARACTER DUMMYIl 

COMMON/ S I MP_ INT/UN J ( 40 . 60 ) . SI MS_U ( 40 ) » SI MS_L (40) 

COMMON/LOCAT /X.LOC ( 100 ) • Y.LOC (30) » XSI _LOC ( 40 ) . XS.LOC ( 40 ) 
COMMON/W.HATS/W.D I SPL ( 40 . 40 ) • W..D l FF < 1 00 . 30 . 40 > . USTARL (40.40). 
1 WSTARU( 40.401 

DIMENSION WORKSP ACE ( 1 360 ; * RE S I B J. < 1 00 * 30 > » RE S I D.U ( 1 00 > 30 ) » 
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1 DI SPL-L ( 40 > » DI 3PL_U ' 40 ) . UX_L <1 00 » 30 > .. WX_U < IOC . 30) .. VMO f 1 00# 40) j- 

2 PHINDAT<4(M00) 

CALL GETCPU< NTIrtEl > 


C Reading info nation about jack iootiens [7NJ.0AT1 
0PEN< 2 i NAHE= ' VN J G STATUS* ' OLD ' > 

R£AD(2j '(T 18iI3)')NXSI 
READ ( 2 » ' ( A 1 ) ' ) DUMMY 
READ(2»*) ( XSI_L0C( I ) * 1=1 *NXSI > 
READ(2»'(ADMfUNHY 
REAl!<2.*)NY»« 

IF(H0D'NY«2!' .EG.OTHEN 
1ANY+NZ 
HY..2=NY/2 
ELSE 

I/-«Y+NZt1 
(U.2=«r,2kJ 
END IF 
DO 1*1. S 

READ ( 2 • < At > > DUhrtY 
END DO 
DO I=1»NXSI 

PEAD(2»*):;M:a.J).J=MZ) 

END DO 
CLOSE (2) 

PRINTtPVNJ.DAT !S FINISHED READING' 

C SET XS=XSI TO USE NORMAL VELOCITIES AT JACK LOCATIONS 
DO I=1<NXSI 
XS.LOCa)=XSI..LOC!I' 

END DO 


C Reading infaroation about noroal velocities IVNO.IATJ 
0PENC»NA«E*'VN0' .S!ATUS='0LD i 
READ(2* 'TAD DDllHHY 
RE AD ( 2 • t ) RLENGTH ? WIDTH 
READ(2< <22X.Fi.D')U-INF 
READ<2>' (AD DOWNY 
READ(2*'(24X»I5)') NX 
READ(2»'(AD' >DUHNY 
READ(2>I)(X-L0C( J) »J=1»NX) 

READ(2»'(34X»I3>’)NY 
READ (2* ' < 18X. Fb. 3 ) ' ) DELTA. Y 
READ(2.'<AD'»DUNNY 
READ(2*t)<Y_L0C(D*I s DNY/2) 
READ(2»'!34X*I3)'iNZ 
READ (2» '(18X»F6v3) ' )DELTA.Z 
IF(M0D(NY»2) .EO-O'THEN 
IZ=NY+NZ 
NY.2=NY/2 

ELSE 

IZ*NY+N7tl 
NY.2=NY/2tl 
END IF 


ORIGINAL page 13 
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DO 1*1 to 

READGr iAD i DUMMY 
END DO 
DO 1=1. NX 
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REAIi(2.*)(VN0(I»J)«J=M2) 

END DO 
CL0SE(2) 

PRINT!? 'VNO. DAT IS FINISHED READING 

C Calling subroutine to calculate the integrator! normal velocities 
C at each X - position specified? for ‘he upper and !ow<-r ft* If ’f 
C the wind tunnel. 

CALL SIMPSON! NXSI ? NY. NZ? DELTA- » . DELTA.,"- 

C Calling the subroutine to calculate tin? Point, load diselacanent matrices 
C using a summation over an infinite seriest truncating when t-.e 
C values have converged. 

CALL MATRIX { NXSI ? NX ? NY ? RLENGTH? WIDTH ) 

C The matrices have been calculated and the normal velocities have been 
C integrated? the normalized load at each va^ is calculated below using 
C the IMSL subroutine CLEQT2P] to solve the matrix equation fur X 

C A X = B 

C where A = the matrix of the point load displacement integrated wrt Y 
C X = the unknown normalized load vector 

C 8 = the normal velocities integrated over the perimeter using 

C Simpson's rule and spanwise using the trapazoiod rule 


N=1 

IER=3 

IA=40 

IB=5 

C Calling IMSL subroutine to solve for the normalized load vector for lower 
C half of tunnel 

CALL LEQT2F ( WSTARL ? N ? NXSI ? IA? SI MS.L » IF? WORKSPACE . IER) 

F’RINTI.'IER ='?IER ! print error flag to see if 

PRINT!? 'IB ='?IB ! calculations a^e o.k. 

IER=3 ! reset error flag for next calculation 

18=5 

C Calling IMSL subroutine to solve for the normalized load vector for the 
C upper half of tunnel 

CALL LEQT2F ( WSTARU? N? NXSI ? I A ? SIMS.U? IB ? WORKSPACE ? IER ) 

PRINT!? 'IER ='?IER ! print error flag to see if calculations are ok 
PRINT!? 'IB ='?IB 

C The variables SIMS.U()?SIMS_L() are now the normalized load vectors for the 
C jacks 
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OPEN ( 3 ? NAME= ' LOADS ' ? STATUS= ' NEW ' ) 

HRITE(3?!)' JACK FLOOR CEILING' 

WRITE(3?!)' LOCATION LOADS LOADS' 

DO I=1?NXSI 

WRITE(3?I)XSI_L0C(I)?SIHS-L(I)?SIHS.U(I) 

END DO 
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CLOSE ( 3 ) 

OPEN(3»NANE='DISPL' f STATUS='NEW' ) 

WRITE( 3»*> ' JACK FLOOR CEILING' 

WRITE ( 3» * ) ' LOCATION DISPLACEhENl DISPLACEMENT 

C Calculate the displacement VECTORS based on the loads 
DO I=1»NXSI 
DO J-liNXSI 

DISF‘1 L ( I )=DISFI L< I )+SIMS_L ( J 2 *W_DI3PL< I » J) 

H 1 3F L_ U < I ) I SPL_U ( I ) +SI MS.U * J ^ *w_DI SPL 1 1 1 J > 

END DO 

HR I TE : 3. * ) XSI.LOC ( I ) • DI SPL.L < I > . DI SPL.U < I> 

END DO 
CL0SE(3) 

C Calculate the differentiated matrix iwall slope) bs matrix multiplication 
OPEN ( 3 . NAHE= ' DUX_L ' < ST ATUS= ' MEM ' ’ 

OPEN! 4 » NAME* ' DUX.U ' .STATUS* 'NEW' > 

DO 1=1. NX 
DO J=1»NY_2 
DO K=1»NXSI 

WX.L<I»J>=WX_L<I« J)fW_0IFF<I*JpK>*3I«S_L<K) 
UXJja.J)=UX.U(I.J)+U-DIFFa»J.K)*SlNS-U>'K) 

END DO 
END DO 

HRITE(3.*hWX_L(I.j).J=l.NY_2) 

WRITE<4.»)(WX_U(If J)?J=1.NY.2' 

END DO 
CLOSE! 3) 

CLOSE! 4) 

C Calculates the residual normal velocities for each panel on the floor 
C and ceil ins of the wind tunnel 
DO 1=1. NX 

DO J=1.NY_2 

RESID.LU.J > =U. INFIUX.L (I.J) +VN0 ( I . J > 

JX=IZ+1-J 

RES I D_U ( I » J ) =U_ INFtWX.U ( I . J ) tVNO ( I > JX ) 

END DO 
END DO 

OPEN ( 2 . NAME* ' RES I DUAL.U ' . STATUS 2 ' NEW ' ) 

OPEN < 3 . NAME* ' RESI DUAL.L ' rSTATUS= 'NEW' ) 

C Write the residuals to file Cresidual.u.dat] and EresidualJ.dat] 

DO 1=1. NX 

WRITE(2.IMRESID_U(I>J). J=1?NY_2> 

URITE(3.t) (RESID.L(I»K).K=1>MY.2) 

END DO 
CLOSE (2) 

CL0SE(3) 

C Create data for FHIXZN.FOR prosram 
DO I=1.NY_2 
DO J=1 .NX 

PHINDAT(I.J)=RESID-L<J.I) 

END DO 
END DO 
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DO I=NY_2tl.NY.2)NZ 
DO J=1«NX 

PHINDAT ( I » J l=VNO< J» I ) 

END DO 
END DO 

DO I=NY2+N2+1 »NY+NZ 
DO J*1,NX 
K=IZ+i - I 

PHINDAni.J^RESID-UIJflO 
END DC 
END DO 
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C Write date to file i^t for»et to be resd by PHIXZM.FOR prodren 
OPEN ( 3 , NAME= ' PH 1 NWALL ' , STATUS* ' NEW ' ) 
HRITE<3»*MfPHINMTiI»J).[*ttIZ)*j*J'NX) 

CLOSE ( 3 > 


CALL 6ETCPLKNTIME2) 

TIME=(NTIME2-NTIME1)/10C, 

NMIN=INT< TIME/60.) 

NSEC 1 I NT ( AMOB ( T IME , 60 . » ) 

PRINT ' < A* I3*A,J3»A)'« 'CPU RUN TINE *'.NhI».' MINUTES' , NSEC * ' SEC 
10NDS' 

2000 STOP 'EXITING PROGRAM' 

END 


SUBROUTINE NATRI X ( NXSI > NX , NY , RLENGTH , UI DTH > 

C THIS SUBROUTINE CALCULATES THE MATRICES FOR T HE DISPLACEMENT, THE 
C DISPLACEMENT DIFFERENTIATED WITH RESPECT TO X. AND THE 3PANWISE 
C INTEGRATED DISPLACEMENT. THE CALCULATED VALUES ARE STORED 
C IN A COMMON BLOCK LABELED 'W-MATSS AND THE VARIABLES USED IN THE 
C CALCULATIONS ARE DEFINED BELOW l 

C THE FOLLOWING VARIABLES ARE USED IN THE PROGRAMS 
C X.LOCO = THE X - DISTANCE ALONG THE LENGTH OF THE WIND TUNNEL 

C Y.LOCO = THE Y - DISTANCE ALONG THE WIDTH OF THE TUNNEL (FROM CENTER) 

C NX = THE NUMBER OF X LOCATIONS TO BE USED FOR CALCULATION 

C NXSI = THE NUMBER OF XSI (JACK) LOCATIONS TO BE USED FOR CALCULATIONS 

C NY = THE NUMBER OF Y LOCATIONS TO BE USED FOR THE CALCULATION 

C XSI-LOCO = LOCATION OF THE JACKS ALONG THE LENGTH OF THE WIND TUNNEL 

C RLENGTH = TOTAL LENGTH OF THE WIND TUNNEL 

C WIDTH = TOTAL WIDTH OF THE WIND TUNNEL 

C W.DISP.L(») * ELEMENT IN THE 2-0 DISPLACEMENT MATRIX 

C W_DISP_U(») * ELEMENT IN THE 2-D DISPLACEMENT MATRIX 

C W-DIFF(,,) = ELEMENT IN THE 3-D DIFFERENTIATED MATRIX 

C W.INTEG(,) = ELEMENT IN THE 2-D INTEGRATED LOWER MATRIX 

C ALPHA.M = PARAMETER USED IN CALCULATING THE ABOVE 
C = M*PI*WIDTH/(RLENGTH*2.) 

C where: 

C M = INDEXING VALUE FOR THE INFINITE SUM 

C PI = 3.41592654 

C B = A RELATIVE Y - DISTANCE BETWEEN 0 and 1 
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INPLICIT REAL <A-H?0-Z) 

IMPLICIT INTEGER (I-N) 

COMMON/LOCAT/ X.LOC ( 100) ? Y.LOC ( 30 ) ? XSI .LOC (405? XS.LOC ! 40 ) 
COHMON/W.MAT S/U.D ISPI. ( 40 ? 40) ?U JIFF { 100*30 .405? VSTARL ( 40 ? 40 ) ? 
1 WSTARU(40?40) 

PARAMETER <PI=3. 141592454) 

CONST=PIIUIDTH/ (2 . IRLENGTH) 

OPEN ( 3 ? NAME 5 'OISP.MAT. STATUS 5 ' NEW' ) 


C Calculation of the displacement matrix follows below 
FACT0R=2.IPII!3 
DO I=1?NXSI 
DO J 5 1?NXSI 
DO M=1 ? 150 

ALPHA.H s H!CONST 

W1 5 ( ( 1 . TALPHA.MITANH ( ALPHA.M ) ) « T ANN < ALPHA.M ) -ALPHA.M ) * 

1 SIN(M*PI»XS!.L0C(J)/RLENGTM)»SIN(M*PI*XSLL0C<I)/F1EN6TH) 

2 /Mt*3 
W.0ISPL<I»J)=W_DISfin*J)4Ul 

C AN AVERAGING OF T)C LAST TERMS OF THE SERIES IS EMPLOYED FOR SMOOTHING 
T10=T9 
T9=T8 
T8=T7 
T7=T6 
T6=T5 
T5= T 4 
T4 5 '3 
T3 5 f2 
12-T1 

T1-W_DISPL(I»J) 

END DO 

100 W-DISPU I * J > 5 ( T1+T2+T3+T4+T5+T6+T7+TS+T9+T 10 ) /( 10 . IFACTQR ) 

END DO 

WRITE(3?IHW_DISPL(I?J);J s l?NXSI) 

END DO 
CL0SE(3) 

PRINT!? 'DISPLACEMENT MATRIX IS CALCULATED' 


page is 

C Calculation of the differentiated matrix follows below? OF POOR QUALITY 

FACT0R=RLENGTHI2. !PI!I2 
DO M?NX 
DO J=l?NY/2 
DO K 5 1?NXSI 

D 5 <UIDTH-2,fY.L0C<J))/WIDTH 
DO M=l?150 

ALPHA_M 5 M!CONST 

IF ( BIALPHA.M . 3T . 1 0 . ) THEN 

W2 5 ( i 1 . tALPHA.rff 2 . IT.LOC ( J ) /HID TH 5 *EXP< -ALPHA.M 

1 *2.*f J.OC(J)/WIDTH))tSIN(M!PIIxSI.LOC(K)/RLENGTH) 

2 tCOS \ h*P I IX.LOC i I ) /RLENGTH ) /MII2 
ELSE 

W2=( (1 ,+ALPHA.MtTANH' ALPHA.M) >*SINH<ALPHA.N«D)- 

1 ALPHA.MIB4C0SH ; ALFHA..M4D) )!SIN(M!Pl!XSI.LuC(f.i / 

2 RL£NGTH)!COS(M!PI!X.LCC( I >/RLENG7H),- iW«2l 

3 COSH (ALPHA.M)) 
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END IF 

W_DIFF(I;J;K)=W_BIFF(I»J«K>+W2 
END DO 

200 O I FF ( I » J t K ) =W_DI FF (I ; J ; K ) /FACTOR 

END DO 
END DO 
PRINT*; I 
END DO 

PRINTS 'DIFFERENTIATED MATRIX IS CALCULATED' 


ORIGINAL PAGE IS 

QE POOR QUALITY 


C Calculation o? spanwise intearated aatm beams below 
OPEN ( 3 ; NAME= ' USTAR.MAT ' , STATUS* ' NEW ' > 

FACT0R=RLENGTH/<PI**4i 
DO I=1;NXSI 
DO J=1;NXSI 

USTARLGt J)-0. 

DO rt= 1 * 1 20 

ALPHA. M r M*C0NST 

«3--SIN'M*PI*XSI.L0C(J)/RLENGTH)*SIN(N*PI*XS.L0C(I)/ 

1 RLENGTH : *< 2. - ( 2. tALFHA_M*TANH< ALPHA.N) ) / 

2 C0SH(ALPNA.N)t/'W4 

US T ARL < I ; J ) =aST ARl < I • J ; +W3 

T10=T9 

T9=T8 

T8-T7 

T7=T6 

T6=T5 

T5-T4 

T4=T3 

T3-T2 

T2-71 

Tl=N3fARL(I»J) 

END DO 

WSTARL(LJ)=(Tl+r?tTi+T4+T5+T6+T7F;P+T?PT10)-10,*FACT0R 
WST ARU ( I f J ) -UST ARL 1 1 ;J ! 

END DO 

URITE(3;*)(USTARL(I;J).jM;NXSI> 

END DO 

CM Debusaina Diaanostics 
DO 1=1 f NXSI 

WRITE(3f*)I>XSI.L0C(I)»XS_L0C(I) 

END DO 
CLOSE! 3) 

PRINT*; 'USTAR MATRIX IS CALCULATED' 

RETURN 

END 


SUBROUTINE SIMPSONINXSI . NY . NZ - DELTA. Y» DELTA.Z) 


IMPLICIT REAL «A-H.0-Z) 
IMPLICIT INTEGER (I-N> 
INTEGER START 


COMMON/SIMP J N T /VN J < 40 ; 60 ) ; S IMS.U 140); SIMS.L ( 40 > 
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COHMON/LOCAT /X.LOC ( 1 00 i . Y J.OC (30'.. XSI.LOC ( 40 ) . XS_LOC< 40 ) 

DIMENSION VINTXU<40),V!NTXLUO) 

C CHECK TO SEE DELTA-2 AND DELTA. Y ARE I HE SAME. IF THEY ARE. USE SIMPSCN'3 
C INTEGRATION AROUND IZ. IF NOT. INTEGRATE EACH SEGMENT OF NY AND NZ SEPERATLY 
IFIABSYDELTA.Y-DELTA.Z) .GT. .001 1G0T0 10 0 


mmmmmmtnmmmmmmmmmmmmnmttmmm 

C THIS PORTION IS FOR EQUAL PANEL WIDTHS IN THE Y AND Z DIRECTIONS 
C CHECK TO SEE IF NY IS ODD OR EVEN FOR THE APPROPRIATE END CORRECTIONS 
IF (M0D(MY»2) .NE.O) THEN 
NYFLAG*1 
IZ=NY+NZH 
IZ.2-IZ/2 
ELSE 

NYFLAG=0 
IZ=NY+NI 
IZ_2=IZ/2 
END IF 

C SIMPSON'S INTEGRATION ROUTINE FOR DELTA.Z * DELTA.Y 
DO KX s l»NXSI ' EACH X LOCATION 
STARTS 

C CHECK FOR ODD OR EVEN INTERVAL. IF ODD. APPLY SIMPSON'S 3/8 RULE 
C FOR THE FIRST FOUR POINTS 

IF(M0D(IZ_2.2).E0.0)THEN 

SIM$.L(KX> s 3.<DELTA. Z/8.KVM J(KX*1>+3»IVNJ(KX«2)+3.I 
1 VNJ(KX.3)+VNJ<KX.4)> 

STARTS 
END IF 

C START OF SIMPSON'S INTEGRATION FDR EVEN INTERVALS 
IF(M0D(START»2).NE.0)THEN 
IE0FLA6-1 
ELSE 

IE0FLAG=0 
END IF 

DO KZ*START»IZ_2-1 

IF(MOD(KZ>2).EO.O)THEN 

IF(IE0FLAG.EQ.1)THEN 

SUM=SUMF2.*VNJ<KX»KZ) 

ELSE 

SUM=SUM+4 . tVN J < KX> KZ ) 

END IF 

ELSE 

IF (IE0FLAG.EQ.il THEN 
SUM=SUMM.*VNJ(KX»KZ) 

ELSE 

SUH-SU«tMVNJ?KX>KZ) 

END IF 
END IF 
END DO 


SIMS.L{KX)=SIMS-L;r./:t(SnmvNJ'KY.S’ART-l>fVNJ(KX.IZ_2» 
1 IDELTA.Z/3. 

SUM=0. 

end-VN=(VNJ(KX.iz-2)+vnj(kx.;’:h;;/2. 
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C CHECK FOR WHICH END CORRECTION TO USE 
IF i NYFLAG . EQ . 1 ) THEN 

ENDCOR'R=DEL TA.Z/2 . > ( HIN ( UN J ( KX .< 12-2.' .• £NI ! _ *K ! iwBS( 
1 VN J< KX » IZ_2 ) -END _VN ) /2 ► > 

ELSE 

ENDCORR- DEL T A.I < 2 . * ( UN J > K 1 . 1 > +M I N ( VN J ( KX « 1 1 . 2 1 • 

1 END.VN ) TABS ! VN J(KX » IZ..2 : -END.VN ) H , > 

END IF 

C INTEGRATED VALUE FOR THE uOWER HALF OF TUNNEL 
SIMS.L (KX ) =-2 ,*( 3M3-L< KX) +ENDCCRR; 

C START SIMPSON’S ROUTINE FOR THE UPPER HALF OF TUNNEL 
SThRT=IZ.2t2 


C CHECK FOR ODD OR EVEN INTERVAL 

IF ( HOIK IZ.2 » 2 ) . EQ .0) THEN 

S IMS.U < KX ) =3 . IOEL TA. Z/8 . * < VN ... < KX - ! Z.2 1 • • +3 . * 

1 VNJ(KXi IZ-2K2 F3»*VHJ(KXjIZ_2+3)+VNvUKXf IZ_IFO) 
START=IZ_2+5 
END IF 


C START THE INTEGRATION FOR EVEN INTERVAL 
IF(M0D(START»2) .NE.OiTHEN 
IE0FLAG=1 
ELSE 

IE0FLA6=O 
END IF 

DO KZ=START<IZ-1 

IF(rt0B(KZ»2).EG.0)THEN 

1FUE0FLAG.EQ.DTHEN 

3Ui1=SUHF2.*VNJ(KX.KZ) 

ELSE 

SI!N=8UN P 4 , f VN J ( K X » KZ ) 

END IF 

ELSE 

IFdEDFLAG.EO.l 'THEN 
SUM=SUt1+4 . KVN J< KX «KZ 5 

ELSE 

SUM=SUM+2. *VNJ(KX»KZ > 

END IF 
END IF 
END DO 

SIHS_U<KX)=SIMS_U(KX)+DELTA.Zt(VNJ KX*S T ART-l)tSUM+ 

1 VNJ(KX»IZ))/3. 

SUH=0. 

END.VN=(VNJ(KX.IZ.2)+VNJ<KX.IZ-ZF1 • 2. 

C CHECK TO SEE WHICH END CORRECTION APPLIES 
IF < NYFLAG . EO » 1 ) THEN 

ENDC0RR=DELTA_Z/2 . t ( HIN ( VNJ '!_?+; t >FND.vNi+ 

1 ABS(VNJ(KX . IZ.2FI )-EMD..vN)/2 . ) 

ELSE 

ENDCORR-DELTA.Z/2 , t ( VN J < KX > I Z ) + MI N ( W J ? K U Z . . M ' . 
1 END_VNHABS<VNJ(KX» IZ.2F15 -END.VHi .2 * > 

END IF 

C INTEGRATED VALUE FOR THE UPPER HALF OF THE TUNNEL 
SIHS_U(KX)=-2.t( S IHS.U ( KX ) FENDCORR ) 

END DO 


ORIGINAL PAG£ i2>. 
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GOTO 1000 


! CALCULATIONS ARE COMPLETE. WRITE T 0 FILE 
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ttmttmttmtmmmtmtmmnmmmmtnmtmmmtmttt 

C THIS SEGMENT FOR PANELS OF UNEQUAL WIDTHS IN THE T AND ? DIRECTION 
C SIMPSON'S ROUTINE FOR DELTA.Y .NE. DELTA. 2 

100 IF(M0D(NY»2) »EQ.0)THEN 
NYFLAG=0 
NY.2=NY/2 

ELSE 

NYFLAG=1 
NY.2=NY/2H 
END IF 

NZ.2=NZ/2 

C SIMPSON'S ROUTINE FOR DELTA.2 .NE. DELTA.Y 
IE0FLAG=0 
DO KX*1»NXSI 
STARTS 

IF(MOD(NY.2»2>.EO.O>THEN 

SIMS.LY*3,>D£LTA.Y/B.t(MNJ(KX»l)F3.tVNJ(KX*2)+3.l 
1 UNJ(KX»3)FWJ<KX»4) ) 

START=5 
END IF 

IF(M0D(START>2) .NE.O) THEN 
IE0FLAGM 
ELSE 

IE0FLAG*O 
END IF 

DO KZ=START»NY_2-1 

IF(M0D(KZ»2).£0.0)THEN 
IFIIEOFLAG.EQ.WHEN 
SUM=SUH+2 . tVN J < M » KZ ) 

ELSE 

SUM s SUN+4 . <UN J ( KX * KZ ) 

END IF 
ELSE 

IF(IEOFLAG.E0.1)TH£N 
sUM*SUM+4 . tVNJt KX »KZ> 

ELSE 

sum»sum+2.*vnj:kx»kz> 

END IF 
END IF 
END DO 

SIMS-LY*SINS.LY+(SUM+VNJU^$TARM 'tVmx.«Y.2) > 

1 IDELTA.Y/3. 

SUM=0. 

END. VN S ( VN J ( KX » NY.2 ) +MN J I KX » N v . 2+ 1 • 2. 

IF(NYFLAG.EQ,1)THEN 

ENDCORR=DELTA_ Y/2 . * ( MIN ( UN J< KX • YY .2 > - EW.VN ) f 
1 ABS(VNJ<KX»NY_2)-END.W/2.) 

ELSE 

ENDCORR=DELT A. Y/2 . > ( UN J ( KX » 1 ) FM IN W . 

1 END.VN ) +ABS ( VN J ( KX » N Y.2 ) -END.MN* H , 

END IF 

SIHS.LY=SIMS_LY+ENBCORR 
START=NY_2+2 

IF ( MODI NZ.2 * 2 ) . EO . 0 > THEN 


ORIGINAL PAGE IS 
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S IMS-LZ-3 . MELTA.Z/8 . * < UN J < KX , NY.2F1 ) +3 . * 

1 UN J < KX » NY _ 232 Hi . *UNJ < KX , NY .2+3 ) FUN J< KX - NY-234 ) < 

IE0FLAG=1 
START=NY_2t5 
END IF 

IF (MOD ( STAR r • 2 ) . N£ . 0 ) THEN 
IE0FLA6=1 

ELSE US £OOR QUALITY 

IE0FLAG=0 

END IF 

DO KZ=START.NY_?3NZ.2-1 
IF<H0D<KZ*2».EQ.0»THEN 
IF (IEOFlAG.Efl.il THEN 
SUh=SU«32.*UNJ(KX»KZ) 

ELSE 

' •LiH=SUN+ 4 . *UN J . K X » K Z ) 

END IF 

ELSE 

IF«IEOFLAP.EB.l»THEN 

3U«=SUH+4.*UNJ<KX.KZ) 

ELSE 

SON=SUM32.n'?M<KXfKZt 

END IF 
END IF 
END DO 

SIMS_LZ=SI«S_LZ+ ( SUN+VN J ( KX . 3 1 AF r- * HUNJ (KX « NY.2+NZ.2 ) ) 

1 WELTA-Z/3. 

SUH=0. 

END_VN1=(UNJ<KX»NY_2)3UNJ(KX<NY. 2h i >< 2. 

END. VN2= ( VN J ( KX » NY.2+NZ.2) FUN J< KX .«Y .? 3NZ .2+1 ) i n . 
ENDC0RRl=DELTA_Z/2.l(NIN(UNJ(KXiNY.23t Hi-ND.UNtH 
1 ABS< UN J ( KX f NY.2+ 1 ) -END.ON1 ) H . ) 

ENDC0RR2=DELTA.Z/2 . *<HIN(UNJ(KXiNY_23NZ_2) » END.VN2 ) + 

1 ABS(VNJ(KXi NY.2FNZ.2 ) -END.UN2 )/2.) 

SIHS-L(KX)=-2.*(SINS-LY+SI«S_LZ+£NDCORR14ENDCORR2) 

SIHS_LY=0. 

SIMS_LZ=0. 

START=NY_2+NZ_232 

IF<NOD(NZ_2f2).£Q.O)THEN 

SIHS_UZ=3 . *DELTA_Z/8 . t (UN J( KX* NY.2+NZ.2+1 ) 33 . * 

1 UN J ( KX * NY.2+NZ.2+2 ) +3 . tUN J ( KX * NY.23NZ.233 H 

2 UN J ( KX i NY.2+NZ.2+4 ) ) 

START=NY_2+NZ-2+5 

END IF 


IF (NOD (START * 2 ) . NE . 0 ) THEN 
IEOFLAG=l 

ELSE 

IE0FLAG=0 
END IF 

DO KZ=START*NY_2+NZ-1 

IF(NOD(KZ»2).EQ.O)THEN 
IF(IEOFLAG.EQ.l>THEN 
SUN=SUN+ 2 . <UN J ( KX • K Z ) 

ELSE 


SUN=SUH+4 . FUN J ( K > > K Z 
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END IF 

ELSE 

IF ( IEGFLA6 . EG . 1 / THEN 
SUM= SURF 4 . IVN J \ KX i KZ) 

ELSE 

SUH=SUHF2 » IVN J ( KX > K2 ) 
END IF 
END IF 
END DO 


SIMS_UZ=SIMS_U2r < SUM+MN J < KX • 37 AF:T- 1 )FVNJ(KX?NY.2FNZ) ) 

1 KDELTA.Z/3. 

SUN=0 . 

END- VN2= < v'N J ( KX » NY.2FNZ ) F VN J ( KX « N Y_ 2 FNZ H > ■ ,' 2 . 
END_VN1=(VNJ<KX»NY_2+NZ_2>+VNJ(KX.ny_24«2.2U ■ ' 2. 
ENDC0RR2-0EL T A. Z 7 2 . * ( NIN ( VNJ ( KX ? NY .2 FNZ > * END 'M2 ■ F 
1 ABS< VNJ< KX »NY_2+NZ ' -END.VN2 )/2. ) 

ENDCORR1 =DELTA.Z/2. * ( NIN ( VNJ< KX .NY . 2FNZ.2F1 5 . END.VN1 • F 
1 ABS ( VN J ( KX »NY_2?NZ_2F1 ) -END .VNl>/2 . ) 

SINS-UZ=SIHS-UZfENDC0RRlFENDCCRR2 


START=NY_2FNZF2 

IF(M0D(NY_2>2) »E8.0)THEN 

SIHS-UY=3.*0ELTA_Y/8.*<VNJ<KX.NY_2FNIF1)F3.* 

1 VN J < KX » NY.2+NZ+2) F3 . *W J < KX , NY-2FHZF3 ) F 

2 VNJ(KX»NY_2FNZF4); 

START=NY_2FNZF5 

END IF 


IF(H0D(START»2).NE.0)THEN 

IEOFLAGM 

ELSE 

[E0FLA6=0 
END IF 

DO KZ=STARMZ-l 

IF<HOD<KZ»2>.EO.O)THEN 

IFUEOFLAG.EQ.DTHEN 

SUH=8UhF2.tVNJ(KX»KZ) 

ELSE 

SUN=SUNF4 . >VNJ( KX>KZ / 
END IF 

ELSE 

IF(IE0FLAG.E0.1)THEN 

SUH=SUNF4,«VNJ;KXfKZ‘. 

ELSE 

SUN=SUNF2 . 4VN J < KX « KZ ; 
END IF 
END IF 
END DO 


ORIGINAL FAG 
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SINS.UY=SIHS_UYF< SUNFVN . ! K X . START- 1 > FVNJi KX. IZ) ) 

1 IDELTA-Y/3. 

SUN=0. 

END. VN= ( VN J ( KX r NY.2FNZ > t VNJ KXiNY 2+NZ+l > • ‘2. 

IF ( NYFLAG • EQ • 1 ) THEN 

ENDC0RR=DELTA_Y/2 . * ( NIN ( VNJt KX • a Y .2 tNZM » • END -VN ) F 
1 ABS(VNJ(KX»NY-2FNZFl)-END-VN)/2. < 

ELSE 

ENBCDRR=DELTA. Y /2 . * < VN J< KX » IZ> FHIN( VNJfKX . 

1 NY.2FNZF1 ) »END.VM)FAPS< W J(KX>NY.2FNZF1 J-END.'JN) »2 • ) 
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END IF 

SIHS.U ( KX ) =-2 . * ( S I MSJJZFSIMS -UY f ENDCORR ) 

SIMS_UZ=0. 

SIMS_UY=0. 

END DO 

1 OOO OPEN ( 3 , NAME* ' 51 MSINTEG ' , STATUS* ' NEW ' ) 

WRITE ( 3 f * ) NXSI 
DO 1=1. NXSI 

WRITE<3.*)XSI-L0C(I)»3IMS-L(I)»3INS.U(I * 

END DO 

mmmmmtnmmmmmmmmmnm 

C COMPUTE THE STREAMWISE INTEGRAL OF THESE SPAN INTEGRATED VELOCITIES 
C AS A FUNCTION OF X 

C THIS IS USED SOLUTION FOR THE JACK LOADS BY USE OF THE US TAR 
C (SPAN INTEGRATED DISPLACEMENT) MATRIX 


VINTXL(l) = 0.5*SIMS-L(l)t XS-LOC(l) 

VINTXU(l) = 0.5*SIMSJJ(l)t XS.LOCU) 

C THIS IMPLICITLY ASSUMES THAT SIMS(X=0) * 0 

DO I =2. NXSI 

VINTXL < I )=VINTXL ! I -1 ) + 0.5NS1MS.LU )«IHS.L(I*1) ) * 
t (XS-LOC(I)-XS.LOC(I-l)) 

VINTXUl I )=VINTXU( 1-1 ) + 0.5t(SIHS_U(IHSWS_U<I-H) t 
* ( XS.LOC ( I ) -XS-LOC (I-l)i 

END DO 

C DUMP STREAHWISE INTEGRATED VELOCITIES BACK INTO SIMS. ARRAYS 

DO 1=1. NXSI 

SIMS.U(I)=VINTXU(I) 

SIMS_L(I)=VINTXL(I) 

END DO 

WRITE<3»3000> 

3000 FORMAT (X. 'FOLLOWING ARE THE STREAMWISE INTEGRATED VELOCITIES') 

WRITE(3.*)NXSI 
DO 1=1. NXSI 

WRITE(3»*)XSI_L0C(I)»XS_L0C(I)>SIHS_L(I)»SIHS-U(I) 

END DO 

CLOSE (3) 

PRINT*. 'SIMPSON' 'S INTEGRATION COMPLETE' 

RETURN 

END 
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PROGRAH AFHODLJ 

COMPUTES THE WALL- FREE VELOCITIES! VN> NORMAL TO THE 
WALLS OF A SQUARE WIND TUNNEL OF AN AEDC MODEL 
10 PE USED IN A LANGLEY EXPERIMENTAL TEST. 

IT ALSO COMPUTES DATA AT THE 0.3 METER TUNNEL JACKS AND 
WRITES THE FILES 'VNO.PAT' AND 'VNJ.BAT' FOR JACOISPL. 

THIS PROGRAM USES 2INX+1 PLUS DX PANELS OVER THE LENGTH OF 
THE WIND TUNNEL AND NYFNZ PANELS AROUND THE HALF 
PERIMETER. TO INCREASE THE NUMBER OF PANELS 
NEW VALUES MAY BE NEEDED IN THE DIMENSION STAIEHENT 
AND r«E LIMITS NX. DX AND NZ - 

r HIS PROGRAM REPRESENTS THE MODEL BY TWO SOURCES(BODY) 
v’EE LINE-V0RT1CES(W1NGITAIL LIFT) AND +1- LINE SOURCES 
'UINGSTAIL THICKNESS). BODY LIFT BY LINE DOUBLETS 
IT ASSUMES THE CP IS LOCATED AT THE ZERO POSITION OF THE 
X COORDINATES 

HE WING IS ASSUMED TO CARRY 80% OF THE TOTAL LIFT 

DIMENSION X(-80:i20).XJ(-20:0).Y(80)»Z<80).VNS<80.-80:i20>» 
JNA(80.-80!120).VN<80.-80'120).DZXL(~80!120)iDZXU(~80»120). 
ZL(-80:i20).ZU(-80:i20).VNN!80.-80U20>.VNJ(80.-80tl20). 
«.LOC(-80!120) 

XJ(IX) ARE THE 1/3 METER TUNNEL JACK LOCATIONS 
XUO.YUZ) AND Z!IZ> ARE THE CENTERS OF THE PANELS 

THE TRACOR ALGORITHM IS APPLIED IN THIS PR06RAH 
THE NUMBERS BELOW ARE FIXED BY THE DIHENSIONS OF THE 
AEDC MODEL 


PARAMETER <PI=3. 14159, S*0. 10668* SV=. 0837863, AR-3. 5, 

D=0 . 03556 » C= . 060 1 98 » CT=0 .03556 » «T= . 0072238 »WTT= . 004267 » 

ST=. 05588, SVT=. 043888, A=. 33. ART=3.1429) 

REAL MACH 
TYPE 10 

FORMAT <10X» 'ENTER CL MACH NX HZ') 

REAB(5,t) CL. MACH. NX. NZ 
BETA-SORT! 1-MACHIMACH) 

THE FOLLOWING RESULT FROM THE MODEL DIHENSIONS 

XRVC--. 03759 ! LOCATION OF VORTEX tC/L INTERSECTION(NETERS) 

XRV=XRVC/BETA 

XRVTC=0.0889 'SAME FOR THE TAIL 
XRVT=XRVTC/BETA 

XSL=-. 10795 ! LOCATION OF FOREWARD SOURCE 

XST=0.1524 ! METERS. SAME FOR THE REAR SOURCE 

IG S . 57735/BETA ! TANGENT OF THE STRETCHED SWEEPAN6LE 
REAL SWEEP IS 30 DEGREES! 

TS=lfT6*TG 


BX=0 

NZ2=NZ/2 

NY-NZ 

NY2 s NY/2 


ORIGINAL PAGE iS 
OF POOR QUALITY 



ORIGINAL page is 

OE POOR QUALITY 

NOW SET UP the COORDINATES OF THE POINTS AT WHICH 
VELOCITIES ARE TO BE CALCULATED. 


DO 12=1 iN22 
Z(IZ'=A/2 

Y(lZ>=A/NZ4<IZ-0.5) 

END DO 

DO : Z=NZ2+ UNY2+NZ 

ran - a /2 

.:H2) A/NZ4(NZF0.5-IZ> 

END no 

DO IZ=NZ+NY2+1»NZ+NY 
2<IZ,=*A/2 

T( IZ'=(NZfNY+0.5-IZ)IA/NZ 

END DO 

20 IX=-NX»NX+DX 

*(IXi=A/NZtIX ! FOR PANEL CENTERS RELATIVE TO THE CP POINT 
END DO 

X.J< -20 >=-23.251.0254 

XJ<-19)=-17, 54.0254 

<J(-18)=-12, 54.0254 

XJ<-17)=-8,54,0254 

XJ(-16)=-5.54,0254 

XJ<-15)=-3,54.0254 

XJ(-14>=-2.04,0254 

XJ<-13)=-0.5*,0254 

XJ(-12)=1.04.0254 

XJ(-11)=2. 51.0254 

XJ(-10)=4,OI.O254 

XJ(-9)=5,54.0254 

*J(-8)=7,54.0254 

XJ(-7)=9,5*.0254 

*J(-6)=11.5I.0254 

XJ<-5>=14,54.0254 

XJ(-4)=18. 54.0254 

XJ(-3)=23, 54.0254 

XJ< -2 >=28. 54. 0254 

XJ<-1)=33.54.0254 

XJ(0)=:9.04.0254 


C I RC=0 . 84CL4S4S/ ( SV4AR ) 

ALPHA=0 . 84CL/ (3.2+1. 7554MACH443 ) !IN RADIANS 

THIS ASSUMES AN APPROXIMATE MACH BEPENCENCE FUR CL/ALPHA 


Q=WT 

QT=WTT 

C I RCT= . 24CL4ST 1ST / ( SVT 4 ART ) 'TAIL CIRCULATION 

BDD= 1+ . 54D4D4BETA4BET A/ ( < XST-XSL ) 4 < XST-XSL ) ) IAPPROX STRENGTH CORRECTOR 

DO IX=-NX»HX+DX 

STRETCHED DISTANCES FROM FIELD POINT TO INTERSECTION OF C/L AND! 

XV= ( XRVC-X ( IX ) ) /BETA ! WING VORTEX 

XL= ( XRVC-C/4-X ( IX) ) /BETA ' WING LEADING EDGE 

XT=(XRVCP34C/4-X(IX))/BETA ' WING TRAILING EDGE 
XLS= OCSL- < ( I X ) ) /BETA ! BODY SOURCE 

XTS=< XST-X < IX) > /BETA ! BODY SINMAFT) 

XVT=(XRVTC-X(IX>)/BETA ! TAIL VORTEX 

XLT=XVT-CT/4/BETA ' TAIL LEADING EDGE 



<T'=XVT«*CT/4/BETA 


< TAIL TRAILING EDGE 


RH IS THE ARM LENGTH FROM THE POINT OF ROTATION OF 
THE BALANCE/STING SYSTEM TO THE CP ORIGIN OF COORDS. 
THE NOBEL IS DISPLACED UPWARDS BY A DISTANCE ALPHABRH 

RM-O. 50126? ! METERS 


DO 12=1. «<2 
F=<mr<IZ5*TG) 

FJ=<H'-Y(IZ)>TG> 

FT«(XVTFY.I2)ITG) 

FT1=(XUT-Y(IZ)<TG) 

*r =Y(fZ)*»2 

ZZ= l Z U Z ) ♦ ALPHAtRH ) 112 
01=4«(FIF+TStZZ) 

Q2=4*(FMFHTSBZZ) 

:3=4*<<.<LWIZ)»TG)tt2+TS*ZZ) 

Q4=4*<(XL-Y(IZ>*TG)t*2+TS*ZZ) 

Q5=4*< (XT+Y(IZ)*TG)«2+TSIZZ) 

96=41 ( ( XT- Y ( IZ ) *TG ) W2+TSIZZ ) 

01T=4»(FT*FT+TS*ZZ) 

Q2T=4*(F1T*F1T+TS*ZZ) 

Q3T=4*<<XLW<IZ)tTG)**2+TSIZZ) 

Q4T=4*( ( XLT- Y ( IZ ) «TG) **2+TS*ZZ) 

G5T=4I ( ( XTT+ Y ( IZ ) »TG ) W2+TSIZZ ) 

06T=4*(<XTT-Y(IZ)*TG)«2+TS*ZZ) 

R1=S8RT ( XVW2+ZZ+Y Y ) 

R2=SQRT<(XV+TG*SV)**2+ZZ+<Y<IZ)-SV)«2) 

R21=S0RT((XV+TG*SV)«2+ZZ+(Y(IZ)+SV)**2) 

R3=5GRT< <XLFTGID/2)»*2+ZZ+<Y( IZ)-D/2)»I2) 

R31=S0RT( ( XL+ TGBD/2 ) W2FZZ+ ( Y ( IZ ) +D/2) M2 ) 

R4=SQRT((XL+TG*S)**2+ZZHY(IZ)-S)«2) ORIGINAL PAGE IS 

R4 1 =SORT ( ( XL+TGIS) M2+ZZ+ ( Y( IZ)+S) t<2 ) OF POOR QUALITY 

R5=S0RT ( ( XT +TGID/2 ) W2+ZZ+ ( Y ( IZ)-D/2 ) M2 ) 
R51=SORT<<XT+TG*D/2)M2+ZZ+(YUZ)+D/2)M2) 

R6=SQRT<(XWGIS)M2+ZZW(IZ)-S)M2) 

Ro 1=S9RT ( ( XT FTGtS ) M2+ZZ+ ( Y ( IZ ) tS) M2) 

R7=ZZW4XLS*XLS 

F;8=ZZFYY+XTSmS 

R1 T=SQRT ( XVTM2FZZF YY ) 

R2T=SQRT<<XVWG*SVT)M2+ZZ+(Y(IZ)-SVT)M2) 

R21T=S0RT( <XVWG*SVT)M2+ZZHY( IZHSVTJM2) 

R3T =SQRT ( ( XLTFTGID/2 ) M2+ZZ+ ( Y ( I Z ) -D/2 ) >12 ) 

R31T=SQRT ( ( XLTfTGtD/2 ) M2+ZZ4 ( Y ( IZ HD/2 ) 1*2 ) 

R4 T -SORT ( ( XLT+TGIST ) M2+ZZT ( Y ( IZ ) -ST ) M2) 

R41T=SQRT( ( XLT+TGIST ) M2FZZH Y l IZ ) +ST) l»2 ) 

R5T=S0RT < ( XTTf TGIB/2 ) M2+ZZ4 ( Y ( I Z )-D/2 ) M2 ) 
R51T=SQRT<<XTWG*D/2)M2+ZZHY<IZHD/2)M2> 

R6T=S8RT( (XTWG«ST)M2+ZZ+< Y< IZ)-ST)M2) 

R61 T =S0RT ( < XTT +TGIST ) M2FZZ+ ( Y ( I Z ) +ST ) M2 ) 


VNS< I Z • I X ) =-0*D*DDD* ( Z < I Z > + ALPHA*RM ) /16/BETA* < 1/ ( R7*SQRT ( R7 ) ) 
-0.61/<R8>SQRT(R8) > ) 

THESE ABOVE ARE THE BODY SOURCE TERNS 

VNS(lZrIX)=-Q*(Z<IZ)+ALPHA*RN)/PI/BETAK(TSIS+XL*TG-Y<IZ))/(03*R4)- 

<TS*D/2+XUTG-Y(IZ))/(Q3IR3)+(TStS«L*TG+Y(IZ))/(04«R4i)- 

(TS*D/2+XLITG+Y<IZ))/(Q4*R31)-(TS$Sm*TG-Y(IZ))/(05«6)4 



<TS*D/2+XTITG-Y(IZ))/(Q5IR5)-(TSIS+XT*TG+Y(I2))/(Q6*Wl)+ 

(T3tD/2W*TGfY(IZ))/<Q6*R51))+VNS(IZ»IX) 

THESE WERE THE WING LINE SOURCE TERNS 

VNS!IZ>IX)--QT*(Z(IZ) +ALPHA*RN ) /PI /BETA* ( ( TS*ST4XLT*TG- Y ( I Z ) ) / 
<G3T*fi4T)-iTS*D/2+XLT*TG-Y(IZ))/(Q3T*R3m<TS*ST+XLT*TG+Y<IZ))/ 
<Q4T*R4!T)~(TS*D/2+XLT*TG+Y(IZ))/(Q4T*R31T)-(TS*ST*XTT*TG- 
Y - 1 2 ) : / ( 05T*R6T ) + < TS*D/2+XTT*TG-Y ( IZ» ) / (Q5T*R5T ) - 
< TS*ST+XTT*TG+Y <12 ) ) / ( Q6T*R61 T)+< TS*B/2+XTT*TG+Y< IZ) ) / 
UT6T*R51T)) + UNS(IZ.IX) 

THESE WERE THE TAIL LINE SOURCE TERNS 

VNA < I Z < I X > =fl*B* ALPHA/8* < ( Y Y-ZZ ) / < Y Y+ZZ ) **2t ( 1 -XLS/SQRT (R7) ) 
+ZZ*X!.S/( YYFZZ)/(R7*SQRT(R7))-0.61*((YY-ZZ)/(YYEZZ)*I2* 

■■ 1 - X T S/SQRT ( R8 ) ) +ZZ*XTS/ ( YY+ZZ)/(R8*S0RT(R8) ) ) ) 

THESE ARE THE DOUBLET (BODY LIFT) TERNS 


TT=4 / 3 1 * i ( TStSVTX V* T6 - Y ( IZ ) )/R2- ( XV*TG- Y ( IZ) )/Rl ) 

TC=4/02* ( ( T S*SV+XM*TG+Y ( IZ ) )/R21 - ( XV*TGtY (IZ) )/Rl) 

VNA( IZ ■ IX)=VNA< IZ» IX )fCIRC/4/PI*(F*TT+Fl*T0+(Y(IZ)-SV)/ 
iZZ+< YUZ)-SV>**2)t(i-<XV+SV*TG)/R2)-(Y<IZ)+SM)/(ZZ+ 

( Y < I Z > +SM) *12) < ( 1 - ( XV+SV*TG) /R21 ) ) 

THE ABOVE ARE THE WING LINE VORTEX TERNS 

ITT=4/01T*((TS*SVT«VT*TG-Y(IZ))/R2T-(XVT*TB-Y(IZ))/R1T) 
T0T=4/22T*< ( TS*SVT+XVT*TG+Y(IZ) )/R21T-(XVT*TG+Y(IZ) )/RlT) 

VNA( IZ) IX >-VNA{IZ>IX)+CIRCT/4/PI*(FTtTTTTFlT*TOT+(Y(IZ)-SVT)/ 
;ZZ+< Y( IZ)-SVT'**2)*(1-(XVT+SVT*TG)/R2T)-<Y(IZ)4SVT)/ 

*. ZZ+ •; Y ( IZ)+SVT)**2) *< l-< XVT+SVT*T6 ) /R21 T ) ) 

THE ABOVE ARE THE TAIL VORTEX TERNS 

VN(IZ* !X)=VNA(IZi IX) IVNS<IZ» IX) 

£i)D HO 


00 IZ=NY2+NZ+1>NY+NZ 'UPPER HALL 
F-(XV+Y(IZ)*TG) 

FldXV-YXIZ)*TG) 

FT=(XVT+Y(IZ)*TG) 

FT1-(XVT-Y(IZ)*TG) 

YY=YdZ)**2 

ZZ= ( Z ( I Z ) + ALPHA4RH ) t»2 
Ql=4*(F*FfFS*ZZ) 

02=4*(F1»F1FTS*ZZ) 

Q3=4*( (XL+Y(IZ)*TG)**2+TS*ZZ) 
P4=4*i(XL-Y(IZ)*TG)**2+TS*ZZ) 

05=4* < < XT+Y (IZ ) *TG ) **2+TS*ZZ ) 

Ge=4*( (XT-Y( IZ)*TG)**2+TS*ZZ) 
Q1TM*(FT*FT4TS*ZZ) 

G:2T-4*(FlT*FlTtTS*ZZ) 

Q3T'4*( <XLT fY < IZ)*TG)**2+TS*ZZ) 
04T=4*(<XLT-Y(IZ)*TG)**2+TS*ZZ) 

05T -*4* < XTT+ Y < IZ ) *TG ) **2+TS*ZZ ) 

06T ; 4*( ( XTT-Y ( IZ)*TG)**2+TS*ZZ) 
Rl-SQRT(XV**2+ZZ+YY) 

R2=SQRT < ( XV+TG*SV ) **2+ZZ+( Y ( IZ)-SV )**2 ) 
R21=S0RT < ( XV+TGISV ) **2+ZZ+ ( Y (I Z ) +SV ) **2 ) 
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R3=S0RT ( < XL+TG*D/2 ) **2+ZZ+ ( Y( IZ )-D/2) *12) 

F;3 1 =SGRT < i XL+TG*D/2 ) **2+ZZ+ ( Y ( IZ ) +D/2 ) t*2 ) 
R4=SQRT ( ( XL f TG*S ) **2+ZZ+ ( Y ( IZ )-S) **2) 
R4l=SQRT ( ( XL+TG*S ) **2+ZZ+( Y ( IZ )+S) **2) 
R5=SQRT ( ( XT+TG*B/2 ) **2+ZZ + ( Y ( IZ ) -D/2 ) M2 ) 

R5 1 = SORT < ( XT+TGID/2 ) W2+ZZ+ < Y ( I Z ) +B/2 ) 1*2 ) 
Ro=5QRT ( ( XT+TG*S )**2+ZZ+< Y < IZ )-S) t*2) 

R6 1 =SQRT < ( XT+ TG*S ) **2+ZZ+ ( Y ( I Z ) +S ) **2 ) 

R?=ZZ+YY+XLS*XLS 

R8=ZZ+YY+XTS*XTS 

R1 T--SQRT ( XVT**2+ZZ+YY ) 

R2T=SORT < ( XVT +TG*SVT ) **2+ZZ+ ( Y ( IZ )-SMT ) t*2) 
R21T=SQRT<<XVT+TG*SVT)**2+ZZ+(Y(IZ)+SVT)**2) 
R3T=SQRT((XLT+TG*D/2)**2+ZZ+(Y(IZ)-D/2>**2) 
R31T=3QRT< (XLI+TG*D/2)**2+ZZ+(Y(IZ)+D/2)<*2) 
R4T=SQRT < ( XLT+TGIST ) **2+ZZ+< Y( IZ)-ST>**2) 

R4 1 T=SGRT ( ( XLT+TG*ST > **2+ZZ+ ( Y < IZ >+ST ) **2 > 
R5T=SQRT ( ( XTT+TG4D/2 ) **2+ZZ+<Y< IZ )-D/2) **2) 
R51 I=SQRT( ( XTT+TG*D/2 > **2+ZZ+( Y< IZ)+D/2)**2) 
R6T=$0RT ( < XTT+TG*ST ) W2+ZZ+ < Y ( I Z ) -ST )**2 ) 
RolT-SQRT ( ( XTT+TG*ST ) **2+Z2+ ( Y < IZHST ) **2> 


fGINAL 


IS 


° F p OOR QUALITY 


VNS<IZdX)=+MMDDIit<Z<IZ)+AlPHA*RN)/16/BETA*(l/<R7*SQRT<R7)) 
i -0.61/<R8*SQRT(R8))) 

THESE ABOVE ARE THE BODY SOURCE TERNS 

VNS( IZ*IX)=0<(Z(IZ)+ALPHA*RN)/PI/BETA<((TS*S+XL*T6-Y(IZ))/(Q3<R4)- 
(TS*D/2+XL*TG-Y(IZ))/<Q3*R3)+<TS*S+XL*TG+Y(IZ))/(Q4*R41)- 
1 < TS*D/2+XL*TG+Y ( IZ ) )/ ( Q4IR31 ) - ( TS*S+XT*TG-Y ( IZ) )/( Q5*R6 )♦ 

t ( T S*B/2+XT *TG- Y ( IZ ) ) / < 05IR5 ) - < TS*S+XT*TG+Y ( I Z ) ) / ( Q6*R61 ) + 

t (TS*D/2+XT*TG+Y(IZ))/(96*R51))+VNS(IZ>IX) 

THESE WERE THE WIHG LIRE SOURCE TERNS 

VNS<IZ>IX)=QT*<Z(IZ)+ALPHA*RH)/PI/BETA*((TS*ST+XLT*T6-Y(IZ))/ 

1 <Q3T*R4T>-<TS*D/2+XLT*TG-Y(IZ))/(Q3T*R3T)+<TS*ST+XLT*TG+Y(IZ))/ 

1 <Q4T*R41T)-(T3*D/2+XLT*TG+Y(IZ))/(Q4T*R31T)-<TS*ST*XTT*TG- 

1 Y < I Z ) ) / < Q5T*R6T ) + < TS*D/2+XTT*TG- Y < IZ) )/( QSTtRST ) - 

1 (TS*ST+XTT*TG+Y(IZ))/<06T*R31T)+(TS*D/2+XTT*TG+Y(IZ))/ 

1 (Q6T*R51T) )+VNS( IZ» IX) 

THESE WERE THE TAIL LINE SOURCE TERNS 

VNA ( IZ > I X ) =-D*D*ALPHA/8* ( ( Y Y-ZZ )/< YY+ZZ ) *121 ( 1-XLS/SQRT ( R7) ) 

1 +ZZ*XLS/W+ZZ)/(R7*SQRT(R7))-0,A1*(<YY-ZZ)/(YY+ZZ)**2* 

1 ( 1 -XTS/SQRT ( R8 )) +ZZ*XTS/ ( YY+ZZ )/ ( R8*SQRT ( R8 ) ) ) ) 

THESE ARE THE DOUBLET (BODY LIFT) TERNS 


TT=4/Q1*( < T S*SV+X V*TG-Y ( I Z ) ) /R2- ( XV*TG-Y (IZ) )/Rl ) 

T 0M/02N ( T S*SV+XV*TG+Y ( IZ) ) /R21- ( XV*TG+ Y ( IZ) )/Rl ) 

VNA<IZ»IX)=VNA(IZ»IX)-CIRC/4/PI*(F*TT+Fl*TO+(Y(IZ)-SV)/ 

1 (ZZ+( Y(IZ)-SV)**2)*(1-(XV+SV*TG)/R2)-(Y(IZ)+SV)/(ZZ+ 

1 (Y(IZ)+SV)**2)*(1-(XV+SV*TG)/R21)) 

THE ABOVE ARE THE WING LINE VORTEX TERNS 

TTT=4/QlT*((TS*SVT+XVT*TG-Y(IZ))/R2T-(XVT*TG-Y(IZ))/RlT) 
TOT=4/G2T*( (TS*SVT+XVT*TG+Y(IZ))/R21T-(XVT*TG+Y(IZ))/R1T) 


VNA(IZrIX)=VNA<IZ»IX)-CIRCT/4/PI*(FT*TTT+FlT*TOT+(Y(IZ)-SVT)/ 


GZHYUZ>-SVm*2)*(l-<XVT+SVTtTG)/R2T)-(Y<IZ)-fSVT)/ 

!ZZUYdZHSVT>**2)*(l-<XMnSVWG)/R21T)> 

THE ARE THE TAIL VORTEX TERMS 

VNdZ<IX>=4VNAdZ>IX)+VNS<IZ»IX) ORIGINAL PAGE IS 

QE EOOR QUALITY 

END DO 

DO 1Z-NY2+1 iNYO+NZ 'COMPLETE SIOEWAU 

fY- v (IZ)»*2 

L-.:UZ)+ALF‘HA*RM)«2 
F=(XV+Y< IZ>*TG) 

El-</V-fdZ)*TG> 

ft=ixvhygz;*tg> 

FT1*<XVT-Y< IZ)*TCi) 


Qi««(r*ftrs*zz) 

w2=4*<fltfltTS*ZZ> 

93=«<:XL + r(IZ)*TG)tt2+TS«Z) 

04=4*< «XL-YdZ»*TG)tt2+TS*ZZ) 

G5= j* c < XT+ Y ( I Z ) *TG ) «2+TS*ZZ ) 
s a * ' < XT- Y ( IZ ) >TG ) 112+TStZZ ) 

01T=4»(FT*FT+TS*ZZ) 

l£T»4*(FlT*FlWStZZ> 

Q3T=4* < < XLT+Y (IZ ) *TG) »*2+TS*ZZ ) 
G4T=4*<<XLT-Y(IZ>*TG)**2+TS*ZZ) 
C<5T=4*<<XTW<IZ)*TG)tt2+TS*ZZ> 
ILT=4K(XTT-YUZ)*TG)«2+TS*ZZ) 

Rl=SQRT(XVtt2+ZZ+Y(IZ)«2) 

R2-SQRT ( ( XV+TGISV ) W2+ZZ+ ( Y ( IZ ) -SV) t<2 ) 

R2 1 =SQRT < ( X V+TG*SV ) M2+ZZK Y ( IZ ) +SV )*»2 ) 

R3=SQRT( ( XL t TGtt/2) «2EZZ+( Y( IZ) -D/2 ) **2) 

R31-SQRK <XL+TG*D/2)«2+ZZHY< IZ)+D/2)«2) 

R4=SQRT ( < XL+TGf S)tt2+ZZ+ ( Y< IZ )-S) «2) 
R4t=S0RU(XL+TG*S)«2+ZZHY(IZ)+S)»2) 

R5*S8RT ( ( XT ETGID/2 ) **2+ZZ+ ( Y ( IZ ) -D/2) 112) 

R51=SQRT((XT+TG*D/2)«2+ZZ+(Y(IZ)+D/2)*l2) 

R6=SC)RT<<XT4TG*S)tt2+ZZW(IZ>-S)tt2> 

R6 1 -SORT ( ( XT+TGIS ) W2+ZZ+ ( Y ( IZ ) +S) «2) 

R7=ZZ+YYEXIS*XLS 

R8=ZZW4XTS*XTS 

RlT=SGRT(XVTtt2«ZW) 

R2T =SQR T ( ( X VT +TG*S VT ) **2+ZZ+ < Y ( I Z ) -SVT ) M2 ) 

R21 T -SORT ( ( XV T4TGISVT ) **2+22+ ( Y ( 12 ) +SVT ) 1*2 ) 

R3T=SQRT ( ( XLT+TGID/2 ) M2+ZZ+ < Y ( IZ ) -D/2 ) M2) 

R31T=SQRT< <XLT+T6*D/2>**2+2Z+(Y(I2>+D/2)**2) 

R4T=SQRT< ( XLT+TGtST ) *S2+ZZ+ < Y < IZ ) -ST ) M2 > 

S4 1 T=SORT< ( XLT+TG4ST )«2+ZZ+( Y( IZ)+ST)M2) 

R5T=S0RT ( <XTT+TG*D/2 ) **2+ZZ+< Y < IZ)-D/2 )**2) 
R51Tr$GRT<(XTT+TG*D/2>«2+ZZ+(Y<IZ)+D/2)tt2) 

R6T=SGRT ( ( XTTETGIST ) **2+ZZ+ ( Y ( IZ) -ST) **2) 
RolT=SGRT((XTT+TG*Sm*2+ZZ+<Y(IZ)+ST)«2> 

VNSUZ.IX >=-[i*P*DDD/ 1 6/BET A*Y ( I Z ) * ( 1/ ( R7*SQRT<R7) ) - 
V.6t/(R8*S0RTiR8))) 

THESE ARE THE BODY SOURCE TERMS ON THE WALL, 

VNS:IZ.IX)=VNS(IZfIX)-fl/PI/FETA*(((XL+Y(IZ)tTG)*(XL«*TG)+ZZ)/ 

Q3/R4-((XLEY(IZ)*TG)«XL+D/2tT6)+ZZ)/Q3/R3-((XL-Y(IZ)*TG)* 
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! aLfSITG)fZZ>/04/R41+<<XL-Y(IZ>*TG)*<XLM/2*TGHZZ>/ 

1 G4/R3l-( !XT+Y(IZ)tTG)t(XT+StTGHZZ)/Q5/R6+((XT4Y<IZ)tTG)l 

l < XT+0/2*TB ) +ZZ ) /05/R5+ ( < XT-Y< IZ)*TG)f(XT+S*TG)+ZZ)/Q6/R61 

1 - ( ( XT-Y <IZ ) *TG ) * ( XT+D/2XTGHZZ ) /Q4/R51 ) 

THE ABOVE ARE THE KING THICKNESS (LINE SOURCE) TERNS 

VNS< IZ> IX!=VNS< IZiIX)-QT/Pl/8ETAt< < <XLT+Y< IZ)tTG)*(XLT+ST*TG> 
TZZ ) /03T/R4T - ( ( XLT+Y (I Z ) * TG ) * < XLT+D/2*TG ) +Z2 ) /Q3T/R3T- 
((XLT-Y(IZ)*TG)*(XLTtSTXTG)+ZZ)/Q4T/R41T+<(XLT-Y<IZ)«TG)* 
<XLT+»/2iTG)+ZZ)/Q4T/R31T-((XTm(IZ)ITG)*(XmSTITG)+ZZ)/ 

Q5T /RfcT + < ( XTT+Y < IZ ) *T6 > $ ( XTTL D/2ITG ) +ZZ ) /05T/R5T + ( < XTT- 
r<IZ)*rG)*(XTT+STtTG)+ZZ)/06T/R61T-((XTT-Y(IZ)ITG)*(Xm 
D/2*TGHZZ)/06T/R51T) 

THESE ARE THE TAIL THICKNESS TERNS 

VNA(IZ»IX)=-PID>ALPHA/84Y(IZ)t<Z(IZHALPHAtRN)/(YY+ZZ)t<(l-XLS/ 
SORT (R7) >* 2 / ( YY+ZZ) -XtS/R7/SQRT(R7>-0.Mt< ( I-XTS/S0RKR8) )*2/ 

< YY f ZZ > -XTS/R8/SQRT ( R8 ) ) ) 

THE ABOVE ARE THE DOUBLET TERNS (BODY LIFT) 

TT-4/0 1*< ( TSISV+XVBTG-Y ( I Z ) ) /R2- ( XV»TG-Y( IZ ) ) /R1 ) 
t 0=4/021 ( ( TSISV+XVITS+Y ( I Z ) ) /R21 - ( XVITG4 Y< IZ))/R1> 

TTT=4/01T*((TS»SVTtXVT»TG-Y(IZ))/R2T-(XVmG-Y(IZ))/RlT) 
r0T=4/02TI< ( TS*SVT+XVTITG+Y( IZ ) ) /R2I T-<XVT*TG+Y( IZ ) )/RlT) 

VNAdZf IX)=VNA(IZ»IX)4CIRC/4/PI>(Z<IZ)4ALPHAtRN)t(TGt(TT”T0) 
t -l/(ZZ+(Y<IZ)-SV)«2)*U-(XV+SVITG)/R2)n/(ZZ+(Y(IZ)4SV)*t2) 

I *(i-(XVfSV*TG)/R21)> 

C THESE ARE THE TERNS FROM THE SWEPT LINE VORTEX 

VNA<IZ*IX)=+CIRCT/4/PI*(Z(IZ)+ALPHA*RN)l(TGt(TTT-T0T)-l/ 

1 <ZZT(YnZ)-SVT)*l2)*(l-(XVT+SVT*TG)/R2T)+l/(ZZ4(Y<IZ)«MT) 

1 «2)*(1-(XVT+SVT*TG)/R21T))4VNA(IZ»IX) 

C THESE ARE FOR THE TAIL VORTEX SYSTEN 

VN(IZ»IX) a VNA(IZ»IX)+VNS(IZ>IX) 

END 00 'COMPLETE TUNNEL IS NON DONE 

END DO 


DO IX=-20>0 


1 

1 

1 

t 

4 

4 

C 


t 

l 

c 


C STRETCHED DISTANCES FRON JACK POINT TO INTERSECTION OF C/L AND* 

XV=(XRVC-XJ(IX))/BETA ! WING VORTEX 

XL=(XRVC-C/4-XJ(IX))/BETA ! WING LEADING EDGE 

XT=(XRVC+3tC/4-XJUX))/B£TA ! WING TRAILING EDGE 

XLS=(XSL-XJ(IX))/BETA ! BODY SOURCE 

XTS=(XST-XJ(IX) )/BETA ! BODY SINMAFT) 

XVT=(XRVTC-XJ(IX))/BETA ! TAIL VORTEX 

XLT=XVT-CT/4/BETA ' TAIL LEADIN6 EDGE 

X TT=X VTF34CT /4/BET A » TAIL TRAILING ED6E 


1 M'-iiivAL PAGE 13 
OF POOR QUALITY 


C RN IS THE ARH LENGTH FRON THE POINT OF ROTATION OF 

C THE BALANCE/STING SYSTEN TO TIE CP ORIGIN OF COORDS. 

C THE NODEL IS DISPLACED UPWARDS BY A DISTANCE ALPHAIRN 


RH=0.501269 ! METERS 
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DO IZ=1»NY2 
F= ( XV+ Y < I Z > * TS > 

F1=UV-YUZ>*TG) :v;V.JLr V AO!5 IS 

Ff=( <UT+Y( IZJITG) OF POOR QUALITY 

FTl = OfVT-Y< IZXTG> 

YY-Y<IZ>M2 

ZZ=<ZOZ)+ALPHA*RM)M2 

Gl=4*<t:«F+TS*ZZ> 

a z -4* < F 1 *F 1 + T 5»Z2 ) 

< i XlHMZ)*TG)M2+TS*ZZ) 

■'M-o*- X|_-r<IZ)»TG)M2tTS*ZZ) 
i;S=4« <r^"£i»T6)*l2+TStZZ) 

QA=4*( (XT-v<IZXTG)**2+TS*ZZ) 

Rir=4t(FT*F7+T5*ZZ) 

02I-4*<F1T*F1T+TS*ZZ) 

UTM*i(XLT+Y(IZ'iTG)»*2+TStZZ) 

34T=«( (XLT-r< IZ)*TG)M2+TS*ZZ) 

G5T=4* ( ( X rm ( I Z ) *TG ) M2+TSYZZ ) 

C(6T=4*((XTT-Y(IZ)*TG)**2tTS*ZZ) 

Rl=SGRT(xVM2+ZZW) 

R2=SQRT ( 4 XV+TG*SV)«2+ZZ+ < Y < IZ > -SV)M2 ) 
R21=SGRT<(XV+TG*SV)M2+ZZ+<Y<IZ)+SV)M2) 

R3 s SGRT ( < XL+TG*D/2 ) M2+ZZ+ ( Y ( IZ ) -D/2 ) M2 ) 
R31=SGRT<<XL+TG»D/2>M2+ZZ+(Y(IZ)+D/2>M2> 

R4-SQRT ( ( XL+TGCS ) **2+ZZ+ < Y ( IZ ) -S) **2) 
R4l=SGRT<<XL+TG*S)M2+ZZ+(Y<IZ)+S)M2) 

R5=SQRT ( ( XT+TGID/2 >**2+ZZ+ <Y(IZ> -D/2)M2) 

R51 =S0RT < < XT+TGtD/2 > M2+ZZ+ < Y < IZ > +D/2) M2 ) 

R6=SGRT ( ( XT+TGIS ) M2+ZZ+ ( Y ( IZ )■ -S ) M2 ) 
R61-SQRT(<XT+TG*StM2+ZZ+<Y(IZ)+S>M2> 

R7=ZZ+YY+XLS*XLS 
R8=ZZ+ YY+XTS4XTS 

R1 T=S0RT < XMTM2+ZZFYY ) 

R2T*SGRT ( < X V WGISVT ) M2+ZZ+ < Y ( IZ ) -SVT ) M2) 

R2 1 T = SORT ( < XO WG*Sy T ) M2+ZZH Y < IZ HSVT ) M2) 
R3T=SQRT((XLTFTG*D/2)M2+ZZ+(Y(IZ)-D/2)M2) 

R31 T=SQRT( ( XLT+TGID/2 ) M2+ZZ+ ( Y( IZHD/2) M2) 
MT=SGRT<<XLT+TGtST)M2«ZHY(IZ)-ST)M2) 

R4 1 T=SQRT< < XLT+TGIST ) M2+ZZ+ ( Y ( IZ ) +ST ) M2 ) 
R5T=SQRT(<xmTGtB/2)M2+ZZ+(Y(IZ)-D/2)M2) 

R5 1 T=SQRT ( < XTT +TGID/2 ) M2+ZZ+ ( Y ( IZ ) +D/2 ) M2 ) 

R6T*SQRT ( ( XTT +TGYST ) M2+ZZ+ ( Y ( IZ ) -ST ) M2 ) 

R6 1 T=SQRT ( ( XTT+TGIST ) M2+ZZ+ ( Y ( IZ ) 1ST ) M2 ) 


VNS ( I Z • I X ) =-D<D*DDDM Z ( I Z ) + ALPHAIRM ) / 1 6/BET A* ( 1 / ( R7ISGRT ( R7 ) ) 

1 -0,61/<R8tSGRT(R8))) 

C THESE ABOVE ARE THE BODY SOURCE TERMS 

VMS < I Z - IX » = -G*< Z ( I Z ) TALPHAYRM ) /PI/BETAt ( < TS*S+XL*TG-Y ( IZ ) ) / ( Q3IR4 ) - 
1 <, TStO/2+XLtTG-Y ( IZ ) )/(G3tR3)+(TS<S+XL<TG+Y( IZ) )/(Q4tR41)- 

l <TS*D/2fXL*TG+Y(IZ>)/(Q4*R31)-(TS>S+XT*TG-Y(IZ))/<Q5*R6)+ 

t (TS*D/2+XTITG-Y(IZ))/<Q5tt5)-(TS*S+XT*TG+Y(IZ))/(06IR«)+ 

1 <TS«i/2+XT»TG+Y( IZ) )/< Q6IR51 )HVNS( IZ»IX) 

C THESE WERE THE WING LINE SOURCE TERMS 

MNS(IZiIX)=-QT*(Z(IZHALPHAtRH)/PI/8ETA>( (TStST+XLTtTG-YdZ))/ 
<G3T*R4T>-<TS*0/2+XLT*TG-Y<IZ))/<Q3T*R3T)WS*ST+XLTITG4Y(IZ))/ 

( Q4TYR4 1 T ) ■ - ( TSHD/2+XL T4TG+Y ( I Z ) ) / ( G4TIR31T ) - ( TSIST *XTT*TG- 


I 

1 
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1 Y(IZ) > / < Q5TIR6T ) + » T S4D/2+XTTITG- Y ( IZ ) )/<05TIR5T ) - 

1 nS*STWT*TG+Y(IZ))/<Q6T*Mim<TS*D/2«TTtTG*r<IZ)>/ 

1 <06T*R51T))+MNS(IZ.IX) 

THESE HERE THE TAIL LINE SOURCE TERNS 

VNAi I Z » I X ) =DIDIALPHA/8t < ( Y Y-ZZ )/ ( YY+ZZ ) «2> ( 1 -XLS/SORT (R7 ) ) 
1 +ZZ»XL5/(YY+ZZ)/(R7*SQRT(R7))-0.61*((VY-ZZ)/(YY+ZZ)«2t 

1 < 1-XTS/S0RT(R8>)+ZZ$XTS/(YY+ZZ)/(R8*S(IRT(R8)))) 

THESE ARE THE DOUBLET (BODY LIFT) TERNS 


f T-4/01K ! TStSVfXV»TG-Y(IZ))/R2-(XVtT6-Y(IZ))/Rl) 

T0=4/Q2*( ( rS*SV+XV*TGFY< IZ) )/R2l-<XVtTG4Y< IZ) )/Rl ) 

VNAi LM X <=VNA( 12, UHCIRC/4/PI*<F*mFl*T0+<Y(IZ)-SV)/ 

1 '.ZZ-K Y<IZ)-SV>W2>*( l-(XV+SVtTG)/R2)-(Y(IZ)FSV)/(ZZt 

1 ( Y ( IZ } 4SV > t*2 ) * < 1 - < XV4SVtT6 ) /R21 ) ) 

THF ABOVE ARE THE UING LINE VORTEX TERNS 

TTTM/Ol T*( ( TStSMTFXVTtTG-Y(IZ) ) /R2T-(XVT*T6-Y< IZ) )/RlT) 

T0TM/I12TK <TS*SVTtXVT*TGf Y( IZ) )/R21T-(XVT»TGFY(IZ) )/RIT) 

VNA<IZ*lX>=VNA(lZ*IX)+CIRCT/4/PI$(FTtTTT4FlTIT0T+(Y<IZ)-SVT)/ 

1 (ZZ+(Y(IZ)-SVT)W2)*a-(XVT«VTITG)/R2T)-(Y(IZ)FSVT)/ 

1 (ZZt<mZ)+SVT)«2)*<l-<XVT+SVT*TG)/R21T)) 

THE ABOVE ARE THE TAIL VORTEX TERNS 

VNJ(IZ*IX) s VNA( IZi IX)+VNS< IZ* IX) 

END DO 

00 IZ=NY2+NZ41*NY4NZ 'UPPER UALL 
F=(XV+Y(IZ)*TG) 

F1=<XV-Y(IZ)*TG> 

FT='XVT+Y(IZ)ITG) 

FTl-(XVT-YaZ)*TG) 

YY=YdZ)«2 

ZZ=(Z(IZ)+ALPHA*RN)t*2 

01Mf(F*FtTS*ZZ) s ‘ ZZ IS 

U2=4*<F1*F1KS*ZZ) Uf' i'OQR OU ATJT^ j 

03=4*< <XLFY(IZ)tTG)»*2+TS*ZZ) 

U4=4*(\XL-Y<IZ)fTG)*t2+TS*ZZ) 

05M*< (XT+Y(IZ)*TG)«2+TS*ZZ) 

Qd=4*< <XT-Y(IZ)»TG)**2*TS*ZZ) 

01T=4*<FT*FT4TS*ZZ) 

Q2TM*(F1T*F1T4TS*ZZ) 

Q3T=4*( IXLT+YI IZ)*TG)»*2+TS*ZZ) 

04T=4*( (XLT-YI IZ)*TG)«2+TSIZZ) 

Q5T=4i ( ( XTT4Y ( IZ ) *TG ) 112+TStZZ ) 

0oT=4>( ( XTT-Y ( IZ ) ITG ) W24TSIZZ) 

Rl=3QRT(XV*$2+ZZ+YY) 

R2=SQRT ( ( XV+TGBSV ) < t 2+ZZ+ ( Y ( IZ ) -SV)$*2) 

R21=SQRTUXVFTG*SV)**2+ZZW<IZ)4SV>«2) 

R3=SBRT(<XLFTG*D/2)«2«Z4(Y(IZ)-D/2)«2) 

R31=S0RT < ( XL4TGtD/2)»2+ZZ+ ( Y ( IZ)+D/2)«2) 

R4=S0RT((XL4TGtS)W2FZZHY(IZ)-S)«2) 

R4l=S0RT(<XL+TGtS)«24ZZHYIIZ)+S)t*2) 

R5=S0RT((XT+TG*D/2)«24ZZ+(Y(IZ)-D/2)«2) 

R51=S0RT((XTtTG*D/2)W2«ZF(Y(IZ)4D/2)«2) 

R6=SGRT(<XWG*S)tt2+ZZHY<IZ)-S)«2) 
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RA1=SQRT ( ( XT fl GIS ) **2+ZZ+<Y <IZ)+S)**2) 

R7=ZZW+XL3tXLS 

R8=ZZ+YY+XTS*XTS 


i’AGS IS 

,0E POOR QUALITY 


R1T=SQRT(XVTM2IZZW> 

R2T=SQRT ( ( XV WGBSVT ) M2FZZ+ ( Y ( IZ)-SVT ) M2 ) 
f'2 1 r=3GRT < ( XVT+T GtSVT > M2+ZZ+ ( Y< I Z ) +SVT ) M2 ) 
R3T-SQRT( (XLTFTG*D/2)M2+ZZHY(IZ)-D/2)M2) 
K31 T-5QRT t ( XLT+TG*D/2 ) M2+ZZ+ ( Y< IZ J+D/2 )**2) 
R4T=SGRT(<XLT+TG*ST>M2+ZZW(IZ)-ST)M2) 
RUT-SQRT ( XLT+TG*ST)M2+ZZ+<Y(IZ)+ST)M2> 
R5T=SQRT(( XTT+TGM/2 ) M2+ZZ+ < Y ( I Z ) -D/2 ) M2 ) 
R5 1 T -SORT < ( XTT f TG4D/ 2 ) M2IZZH Y ( IZ ) +D/2 ) M2 > 
Ri T =SQRT ( ( X TT+TGBST ) M2+ZZ-K Y ( IZ ) -ST ) M2) 

Rot T =SQRT ( ( XT WG*ST > M2FZZ+ ( Y ( IZ HST ) M2 ) 


WS:lZ.IX)=+[itD*DDD*(Z(IZ)+ALPHAIRM)/16/BETA*a/(R7tSQRT<R7)) 

1 -0 » 61 / < R 8ISGRT < R8 ) ) ) 

THESE ABOVE ARE THE BODY SOURCE TERMS 

t'NS(:Z,IYi=G*!ZUZ)+ALPHA»RM)/PI/BETAK(TS*S+XL*TG-Y(IZ))/(Q3*R4)- 
1 <TS»0/2+XL*TG-Y(IZ))/(Q3*R3)HTS*S+XL*TG+Y(IZ))/<Q4*R41>- 

l aS*D/2fXL*TG+Y(IZ))/(04*R31)-(TS<S+XTITG-Y(IZ))/(Q5IR6)+ 

1 (TS*D/2+XTITG-Y(IZ))/(G5*R5)-(TSIS+XTITG+Y(IZ))/(Q6tR61)+ 

1 ( TS4D/2+XTITG+Y ( IZ) )/< Q6IR51 ) )+VNS< IZ > IX) 

THESE WERE THE WING LINE SOURCE TERMS 

WS(IZUX)=QT*(Z(IZ)+ALPHA*RM)/PI/BETAt((TS»ST+XLT*TG-Y(IZ))/ 

1 <Q3T*R4T) - 1 TS*D/2+XLT*TG-Y( IZ ) >/(Q3T*R3T)+< TS*ST+XLT*TG+Y< IZ) )/ 

1 ( 04TIR4 1 T ) - ( TS4D/2+XLT* TG+Y ( IZ ) ) / ( Q4TIR31 T)-( TSBSTtXT T»TG- 

1 Y < I Z ; )/ ( Q5T4R6T ) + < TS*D/2+XTT*TG-Y < IZ ) ) / < Q5TBR5T )- 

1 ( TSISTIXTT ITG+Y (IZ ) » / ( Q6T4R61 T) + ( TS*D/2+XTT*TG+Y( IZ) )/ 

1 <06T*R51T)HVNS(IZ f IX> 

THESE WERE THE TAIL LINE SOURCE TERMS 

VNA( i:..IX)=-D*D*ALPHA/8*( < YY -ZZ ) / ( Y YIZZ ) M2* ( 1 -XLS/SQRT ( R7 ) ) 

I +ZZ*XLS/(YY+ZZ)/(R7*SORT<R7) )-0.61*( (YY-ZZ)/(YY+ZZ)M2» 

1 ( 1-XTS/SGRT (R8) )+ZZ*XTS/< YY+ZZ)/(R8*SQRT(R8) ) ) > 

THESE ARE THE DOUBLEKBOBY LIFT) TERMS 


TT=4/Q 1 * ( ( TSISV+XVtTG- Y <IZ))/R2-( XVtTG-Y ( IZ) ) /Rl ) 
T0=4/02I((TS*SV+XVITG+Y(IZ))/R21-(XV*TG+Y(IZ))/R1) 

VNA < I Z » I X )=gNA ( I Z f I X ) -C I RC/4/PI * ( F*TT+Fl*TO+ ( Y (IZ ) -SV) / 

1 (ZZ+(Y<IZ>-SV)M2)*<1-<XV+SV*TG)/R2)-<Y<IZ)+SV)/(ZZ+ 

1 <Y(IZ)+SV)M2)*<1-(XVISV*TG)/R21)) 

THE ABOVE ARE THE WING LINE VORTEX TERMS 

TTT=4/Q1T*((TS*SVT+XVT*TG-Y(IZ))/R2T-(XVT>TG-Y(IZ))/R1T) 

T0T=4/02T*((TS*SVT+XVTXTG+Y(IZ))/R21T-(XVT*TG+Y(IZ))/R1T) 

VNA(IZ» IX)=VNA(IZ»IX)-CIRCT/4/PI<(FT4TTT+FlTIT0T+(Y(IZ)-SVT)/ 
1 (ZZ+<Y(IZ)-SVT)M2)4tl-(XVT+SVTtTG)/R2T)-(Y(IZ)+SVT)/ 

1 (ZZKY(IZ)+SVT)M2)»(1-(XVT+SVT*TG)/R21T)) 

THE ABOVE ARE THE TAIL VORTEX TERMS 

VNJ<IZ»IX)=+VNA'IZfIXHVNS<IZ»IX> 


END DO 
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DO I2=^Y2+1 -■> MY2+NZ ! COMPLETE SIDEWALL 

YY=Y(IZ)t*2 

ZZ= ( Z ( I Z ) 4 ALPHAIRH >1*2 
F=(XV4Y<IZ>*TG> 

Fl=(XV-Y(IZ)*TG) 

FT=(XVT4Y(IZ)*TG) 

FTl»iXVT-Y<lZ)lTG) 

01 -4»<F*F+TS*ZZ) 

()2=4*(FltF14TStZZ) 

1)3=1*' :XL+Y<IZ)*TG)**2+TStZZ) 
l!4=4*((XL-Y<IZ)fTG)«24TS*ZZ) 

15=4*< (XT+Y( IZ)tTG)lt24TSIZZ) 

«**»*< (XT-Y(IZ)tT6>tt24TS*ZZ> 
M1T=4*<FTIFT+TS*ZZ> 

.‘.?T=4*(F1T*F1T4TS»ZZ) 

Q3T=4*(<XLT+Y(IZ)*TG)«2+TS*ZZ) 

U4T=4*<(XLT-Y(IZ)*TG)»*24TS*ZZ) 

QST-4* ( ( XTT4Y ( IZ) ITG) II24TSIZZ) 
l)6r=4*((XTT-Y(IZ)tTG)W2+TS*ZZ) 
rt 1 =SGR T ( XVII24ZZ4 Y ( I Z ) 1*2 ) 

R2=S9RT l < XV4TGISV ) W24ZZ4 ( Y< IZ )-SV ) M2 ) 
R2l=SQRT<<XV4TG*SV)**24ZZ4(Y<IZ)4SV)«2) 
R3=SQRT ( ( XL4 TGID/2 ) **24ZZ4 ( Y ( IZ ) -D/2) **2 ) 
R3l=SQRT < < XL4TGID/2 ) W24ZZ4 ( Y ( IZ ) +D/2) «2) 
R4=SORT( ( XL4TG*S)**24ZZ4<Y< IZ)-S)**2) 

R4 1 =$QRT < ( XL+T6IS ) M24ZZ4 ( Y ( I Z ) 4S ) If 2 ) 

R5=S0RT ( ( XT4TG*D/2) M24ZZ4 ( Y ( IZ)-D/2)«2) 

R51=SQRT( < XT 4TGID/2 ) t*24ZZ4 ( Y < IZ ) 4D/2) *12 ) 

R6=S8RT ( < XT4TG*S ) **24ZZ4 ( YIIZ )-S ) 112) 

R6l=SQRT((XT4TG*S)**24ZZ4<Y(IZ)4S)**2) 

R?=ZZ4YY4XLS*XLS 

R8=Z2+rY+XTS*XTS 


ORIGINAL PAGE IS 
OF POOR QUALITY 


R1 T =SQRT ( XVT**24ZZ4YY ) 

R2 T=SGRT ( < X VT 4 TG*SVT ) **24ZZ4 ( Y ( I Z ) -SVT ) **2 ) 

K2lT=SGRT(<XVT4TG*SVT)**24ZZ4<Y<IZ)4SVT)**2) 

F'3T=SQRT((XLT4TG*D/2)**24ZZ4<Y(IZ)-D/2)**2) 

R3 ! T=SQRT< < XLT4TG*D/2) **24ZZ4 ( Y< IZ ) 4D/2 ) *<2 ) 

R4 T =SQRT ( < XLT +TGSST ) <*24ZZ4 < Y ( IZ ) -ST ) **2 ) QUALITY 

R41T=SQRT<(XIT4TG*ST)**24ZZ4(Y(IZ)4ST)**2) 

R5T=SQRT ( < XTT4TG*D/2 ) W24ZZ4 < Y ( IZ ) -D/2 > **2 ) 

R51 T =SGRT < ( XTT4TG*D/2 ) **24ZZ4 ( Y ( IZ ) 4D/2 ) **2 ) 

R6T=SQRT < ( XTT4TG*ST )**24ZZ4 < Y ( IZ )-ST ) **2 ) 

R61T=SQRT((XTT4TG*ST)**24ZZ4<Y(IZ)4ST)**2) 


VNS ( I Z * IX ) =-D*D*DDD/ 1 A/BET A*Y< IZ ) *< 1/ < R7ISQRT ( R7 ) ) - 
I <>.61/<R8*SGRT(R8))) 

C THESE ARE THE BODY SOURCE TERMS ON THE WALL. 

VNS<IZ»IX)=VNS<IZ»IX>-Q/PI/BETA*(<(XL4Y(IZ)*T6)*<XL4S*T6>+ZZ)/ 
1 Q3/R4-) (XL4Y(IZ)*TG)*(XL40/2*TG)4ZZ)/Q3/R3-((XL-Y(IZ)*TG)t 

1 ( XL+SITG ) 4ZZ ) /Q4/R41 4 < ( XL-Y ( I Z ) ITG ) $ ( XL4D/2*T6)4ZZ)/ 

1 Q4/R31-((XT4Y(IZ)*TG)*(XT4S*TG)4ZZ)/Q5/RA4<<XT4Y<IZ)*TG)* 

1 ( XT4D/2*TG ) +ZZ ) /05/R54 ( ( XT-Y ( IZ ) ITG ) * ( XT4S*TG)4ZZ ) /QA/RA1 

1 -<<XT-Y(IZ)*TG)*<XT4D/2*TG)4ZZ)/QA/R51) 

C THE ABOVE ARE THE WING THICKNESS (LINE SOURCE) TERNS 


VHS(IZfIX)=VNS(IZfIX)-QT/F'l'BETAK((XLT+Y(IZ)*TG)l(XLT+STtTG) 



+ZZVQ3T/R4T-' (XLT+l'f IZ)*TG)t(XLT+D/2tTGHZZ)/Q3T/R3T- 

< <XLT-Y< IZ )*TGi*(XLTfST*TG)+ZZ)/04T/R41T+< (XLT-Y( IZ)*TG)* 

(XLT+D/2*TG)IZZ)/Q4T/R31T-((XTT+Y(IZ)*TG)*(XmSTITG)+ZZ)/ 

05T/R6T+((XTTtY(IZ)*TG)*(XmD/2*TG)«Z)/Q5T/R5T+<(XTT- 

Y(IZ)ITG>*(XTT+ST*TG)+ZZ)/0AT/R61T-UXfT-Y(IZ)*TG)*(Xm 

D/2*TG)+ZZ)/Q6T/R51T) 

THESE ARE THE TAIL THICKNESS TERHS 

VNA(!Z,IX)=-D*D*ALPHA/8*Y(IZ)*<Z<IZ)+ALPHA*RN)/(YY+ZZ)*(<1-XLS/ 

sGRT(R7))*2/(YY+ZZ)-XLS/R7/SQRT<R7)-0.6U((l-XTS/SGRT(R8))*2/ 

< YY +ZZ - -XTS/R8/SQRT <R8 ) ) ) 

THE ABOVE ARE THE DOUBLET TERHS (BODY LIFT) 

TT=4/01 *( < TS*SV+XU*TG"Y( IZ) )/R2-(XV*TG-Y( IZ) )/Rl ) 

TO=4/Q2*( (TS<SV+XVtTG+Y(IZ))/R21-(XVITG+Y(IZ))/Rl) 

TTT=4/QlT*<(TS*SVT+XVr*TG-Y(IZ))/R2T-<XVT*TG-Y(IZ) )/RlT) 
ri)T=4/i?2TI< (TStSVT+XVTtTG+Y(IZ))/R21T-(XVTITG+Y(IZ))/RlT) 

v'NA(IZ.IX)=WA(IZ>IX)+CIRC/4/PI*(ZUZ)+ALPHA*RH)XTG*(TT-T0) 

-i/(ZZ+(Y(IZ)-Sy>**2)Kl-(XV+SV*TG)/R2)+l/(ZZ+(Y(IZ)+SV)«2) 

*<1-<XV+SV*TG)/R21)) 

THESE ARE THE TERHS FROM THE SWEPT LINE VORTEX 

VNA ( I Z . IX ) =+C I RCT/4/P r * < Z < I Z ) +ALPHA*RM ) * ( TG« < TTT-TOT) -1/ 
(ZZ+(Y(IZ)-SVT)W2)*(l-<XmSVT*TG)/R2T)+l/<ZZ+(Y(IZ)+SVT) 
**2)*(l-<XVT+SVWG)/R21TmVNA<IZ,IX) 

THESE ARE FOR THE TAIL VORTEX SYSTEM 

VNJ( IZ « IX >=VNA( IZ r IX ) tVNS< IZr IX) 

END DO (COMPLETE TUNNEL IS NOW DONE 

END DO 

WE NOW WRITE A FILE FOR VNO TO FEED JACK.DISPL! 

FIRST WE MUST DEFINE A NEW X-WISE COORDINATE MEASURED 
FROM THE HALF INCH POINT ON THE WALL PLATES. 


DO IX=-NX*NX 

X-L0C(IX)=X(IX)+(27. 510.0254) !CP IS LOCATED AT 28 INCH POINT 
END DO 

0PEN<UNIT=8*NAME='VN0,DAT'*STATUS='NEW') 

WRITE (8. 40) CL* MACH 

FORMAT ( 'THIS DATA IS FOR CL='F4.2*3X* 'HACH='*F4,2) 

WRITE(8*t) 1.79959, 0.3302 
WRITE ( 8»42) 1.0 

FORMAT ( 'FREESTREAH VELOCITY = '*F6.1) 

WRITE(8<43) 

FORMAT (UNIT OF MEASUREMENT IS METERS') 

WRITE (8.44) ( 2*NX+1 ) 

FORMAT (NUMBER OF X LOCATIONS = '*15) 

WRITE < 8f 45 ) 

FORMAT* 'LOCATIONS') 

WRITE(8*D (X.LOC(IX).IX=-NX,NX) 

WRITE(8*46) NY 

FORMAT ( NUMBER OF PANELS IN Y DIRECTION = ',13) 

WRITE<8*47) (A/NY) 

FORMAT ( 'WIDTH OF PANELS ^ '*F6.5) 

URITE(8*48) 



FORMAT < 'LOCATIONS} ' ) 

NRITET8i4) (Y(IZ)i IZ=1»NY2) 

URITE(8»49) HZ 

FORMAT ('NUMBER OF PANELS IN Z DIRECTION = M3) 
URITE(8»50) TA/NZ) 

FORMAT < 'WIDTH OF PANELS = OF6.5) 

WRITET8.51) 

FORMAT (2X >//////) 

HO IX=-NX*NX 

WRITEi3?4) (VN(IZ»IXMZ=1>NY+NZ) 

END DO 
CLOSETS) 

WE NON WRITE A FILE FOR VNJ TO FEED JACOISPL! 

REDEFINE THE JACK LOCATIONS NON FROM THE i/2 INCH STATION 

XJ(-20)=4, 251.0254 
XJ(-19)*10. 4,0254 
XJT-18M5.4.0254 
XJT-17M9, 4.0254 
XJ(-16)=22,4.0254 
XJ(-15)=24.4.0254 
XJ(-14)=25. 54.0254 
X J < - 1 3 ) =27 » 4 , 0254 
XJ(-12) s 28.5t,0254 
XJ(-11)=30. 4.0254 
XJ(-10)=31.54,0254 
XJ< -? % =33,4,0254 
X.K -8)=35. 4.0254 
XJT-7)=37.4,0254 
XJ(-6A39.4,0254 
< J(-5)-42 .4.0254 
XJ{-4)=46.4,0254 
XJ(-3)=51.4.0254 
XJ(-2)»54. 4.0254 
XJ<-1)=61. 4,0254 
>:J<0)=66, 54,0254 


ORIGINAL PAGE IS 
OF POOR QUALITY 


0PENTUNIT=7»NAME='VNJ.BAT'fSTATUS*'NEH') 

WRITET7r61) 17 

FORMAT ('NUMBER OF JACKS -M3) 

WRITE(7»62) 

FORMATT 'DISTANCES IN METERS') 

URIT£(7rf) (XJ(IX)»IX=-20>-4) 

WRITE(7,63) 

FORMAT! 'NUMBER OF Y PANELS '5X. 'NUMBER OF Z PANELS') 

WRITE<7»4> NY» NZ 

WRITE<7»64) 

FORMAT 1 2X » //// / ) 


DO IX=-20»-4 

WRITE(7»4) (VNJ<IZ»IX)>IZ S 1»NY+NZ) 
END BO 
CLOSET 7) 
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APPENDIX E - VEEXPHINO FORTRAN LISTING 


Primary Symbols 


A 

AO, AX 


D 

DYDZ 

GO 


MACH 

MEASVN , MEASVX 

PHINO 

VX 


VN 


Tunnel height 

Normal velocities at panel centers due 
to horseshoe vortices on a panel in a 
ring at IX * -NX, fore and aft symmetric 
and unsymmetric contributions 
Tunnel breadth 

Ratio of floor panel width to wall panel 
height 

Axial disturbance velocities at panel 
center due to a Green's source panel in 
a ring at IX « -NX 
Mach number 

Measured wall values of normal and axial 
disturbance velocities 
Computed wall-free normal velocities at 
panel centers 

Calculated axial velocity increment due 
to measured normal velocities (see 
report) 

Residual normal velocity field to be 
nulled by the walls. Input into 
PHIXZM provides flow distortion field at 
model 
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PROGRAM VEEX-PHINO 


C 

C 

C 

C 

C 

C 

c 

c 


THIS PROGRAM COMPUTES A FIRST GRDER APPROXIMATION TO 
WALL-FREE NORMAL VELOCITIES AT RECTANGULAR TUNNEL WALLS 
PRODUCED BY A TEST MODEL. INPUTS REQUIRED ARE MEASURED 
NORMAL AND AXIAL DISTURBANCE VELOCITIES AT the WALLS OR 
RECTANGULAR CONTROL SURFACES DURING A TEST. 

THE FIRST PART OF THE COMPUTATION GETS THE AXIAL DISTURBANCE 
VELOCITY FIELD AT THE WALLS CAUSED BY THE MEASURED TEST 
NORMAL VELOCITIES AT THE WALLS. 


A RECTANGULAR TUNNEL IS ASSUMED OF HEIGHTS AND 
WIDTHS 


SELECT AN EVEN!! NUMBER OF PANELS ON THE VERTICAL 
UALL(NZ). THIS FIXES THE HEIGHT OF T HE PANELS WHICH IS 
ALSO TAKEN AS -THE LENGTH OF THE PANELS. THE LENGTH OF 
TUNNEL CONSIDERED CAN BE SET BY SELECTING THE NUMBER OF 
PANELS FOREWARD AND AFT OF THE ORIGIN. NX. 


PRANDTL-GLAUERT COMPRESS I BILTY CORRECTIONS ARE USED 


C 

C 

c 


A PARAMETER 'M' IS USED TO FIX THE NUMBER OF SOURCE 
IMAGES USED BEFORE SMEARING TO OBTAIN A CLOSED FORMULA 
FOR THE REMAINING SOURCES OUT TO INFINITY. 


DIMENSION X(-80:80). r<80).Z<80). 

AO<10»80)»YBN<5»80)>TOP<5.80MSM<5.80MSP<5.80)»YD<5>. 

YP<5).YD6<16580).YS6U6:80)»ZD<6115)»ZDP2(6510.6:i3). 

ZBN2<65l0.6515).ZBF'<6510MZBM<6510MZl<5»o5lS).rZ2?5.4515>. 

YZ3<5»6515)»YZ4<5»65i5).ZYl(6510.5).ZY2(4:iC.5>.?T3i6:iO»5>. 

ZY4(6!10»5)>ZBP3(6510.6515).ZDH3(6!10«6515>.TD1!5»16!P0). 

YD2(5» 16580). YB3<5» 16580 .YD4( 5. 16530) »TRl <6510* 16560) • 

YR2(6! 10.16580) »YR3<65 10. 16580 >«YR4(651M6!90>f 
AX(10.80»-80580).PHIN0<80. -80580). VX(8C.-90580/. 
MEASVN(80.-80580).M£ASVX(80.-80580).VN(80. -80580). W(80.-H0580). 
PT < 20 5 20 » 40 . 80 ) . P J ( -20 5 20 • 40 . 80 > • GO ( 40 . 80 . -80 5 80 » • 
0(0520.40.80). R(0520.40.80)»S(0520.40.80M(0520.40.60). 
YPP(40.80)»YDD(40»80) 


10 


PARAMETER <PI=3.1415?) 

REAL MACH 
TYPE 10 

FORMATdOX. 'ENTER NX NZ MACH A D AND M'> 
READ(5»t) NX. NZ. MACH. A. D.M 
BET A=SORT ( 1 -MACHIMACH ) 

NZ2=NZ/'2 

NY2=NINT(D*NZ2/A+0,1) 

NY=NY2»2 

now set up the coordinates of the panels. 

DO IX=-NX?NX 
X(IX)=IXtA/NZ'BE T A 
END DO 
DO IZ-1 *NY2 
Y(IZ)=(IZ-0.5)ID/NY 
Z(IZ)=A/2 


< J> i’l . % .'ll l\L 1 J i *A.J ]-V A KJ 

OF; P(X)R QUALITX 



I>0 IZ=(NY2+1) • <WY2+VZ> : -* 

Y<IZ)=D/2 OF POOR 

Z(IZ)*((MZWm)/2.00-!Z>«-W 

END DO 

DO IZ=(NZ+NY2tl inNZ+N'f ■ 

Y < IZ ) =< ( NZ+N* +0 . 5 > - ! • ■ #?«/*•* 

Z(IZt--A/2 
END DO 


DYDZ=D/A*NZ/NY *WT:.: ;,F f ANEL FLOOR Af-D WALL WIDTHS 

SET UP COMMON REPEALING TEPfiS 'V SAVE TI.iE 1 

DO IZ=1.NY+NZ 
DO JZ r l»NZ2PNY2 
YDD<JZ.IZ»=Y(JZ>-r>lZ 
YPPUZ.IZ»=YfJ2)+Y;I!» 

DO 1=-M»M 

pt< i .JZriz>=:2*A*i+zijz»-z rz i>«: 

P J ( I • JZ . I Z » = ; 2*A»I+A-Z< JZ >-? ' : Z :■ i t X; 

END DO 
DO L=OjM 
0L=D*L 
DP=D*(LH> 

CKL«JZjIZ)=(DL+YD[i(.IZ,rz>)**2 

P(LiJZ»IZ)stOP-rPP<JZ»lZ?»*»2 
S < L • JZ • I Z > - « 0 L+yPP< JZ«IZ>)** 2 
T ( L • JZ ■ I Z ' s ( DP- Y DD < JZ < I Z ) ■ » *2 
END DO 
END DO 
END DO 

THE PRIMARY COMPUTATION BEGINS HERE' 

THE FIRST STEP IS TO GET THE PHI-X VALUES FROM A *!.* GE 
SOURCE PANELS AT 

DO KX=-NX+1»NX 
XX=X(-NX)-X(KX) 

XS=XX**2 

FOR TIME SAVING WE SMEAR THE SOURCE IMAGES BEYOND » 
PAIRS IN BOTH HORIZONTAL AND VERTICAL DIRECTIONS. 

E=D/A* NZINZYPIYBETA 

EEE=1/E 

DM=(MF1)*D 

DMS=DM«2 

HM=(M+0.75)*2*A 

HM2=(M+0.25)»2*A 

HMS-HMYY2 

HM2S=HM2**2 

GS-EEE/3ETA*!ATAN<-XX*SORTl *S+IMS*HMS>/tRi*Dfi> 

ATAN< -XXISQRTK XSPDitS+HM2S » .’(H«?*DM i » ) » 

THESE ARE THE PHI-X VALUES FROM THE SMEARED SOURCES 

DO IZ=I f NY+NZ 

DO JZ=1 ?NY2 ! JUST THE FLOOR SOURCES 
DO I=-M,M 
P1=PT(I.JZ»IZHXS 


_ ■ J -_V lit 

QUALM 
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P2=PJ<I»JZiIZ)+XS 
DO L=0.« 

TT1=P1+Q<L» JZ* IZ ) 
TT1=TT1IS0RT<TT1> 
TT2-P2+R<L* JZ» IZ ' 
TTMT2<S<JRT< TT2> 
TT3-P1 +S<L»JZ»I2> 
TT3-TT3IS0RT ( TT3 > 
TT4«P2+T(L* IZ.IZ' 
m=lT4*S«RT(T'4' 


<JQsA<A/<;.’<P!*RETA*NIl»7 ) 
TTT=(1/TT1+1/TT'.!+’ TT?^ 1 'TT4)tQQ 


GO(JZ»IZ»KX)=-XX»!Tr .1? « IZrKX) 

END DO 
END DO 

THE TERN DYDZ APPEARS BELCW TC ACCOUNT FOR THE DIFFER- 
ENCE IN PANEL WIDTHS OF FLOOR AND WALL » IF ANY. 

GO(JZ»IZ»KX)=(60(JZ»IZ>KX)+6S)<IiYOZ 
END DO 

DO JZ S NY2+1 »NZ2 'LONER HALF OF THE wall SOURCES 

DO I=-NfN 
P1*PT<I»JZ>IZ)4XS 
P2=PJ<I»JZ»IZHXS 
DO L-OiN 

tti*pito(l*jz*iz) OEiCn-^r ; 

TTl=TTl*SQRT(TTt ) OE POOR Q 

TT2-P2FR(L.JZ.IZ) 

TT2-TT2*S0RT<TT2) 

TT3=P1+S<L»JZ«1Z> 

TT3=TT3<SQRT(TT3> 

TT4-P2+T <L*JZ*IZ> 

TT4= <T4<S0RT<TT4J 

TTT*<l/TTl+l/TT2+l/rT34l/TT4><00 

G0(JZfIZfKX)=-XXtTTT/BETA+60< JZ»IZ»KX) 

END DO 
END DO 

GO<JZ»IZ»KX)=60UZ>IZ»KX>M?S 
END DO 
END DO 
END DO 


NOTE THAT THE FUNCTION •><) 13 i*E» DELON TO CREATE 
Gi DERIVED IN THE REPORT. USfc OT THE REPEATING GO SAVES TINE 
AND NENORY SPACE WHEN A LARGE NUnBER OF PANELS ARE USED. 

THIS COHPLETES THE NOP\ FOR ALL M yiURC'S I MATED AT JX 
EQUAL TO -NX IN THE LONER HALF TUNNEL 

OPEN< UNIT =1»NA«E=' ilEASVX.DAT .STATUS*' OLD' : 

READdrO (<HEASVX(IZ«IX)« IZ=:<NY+NZ)*IX»-NX.N<) 


CLOSE U) 

0PEN(UNIT=2>NAME='MEASVN.DAT'»STATUS='GLB') r . r 

REAIK2f») < ( ME AS VN ( IZ f IX ) » IZ=1 f NY+NZ ) > IX --NX > NX i *’ 

CLOSE ( 2 ) Of - 


DO IX=-NX?NX 
DO IZ=1.NY+NZ 
IZ=NY+NZ+1-IZ 
DO JX=-NX»IX-1 
KX=-NX-( JX-IX) 
DO J7M,NY2+NZ2 


VX<tZ.IX)=VX<IZ>IXHMEASVN<JZ.JX>*GOiJZ,IZ*KX> 

END DO 

DO JZ=(NY2+NZ2+l«'NYrNZ 
KZ=NY+MZ+1-J/ 

VXaZ.IX)=VX(!<MX}+DEASi/N(JZ)JX)*GO<KZ.LZfKX) 

END DO 
END DO 

DO JX=(IX+1).NX 
LX=-HX+ JX-IX 
DO JZ=1»NY2+NZ2 

VX(IZ»IX)=MX<IZ»IX)-«EA'?UNtJZ*JX)IGO(JZ?!X»LX) 

END DC 

DO JZ=NZ2+NY2+1.NZ+NY 
KZ=NY+NZ+1-JZ 

VX<IZfIXi=VX(IZ>IX)+MEASVN(JZ*JX»*GO<KZ>LZ.KX> 

END DO 
END DO 

DO JX=(IX+1).NX 
LX=-NX+ JX-IX 
DO JZ=1*NY2+NZ2 

VX(IZ*IX)=VN<IZ*IX)-MEA3VN<JZ>JZ)IGO<JZ*IZ»LX) 

END DO 

DO JZ=NY2+NZ2+1 ? NY+NZ 
K2=NY+NZ+1-JZ 

VX<IZ»IX)=VN(IZ>IX)-M£ASVN(JZ»JX)*GO<KZ»iZ»lX) 

END DO 
END DO 

THIS IS THE PHI-X OF THE MEASURED NORMAL VELOCITIES 
DISTRIBUTION* NOW SUBTRACT THIS FROM THE MEASVX VALUES 
FOR INPUT INTO THE FOLLOWING PHINO COMPUTATION. 

VX(IZfIX)=MEASVX(IZ»IX)-VX(IZ>IX) 

END DO 
END DO 

NOW USING THE COMPUTED VX GET WALL-FREE NORMAL VELOCITIES 

THIS PROGRAM COMPUTES THE RESIDUAL NORMAL VELOCITIES 
AT THE TUNNEL WALLS GIVEN THE INPUT FUNCTION VX PROM THE 
PRECEDING WORK. THE OUTPUT CALLED PHINWALL.DAT IS 
THEN USED TO COMPUTE THE INTERFERENCE VELOCITIES 
AT THE POSITION OF THE MODEL USING PROGRAM 'PH|XZM\ 

IT ALSO USES THE INPUT FILE CALLED NEASVN THAT 
IS THE MEASURED WALL SLOPE PLUS BOUNDARY LAYER SLOPE. 

IT DEVELOPS THE PANEL EQUATIONS FOR THE NORMAL VELOCITIES 
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C PRODUCED BY A HORSESHOE VORTEX LYING AT THE CENTER OF A PANEL 

EE=2*A/NZ/BETA 'CIRCULATION FOR PANEL OF UNIT v'X 

SA=A/2/NZ 

SFA=A/2/NZ*DYDZ 

DO IZ=1»NY2 ' FIELD POINT ON THE FLOOR 

DO JZ=lfNY2 ' VORTEX ALSO ON THE FLOOR 

YDP(JZ»IZ)*Y(IZ)-Y(.IZ)«FA 
YDHUZ»IZ)=Y<IZ)-Y(JZ)-SFA 
YSP(JZ»IZ)*Y(IZ)FY(JZ)+SFA 

YSH(JZ.IZ)=T(IZ)+Y(JZi-SFA 

AO< JZ» IZ)=EEK l/YDPI JZf IZ I -1 /Y(W< JZt IZH 
1 l/YSF(JZ*IZ)-l/rSH< ! CASE ONE 

END DO 

YD(IZ)*Y(IZ)-0/2 

YP(IZ)*Y(IZ)«'.'2 

DO JZ=NY2+lfNY2+NZ2 • VORTEX CN LOWER HALF OF WALL 
ZDP(JZ)=A/2-Z(JZ)+SA 
ZD«(JZ)=A/2-Z(JZ)-SA 
ZYl(JZfIZ)=YD(IZ>*YD(IZ>+ZOmJZ)tZD«<JZT 
ZY2<JZ»IZ)*YD(IZ»>YD(IZ)+ZDf< JZ)tZDP(JZ) 
ZY3(JZ*!Z)=YP(IZ)*YP(IZ)«DM(JZ'«ZDNUZi 
ZY4(JZfIZ)*YP<IZ>tYPUZ>+ZDP<JZ)«ZBP<J2> 

A0(JZ»IZ>*EE*(YDaZ)/ZYI(JZ»IZ)-YD<IZ)/ZY2<JZfIZ:- 
1 YP<IZ>/ZY3<JZfIZWP<IZ)/ZY4<JZfIZ» 'CASE TWO 

END DO 
END DO 

DO IZ=NY2+1»NY2+NZ ! FIELD POINT ON THE WALL 
ZD(IZ)=Z<IZ)-A/2 

DO JZ=lfNY2 ! VORTEX ON THE FLOOR 

YDP(JZfIZ) s Y(IZ)-Y(JZH$F 
YDN(JZ»IZ)*Y(IZ)-Y(JZ)“SF 
YSP< JZ» I2)=Y( IZ)+Y< JZ)+SF 
YSN(JZfIZ) s Y(IZ)+Y(JZ)-SF 
YZl(JZiIZ) s ZD(IZ)tZD(IZ)+YDP(JZ»IZ)tYDP(JZ»IZ) 
YZ2(JZ*IZ)*ZD<IZ)*ZD(IZ)+YDN(JZ,IZ)»YDH<JZfIZ) 
YZ3(JZ»IZ)sZD(IZ)lZD(IZ)+YSP!JZ»IZ)*YSPUZ»IZ) 
YZ4(JZ*IZ) s ZD(IZ)tZD(IZ>+YSH<JZ>IZ)tYSN(JZ»IZ) 

AO(JZfIZ) s -EEtZD(IZ)*(l/YZl' ,JZ»IZ)-1/YZ2(JZ.IZ)+ 

1 1/YZ3(JZ*IZ)-1/YZ4(JZ»IZ> > < THREE 

END DO 

DO JZ=NY2+1» NY2+NZ2 ! VORTEX ALSO ON LOWER HALF OF WALL 
ZDP2(JZ*IZ)=Z(IZ)-Z(JZ)+SA 
ZDN2(JZfIZ)*Z(IZ)-Z(JZ)-SA 
ZDP3 ( JZf IZ)*ZDP2< JZ» IZ) tZDP2(JZ» I Z ) +0*0 
ZDN3 ( JZ f I Z ) =ZDH2 (JZ t IZ ) *Z0H2 ( JZ f IZ >+0*0 

A0(JZ»IZ)=EE<(l/ZDP2(JZ»IZ)-rZDN2(JZ.IZ"-EE«(ZDP2'J2.IZ)/ 
1 <ZDP2(JZtIZ)<ZDP2(JZfIZ)+0<0>-ZDH2< JZf IZ)/<ZDN2<JZ*IZ)* 

1 ZDH2(JZfIZ>+0*D)> • FOUR 


original page is 
OF POOR qualitx 


DO IZ-NY2+NZ+1 jNY+NZ OJE POOR 

DO JZ=1*NY2 

YDP(JZ»IZ)=Y<IZ)-Y(JZHSFA 

YDH<JZ.IZi=YfIZ)-Y(JZ)-SFA 

YSP(JZ»IZ)=Y(IZ)+Y(JZ)tSFA 

YSH( JZ»IZ>=f(lZ)+Y(JZ>-SFA 

YD1 < JZ» IZ V = fDP • JZ» I Z ) *YDP < JZ . I Z > +A*A 

YD2 ( JZ- 1 Z • = Y DM JZ » IZ WON < JZ » IZ) +A»A 

YD3iJZ> IZ > = rS?i JZi IZ.UYSP' JZ. FZHAtA 

YD4(JZ»IZ‘=fW:j:.!Z)»rSH(J7i!Z)«*A 


AO ( JZ , I Z > =-EE* ■ . Y DP ( JZ . ■ Z i YD H JZ .. I Z ! - 
YDH(JZ.:Z; ?02(JZ.lZ!t 
YSPUZ.Sr..' Yf.’f 17 * ! Z ' - 
YSHUZ.IZ W fD4( J7.IZ ■ » 


END DO 




YD6(IZ)*Y(IZ)-D/2 

YS6(IZ)=Y(IZ)+D/2 

DO JZ=NY2+1»NY2+NZ2 ! VORTEX ON T HF LOWER HALF WALL 
ZDP ( JZ ) =- A/2-Z ( JZ ) +SA 
ZDN(JZ)=-A/2-Z(JZ)-SA 

YR1(JZ»IZ)=Y06(IZ)*YD6(IZ)+ZDP<JZ)*ZDP<JZ> 

YR2(JZ»IZ)=YD6<IZUYD6(IZ)+ZDH(JZ)*ZDMJZ> 

YR3(JZ>IZ)=Y36(!Z)*YS6(IZ)+ZDP(JZ)>ZDP(JZ; 

YR4(JZ»IZ)=YS6(IZ;»YS<SUZ)+ZDN<JZ)*ZDmJZ) 


AO(JZiIZ)=EEt(YD4(lZ)/YRl(JZ*IZ)- 
YD6(IZ)/YR2( JZ.IZ)- 
YS4(IZ)/YR3(JZ»IZ)F 
YS4(IZ)/YR4(JZ»IZ)> 


END DO 
END DO 

THIS COMPLETES FILLING AO( JZ*IZ).NOU SO ON TO AX(JZ»IZ»KX) 

DO KX=-NX+l»NX 
XX=X(KX)-X(-NX) 

XS=XX*XX 


DO IZ=1»NY2 'FLOOR 
DO JZ=1«NY2 ! FLOOR 

ONE' 

R1=XSFYDP( JZ. IZ)>YDP(JZ>IZ) 

R1=SQRT(R1> 

R2=XS+YDH-,JZ»IZ)*YDH(JZ.IZ) 

R2=S0RT<R2) 

R3=XS+v5P(.MZ)*YSF<JZ»IZ) 

R3=S0RT '■ R3 > 

R4=XS+YS«> JZ.IZXYSmJZ.IZ) 

R4=SQRT<R4) 

AX<JZ»IZiKX)=EEXX>»(l/F;»( YDF< JZ.IZX1/XSH/YDP(JZ»IZ>)- 
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1 l/R2*(YHh(JZ.!Z)»l/XM/YDf1(JZ»IZ))+l/R3*(YSP(JZ»IZ)* 

1 1/XS+1/YSP(JZ.<IZ • i-t .'M*'Y5rff 'Z.IZUl/XSH/YSMJZflZ))) 

END DO 

DO JZ=1*WZ2 ! WALL VORTEX 

R1=XS«Y1(JZ»IZ) 

R2=XS+ZY2(JZ»IZ) 

F3=XS)ZY3UZfIZ> 

R4«XS4ZY4(JZ»IZ) 


AX(JZ»IZ»KX)=EEm*(l/Rl*YWIZ)/ZmjZ»IZ)-:.-P2Mi IZ). 
1 ZY2 ( JZ f I Z > •t/R3lYf > (IZ)/ZY3(JZ*IZ)41 <P4HY'- ll>; 

1 ZY4(JZ>IZ>! 

TWO 


END DO 
END DO 


1 


DO IZ-NY2+1 »NY2+NZ 
DO JZ S 1»NY2 

R1=XS+YZ1(JZ»IZ) 

Rl-SQRT(Rl) 

R2=*S+YZ2< JZjIZ) 

R2*SQRT(R2) 

R3*XS+YZ3(JZfIZ) 

R3=SGRT<R3) 

R4-XS+YZ4(JZ«IZ) 

R4=S0RT:R4) 

AX(JZ»IZ'KX>=-PE*XXiZ0(IZ)t<l/Rl/YZl(JZfIZ)-l/R2/ 

YZ2(JZ.IZ:M'R3/YZ3(JZiIZ)-1/R4/YZ4(JZ»IZ)) 


gplgdiaj: fag® m 

XTAI.TTYi 


THREE! 

END DO 

DO JZ=NY2+1»NY2+NZ2 

R1=XS+ZDP2( JZf IZ)*ZDP2 (£•!!■ 

Rl=SORT(Rl) 

R2=XS+ZDH2( JZr IZ) <ZDH2 (JZ?IZ> 

R2=SQRT(R2) 

R3=XS+ZDP3(JZf IZ) 

R3=SQRT(R3) 

R4=XS+ZDM3< JZ»IZ) 

R4=SORT(R4) 

AX(JZ»IZ>KX>=EE*XXtU/Rl*ZDP2(JZ»IZ)*0'X3+l/ZDP2UZ»IZ))- 
1 1/R24ZDN2 < JZ » I Z 1 1 < 1 /XS+1/ZDM2 ( JZ • l Z ) >-l /R3IZDP2 ( JZ > I Z ) > 

1 (l/ZDP3(JZ«[Z)+l/TKS4D*Di)+l/R4IZDH2(JZ»IZ)t 

1 (l/ZDH3(JZrIZ>+l / (X£*D*D )> 

FOUR 

END DO 
END DO 


XDD=XS+D*D 

DO IZ=NY2+NZ+1»NY+NZ 
DO JZ=1»NY2 


R1=XS4YD1(JZ»IZ) 
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Rl-SGRT(Rl) 
R2=XS+YD2( JZ»IZ) 
R2=SGRT(R2) 


R3=>:3+YH3< JZ.I2) 
R3=SGRT'R3) 
R4=XS+YD4( J7’IZ> 
R4=SQRT(R4) 


Oh 


X' 


AX(JZ»IZ»KX*=-EE*XX*<YIiF<JZ»IZ)/Rl*' 1 XCO+i r?i- JZ,\Z))- 
1 YDM ( JZ » IZ ) /RZ H i /X0D+ 1 /YD2 ( JZ . IZ ) ) + 1 Sf : JZ . 1 2 i R3* 

1 < 1 ■' XDD+l/YSF ■: JZ - [ 7 : > - < 3H< jZ • I Z > / R4* ( 1 /MB* 

1 F1VF' 

ENG DO 

DO jZ=NY2+1 • NY2+NZ2 

R1=XS+YR1(JZ»IZ> 

R1=SQRT(R1) 

R2=XStYR2<JZ?IZ) 

R2=SQRT(R2) 

R3=XS+YR3< JZfIZ) 

R3=SQRT <R3 > 

R4=XS+YR4 ( JZ » I Z > 

R4=SQRT(R4) 


AX(JZ>IZ»KX)=E£fXXI(YD6( IZM(1/R1/YRI(JZ» (Z •-l'R2' v R2(JZ>IZ) )+ 
YS6(IZ)»i;/R3/YR3(JZ»IZ)-l/R4/YR4(JZ.IZt)i 


C I SIX 

END DO 
END DO 
END DO 

C ARRAY AX IS NOU FILLED! THE ARRAY A DISCUSSED IN THE REFORT 

C IS SYNTHESIZED FROH AO AND AX TO SAVE SPACE WHEN A LARGE 

C NUMBER OF PANELS ARE USED, 


DO IX=-NX»NX 
DO IZ=1»NYFNZ 
LZ=NY+NZ+1-IZ 
DO JX=-NX,IX 
KX=-NX-JXPIX 
DO JZ=1,NY2+NZ2 

PHINO(IZ»IX)=PHINO(IZ,IX)m(JZ»JX)*(AO(JZ,IZ.' + 
1 AXUZ<IZ,KX) ) 

END DO 

DO JZ=NY2+NZ2+1 .NY+NZ 
NZ=NYPNZ+l-r 

PHINO(IZ»IX)=PHInO; IZ i I<)+VXi JZ< JX)l(A0(KZiLZ)4 
1 AX(KZ,LZ»KX) ) 

END DO 
END DO 

DO JX=IX+1,NX 


LX=-NX+JX-IX 
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DO JZ=1»NY2+NZ2 


1 


PHINO(IZ» IX ) =f'HINO ( I Z » I X > +VX < JZ . JX > * < JZ r IZ > - 

AX( JZfIZjLX) > 

END DO 


DO JZ=NY2tNZ 2+1 *Nf +NZ 


KZ*WY+MZ+1-J2 

PHiNOUZ*iX3=phiNO(iz.ix)m(jz»Jx»»(.;o:xz.i.:^ 

1 AXOCZfLZfUu 

END DO 
END DO 
END DO 
END DO 

C THE PHI NO ARRAY IS NOW FILLED' 

C REMEMBER THAT PHINO IS THE NALL-FREE NORMAL VELOCITY 

C CALCULATED FROM THE MEASURED WALL FHI-X AND PHI-N 

C VALUES' HE ASVX AND MEASVN). 


DO IX=-NX»NX 
DO IZ=1»NY+NZ 

VN(IZ»IX)=PHINO(IZ*IX)-«EASVN(IZ»IX) 

C THIS IS THE RESIDUAL EFFECT OF THE HALL PRESENCE. IF THESE 

C TWO VALUES UERE EQUAL WE UOULD HAVE PERFECT WALL ADAPTION 

END DO 
END DO 

OPEN(UNIT=3»NAME s 'PHINUALL.DAT'»STATUS='NEW' > 

URITE<3»*) <<VN(I,J)fM,NY+NZ)»J=-NX»NX> 

CLOSE! 3) 


END 


c-^ 
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